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ERRATA  AND  ADDENDA 


Vol.  I 


p.  2a,  1.  3*:  Read  "Vlllat"  for  "Vlllet." 

p.  23,  1.  2:  The  end  of  the  line  should 

read,  "where  y»o  Is," 

p.  88,  formula 

(13:40):  In  the  first  term  read  z-c 

for  z-o. 


Vol.  II  , 

* 

p,  147,  1.  13:  Replace  the  sentence  beginning 
"Weinstein  by  "Com¬ 

pleteness  follows  from  known 
criteria  [see,  e.g.,N.  Levinson, 
Gap  and  density  theorems,  Amer. 
Math.  Soc.  Colloq ,  Publ.  no.  27, 
New  York,  1940,  ch.  I],  However, 
both  orthogonality  and  com¬ 
pleteness  are  consequences  of  the 
general  theory  of  Sturm-Llouvllle 
systems. 

d.  148.  formula 

(16.5):  The  limits  of  the  integrals  with 

respect  to  y  should  be  -m  and 
0  instead  of  0  and  oo  . 

'•in. 

p.  150,  1.  7*:  Read  23a  for  22. 

p.  154,  1.  12:  Read  e  fop  a 

p.  156,  1.  7:  The  line  should  begin  "in  terms 

of  the  wave  length  and  the  con¬ 
formal  mapping . " 

p.  156,  1.  2*:  The  line  should  begin  "Ursell, 
and  also  Mamyanskli  [1954] , 
also  ..... 

p.  157,  1.  9:  Add  to  this  paragraph  the  follow¬ 

ing  sentence:  "Bartholomeusz 
[1958]  treats  the  long-wave 
approximation  for  reflection  at 
a  step  in  the  bottom." 
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. t . "if"  ’ 


. . . . . . . . . . . 


Pc  157,  1,  L3,  14: 


p.  159,  I,  12: 
p.  160,  1.  3: 


p.  161,  1.  5-7: 


p<  161,  l-i  10: 


p.  167,  1.  7*: 
p.  175,  1.  13: 


p.  179,  1.  6: 


p.  ibO.  formula 
(17.30): 


p,  180,  1„  2*: 

p.  181,  1.  4; 
p.  191,  1.  7: 


p.  204,  1.  2*: 
p.  205,  1.  8*- 
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These  two  11ms  should  read, 
"using  the  completeness  of 
the  functions  [cf.  (16.1)]." 


For  "value"  read  "value  (see 
Kreisel  11949]  or  Meyer  [1955])." 


Insert  the  following  sentence 
between  this  formula  and  line  4: 
"Other  general  considerations  will 
be  found  in  Biesel  and  LaM£haut£ 
[1955] 


Read  "Dean  [1945  1946]"  for "Dean 

[1945]."  In  lines  5  and  6  delete 
the  clause  "the  ....  Tr/in."  In 
line  7  read  "[1957]"  for "[1958]." 


Add  the  following  to  this  para¬ 
graph:  "Dean  ana  John  alio  treat 
barriers  inclined  at  an  angle  if /In . 
Levine  and  Rodemich  [1958]  solve 
the  vertical -barrier  problem  by 
several  methods,  including  the 
cited  ones ,  and  then  apply  one  of 
them  to  the  problem  of  waves  in¬ 
cident  upon  two  parallel  vertical 
barriers . " 


Read  ic: 


Complete  the  sentence  by  "when  the 
body  is  completely  submerged."' 


Between  "again"  and  "by"  insert 
"bv  Mac dona 1  i  [ 1896] ,  Pocklington 
[1$21]  and  In  line  7  read  "All" 

for  "Both." 


On  the  right  side  read  J ur  f  (})> 

Read  "Lehman"  for  "Lehmann." 

Read  (17.32)  for  (17,12). 


For  "enough"  read  "enough  and  the 

body  la  submerged," 


Read  "Dcnn"  for  "Dunn." 


Replace  the  third  sentence  of 
this  paragraph  by  the  following: 
"The  solutions  (18.33)  for  Y  sr 
7r/2(2.«a  l)  have  also  been  given 
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p.  208,  1.  1*: 

p#  2X8 1  1  *  7 : 


p.  216,  i.  8*: 


p.  216,  1.  1*: 


p.  218,  formula 
(19.11): 


p.  225,  1.  7: 


by  Macdonald  [1896].  Ac  Chase 
critical  angles  the  solution  (18.33) 
does  not  vanish  as  x -*<*>,  Macdonald 
apparently  discarded  the  other 
solutions  as  being  of  little  interest, 
not  'being  sensible  at  a  distance 
from  the  edge.'  Roseau  [1958]  has 
recently  carried  through  a  systematic 
investigation  of  edge  waves,  in* 
eluding  ones  with  singular  be¬ 
havior  at  the  edge." 

For  "discovered"  read  "states." 

Insert  the  following  sentence: 

"Alblas  [1958]  treats  similar 
three-dimensional  problem  in  which 
the  motion  is  periodic  along  the 
length  of  the  strip." 

For  "approximately"  read  "approxi¬ 
mately;  by  improving  the  approxi¬ 
mation,  Levine  [1958]  has 
clarified  certain  anomalous  re¬ 
sults  of  Kochin  for  an  oscillating 
horizontal  plate." 

Add  the  following  sentence  to  the 
paragraph:  "A  general  survey  of 
methods  of  generating  waves  in  the 
laboratory,  including  some  account 
of  theoretical  results,  may  be 
found  in  a  recent  paper  by  Bieiel 
and  Suquet  [1951,  1952]." 

The  first  integrand  should  read 

5*  ~  ^  ( T)]  . 

In  the  space  following  "breadth" 
insert  "b." 


p.  231,  1.  4:  Delete  the  sentence  beginning  on 

this  line. 


p.  239.  formula 
(19.63): 

p.  241 ,  1 .  6 : 


Insert  "«"  before  "Re," 

i  I'l 

For  "forces"  read  "coefficients," 


p,  242,  1.  11:  Replace  "by  MacCamy  [1954]"  by  "by 

Barakat  [1958]  (in  an  earlier  in¬ 
vestigation  by  MacCamy  [1954]  the 
multipole  teimis  In  the  poten^l^l 
for  the  diffracted  wave  were  omitted). 


’mm  . . 


Page  4 


p.  243,  1.  4: 


Replace  Che  rest  of  the  paragraph, 
starting  with  "MacCam^,"  by^  the  0 


following:  "Baraket  shows  that  cf 
can  be  found  as  a  series  in  functions 
of  the  form  (13.21),  with  ^rOand 
account  taken  of  certain  symmetries, 
and  functions  of  the  form  (13.20)  with 
Jtr  ~  Dand 
Let 
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Then  ^  may  be  expressed  as  follows 

;Z  £v^ 


1  Z.  o.  C-,U  •* 


°0 

z  z: 

1t=  1  ^-O 
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where  the  complex  coefficients  Aim 
AVj  C«  are  to  be  determined 

from  (19.75).  No  numerical  computa¬ 
tions  seem  to  be  available." 


I: 


247,  1.  1,  3,  4, 


Read  for  "  evJ»/V'  In 

line  3  read  "  >i"  for  "  >0-" 


p.  248,  1.  1*: 


To  the  sentence  terminating  on  this 
line  add  "(see  also  Lamb  [1934]).'> 


p.  249,  1.  14: 


In  the  parentheses  Insert  "Kochin 
[1937  ab,  1938  c] before  "Long." 
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p.  250,  1.  7: 

p.  250,  foxmula 
(20.18): 


Insert  "(*o,-W,^.  )"  after  "point." 
Read  "  *-=?*  +  A  "  and  "  1^**  },  " 
for  "f"  .nd  "  i  "  respectively. 


p.  257,  1.  8-9: 
p.  257,  1.  11: 
p.  266,  1.  5: 

p.  268,  1.  6: 

p.  274,  1.  4*: 


Read  "singularities"  for  "slngularlza 
tion." 

Read  "a  single -valued  complex"  for 
"the  complex." 

Read  "[1951,  1958]"  for  "[1951]." 

To  the  end  of  the  paragraph  add  the 
following:  "Cummins  [1954]  finds  the 
additional  effect  of  a  train  of  waves 
on  the  surface." 

Read  "Froude"  for  "Fronde." 


p.  275,  1.  7*,  8*:  The  statement  should  be  weakened  as 

follows:  "'mis  is  confirmed  only 
approximately  in  photographs;  diffi¬ 
culty.... 


p.  276,  1.  6,  12,  Read  "  "  for  In 

1*:  formula  (20.73)  the  lower  limit  of  in¬ 

tegration  is  "-q.." 


p.  288,  1.  9:  Add  the  following  to  tils  paragraph: 

"More  recently  the  problem  has  been 
considered  again  by  authors  unaware  of 
the  earlier  work.  Squire  [1957]  has 
analyzed  a  gliding  flat  plate  by  a 
method  similar  to  that  used  by 
Sretenskii  and  Maxuo.  Certain  integrals 
involved  in  this  mer.hod  have  been 
tabulated  by  Miller  [1957]. 

Cumberbatch  [1958]  has  used  a  method 
similar  to  Sedov's.  Both  authors  add 
new  results  to  the  earlier  work." 


p.  301,  1.  3:  Add  "  f>=  +  j  ,"  to  this  line. 

p.  311,  1.  8:  Add  the  following  to  this  paragraph: 

"Also,  in  the  neighborhood  if  <v  = 
the  expressions  in  (21.36)  are 
inaccurate;  in  this  region  *t  may  be 
expressed  in  terms  of  Airy  functions 
(see  Ursell  [19591)."  In  lines  5  and 
8  read  "Fig,  la"  for"Fig.  1." 


p.  315,  1.  5*:  Replace  this  line  by  "analogues  of 

(13.49)  and  (13.53)  for  pressure  dis¬ 
tributions  may  be  found  in  equations 
(22.48)  and  (22.49)  or  in."  In  line 
4*  read  "these"  for  "this." 


Vol.  Ill 


p.  329, 

1,  7: 

Read  "final"  for  "find." 

p,  32 2, 

1.  11*: 

Read  -for 

p.  341, 

1.  6: 

Insert  "(see  alto  Wurtele  [1955])" 
following  "[1953]." 

p.  367, 

1,  3: 

Add  the  following  to  the  paragraph: 

"A  problem  somewhat  related  to  those  of 
this  section  is  the  motion  of  a  freely 
floating  body  In  a  fixed  bounded  basin 
(there  is  now  no  dissipation  of  energy 
as  in  the  problem  treated  at  the  end  of 
section  22 jO.  This  problem  has  been 
deal  to  with  by  Perzhnyanko  [1956]  and 
Moiseev  [1958]." 

p.  378, 

1.  2*: 

Read  "  ryabi  "  for  "  ryaby 

p.  422, 

1.  11: 

Read  "approximation"  for  "approxima¬ 
tions." 

p.  422, 

1.  1*: 

Read  [1919a]  for  [1918]. 

p.  428. 
(27.21): 

formula 

Read  T '  'P0*'  3t'*v"fW, 

p.  430, 
(27.29): 

formula 

(  1 

Read  77  for  7-7  ■ 

1  u  16 

p,  433, 

1.  3*: 

Insert  "and  by  De  [1955]  to  the  fifth 
order"  following  "third  order." 

(27.43) : 

formula 

Read  U  for  ,  to  conform  with 

the  notation  elsewhere. 
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A.  Introduction 

The  various  problems  of  fluid  motion  treated  in  this 
article  have  in  common  the  property  that  the  fluid  i"  subject 
to  a  gravitational  force.  In  addition,  in  almost  all  cases 
they  also  have  in  common  the  presence  cf  surfaces  separating 
two  fluids  of  different  densities  or,  if  only  one  fluid  is 
present,  of  so-called  free  surfaces.  However,  not  all  fluid 
flows  falling  into  this  category  are  treated  here;  tidal 
motion  is  treated  in  vol.  48  in  the  article  by  A,  Defant. 

The  observed  properties  of  ocean  waves  and  their  generation 
by  wind  are  treated  in  the  article  by  H  U.  Roll,  also  in 
vol.  48.  Closely  related  problems  concerning  flows  with  free 

surfaces  are  treated  in.  the  article  by  b.  Gilbarg . in  this 

volume , 

The  subject  of  water  waves  engaged  many  of  the  mathema¬ 
ticians  and  mathematical  physicists  of  the  last  century. 

Moreover,  the  last  several  years  have  brought  a  renewed 
interest  in  the  theory  of  water  waves  In  addition  to  this 
extensive  literature  on  theoretical  aspects  of  the  subject, 
there  have  also  been  many  experimental  investigations, 
usually  carried  out  by  hydraulic  engineers.  Hydraulic 
engineers  have  also  produced  an  extensive  literature,  both 
theoretical  and  experimental,  on  open  channel  flow,  flow  over 
weirs  and  through  sluice-gates,  etc.;  included  is  a  consider¬ 
able  literature  on  numerical  and  graphical  methods  of  solving 
the  equations  involved.  Oceanographers  have  produced  their 

|L 

own  literature,  usually  emphasizing  different  aspects  of  the 
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subject.  The  theory  of  ship  waves  has  produced  its  own 
literature. 

All  this  material  is  pertinent  to  this  article.  Clearly 
some  selection  has  to  be  made.  We  have  followed  roughly  the 
following  rules:  Fundamental  results  axe  derived  in  full. 

The  treatments  of  various  special  problems  are  selected  so  as 
to  exemplify  particular  methods,  other  methods  being  mentioned 
only  by  literature  citation.  Experimental  results  are  not 
usually  reproduced,  but  references  are  given.  Numerical 
methods  of  solving  equations  are  not  treated  at  all.  The  more 
special  problems  of  hydraulic  engineering  are  also  not  treated. 
Geophysical  aspects  which  are  omitted  have  already  been 
men  tioned . 

Several  excellent  expositions  of  the  theory  of  waves  or 
of  various  parts  of  it  already  exist.  We  mention  the  follow¬ 
ing*:  Lamb  (1932,  chs.  VIII  (pp.  250-362)  and  IX  (pp.  363- 
475)] ;  Basset  [1888,  ch.  XVII  (pp.  144-187)];  Wien  [1900, 
ch.  V  (pp.  166-224)];  Kochin,  Kibel',  and  Roze  [1948,  ch.  8 
(pp.  394-526)];  Milne-Thomson  [1956,  ch.  XIV  (pp.  374-A31)]; 
Airy  [1845];  Bouasse  [1924];  Auerbach  [1931];  Thorade  x9'" 
Sretenskii  [1936];  Khristianovich  [1938];  Keulegan  [l^j./; 
Eckart  [1951];  and  Stoker  [1957].  The  last  cited  book  by 
Stoker  gives  an  up-to-date  account  of  much  of  the  fundamental 
theory.  For  observation  of  waves  of  many  kinds,  Cornish 
[1910,  1934]  and  Mlche  [1954]  should  be  consulted.  Shuleikin 
[1953,  part  3  (pp.  213-292)]  contains  a  general  discussion  of 

^References  are  collected  at  the  end  and  identified  in  the 
text  by  author  and  date. 


topics  of  Interest  In  oceanography .  Russell  end  Macmillan 
[1952]  give  a  rather  nontechnical  discussion  of  ocean  waves. 

A  volume  published  by  the  Society  of  Naval  Architects  of  Japan 
(ZSsen  Kydkei)  contains  expositoty  papers  on  various  aspects 
of  water-wave  theory  related  to  ships  (see  Maruo  [1957], 

Jlnnaka  [1957],  Nlshlyama  [1957],  Bessho  [±.957],  and  Inul 
[1957]). 

For  extensive  bibliographies  one  should  consult  Thorade 
[1931,  pp.  195-211];  Sretenskll  [1936,  pp.  294-303];  Kampe  de 
Feriet  [1932,  pp.  225-229];  and  Stoker  [1957,  pp.  545-560]. 
Sretenskll  [1950,  1951]  In  a  survey  of  the  accomplishments  of 
the  USSR  during  the  years  1917-1' 47  has  given  a  rather  complete 
bibliography  of  Russian  papers  during  those  years.  Takao  Inul 
[1954]  has  Included  a  valuable  bibliography  of  Japanese  papers 
in  a  survey  of  Japanese  contributions  to  the  theory  of  ship 
waves.  An  Interesting  early  history  of  the  subject  may  be 
found  in  a  paper  by  St.-Venant  and  Flamant  [1887],  The 
treatise  by  the  Weber  brothers  [1825]  Is  still  of  interest 
for  its  content,  and  especially  for  Its  many  references  to 
and  summaries  of  the  early  papers  on  water  waves.  The  section 
on  waves  of  the  article  on  hydrodynamics  by  Love  [1914],  as 
modified  by  Appell,  Beghln  and  VI lie.,  In  the  Encyclopedic  des 
sciences  mathfaratiques  gives  brief  indications  of  the  contents 
of  many  of  the  papers  published  up  to  about  1912. 


B.  Mathematical  Formulation 


l  •  Coorj^nate,  sy  a  tiTOf^jR^^^y^nj^png 

In  the  mathematical  description  of  waves  one  may,  as  in 
fluid  mechanics  in  general,  describe  the  motion  by  describing 
either  the  paths  of  individual  fluid  particles  ("Lagrangian" 
description)  or  the  velocity  (and  acceleration)  field  in  the 
region  occupied  by  fluid  at  a  given  moment  ("Eulerian"  de¬ 
scription)  .  Generally,  but  not  always,  the  Eulerian  descrip¬ 
tion  will  ba  used. 

Rectangular  coordinates  may  be  used  conveniently  for  al¬ 
most  all  problems.  The  y-axis  will  be  taken  directed  oppo¬ 
sitely  to  the  force  of  gravity,  "he  x-axis  and  z-axls  so  as 
to  form,  a  rlgnt-handed  system  (i.e  ,  If  the  y-axis  is  toward 
the  top  of  the  page  and  the  x-axis  is  toward  the  right,  the 
z-axls  will  point  toward  the  reader).  This  is  a  somewhat 
unconventional  choice  for  the  z-axis,  but  has  the  obvious  ad¬ 
vantage  that  in  two-dimensional  problems  one  can  delete  z- 


dependent  terms  from  the  equations,  have  conventional  (x,y) 
coordinates,  and  set  z«x+iy  without  ambiguity  when  complex- 
variable  methods  are  convenient . 

It  seems  hardly  worth  while  to  try  to  formulate  rules 
concerning  when  a  moving  coordinate  system  is  preferable  to  a 
fixed  one.  However,  use  of  a  moving  coordinate  system  t«  clearly 
convenient  in  those  cases  where  it  allows  one  to  formulate  a 
problem  in  a  time -independent  manner, 

The  following  well-established  convention  with  rugard  to 
use  of  certain  letters  will  be  adhered  to.  The  compr/nents  of 
the  velocity  vector  ^  will  be  denoted  by  u,v,w  the  pressure  by 
f  and  tha  dansity  by  f .  Ihe  coefficient  of  viscoaJty  of  the 


MfliiHNM 


1  S 


•  1  ||H|INjH||ll||i||t||||i  R||  I  |)|l||iP||)||H|lt|||t.|i|  ||||  Nlpp|||l||l||[|||h  | j?|j 

. .  . . . . . . 


|HI# 


4 


fluid  will  be  dtnottd  byp,  tht  coefficient  of  kln«Mtlc 
vt  *<  osl  ty  ,  »  ft  ,  by  v.  The  acctUratloti  resulting  frc*s 
gravity  (•  denoted  by  %. 

In  the  ftulerten  formulation  one  seeks  v  y  end  f  as 
functions  of  •  )  )  t  l.e,  at  eny  Instant  t  one  seeks  o 
vector  function  and  two  scalar  functions  defined  on  the 
region  accept*  *  t*j  fluid  at  that  instant.  In  the  Lag* 
rattan  system  me  focuses  attention  on  the  trajectories 
of  individual  particles  la  the  fluid  if  a  are  the 
ce**rd1  nates  af  a  particle  at  time  t  *©  ,  than  one  seeks 
the  position  i  «  ( i  t  ;  •  k  ‘  )  *  <  *  f  *f  «bf  a 

paint  at  a  later  tlae  t.  One  soy  peas  free  one  seats* 
to  the  other  by  aeons  of  the  epatlma 


with  »  V  -  (  ’  ‘  of  *■*  -  as  Initial  condl tlons 


2 .  Equation*  of  motion 

Derivations  of  the  fundamental  equations  describing  fluid 
motion  are  available  in  many  places  (e.g.,  vol.  8  of  this  Ency 
clopedla).  The  equations  are  reproduced  here  for  convenience 
of  reference. 

The  equation  of  continuity  in  Eulerian  coordinates  is 


li  ,+  ZjH  4  =0 

I>  t  d  x  3  -j.  v  ^ 


(2.1) 


If  the  fluid  is  incompressible,  but  not  necessarily  homogeneous, 
Jj/Jt  mQ  (but  not.  necessarily  •  0)  and  eq.  (2.1)  becomes 


7)  U  7>  V-  9  !~r 


(2.2) 


In  Lagrangian  coordinates  this  may  be  written 
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For  an  incompressible  fluid  f  (*/,.« 

becomes 

T>r  i 
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f (q/€, C ,0)  and  (2,3) 


(2. A) 


The  dynamical  equations  take  different  forms  according  as 
>ne  does  or  does  not  try  to  take  account  of  viscosity.  The 
Havier "Stokes  equations  for  the  motion  of  an  incompressible 
viacoui  fluid,  when  the  only  externnl  fdtpt  is  that  of 


Irdwrty 


are . as  follows  in  Eulerian  coordinates; 
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If  viscosity  is  neglected,  the  last  two  terns  on  the  right  side 

of  the  equations  are  to  be  deleted  and  one  obtains  the  equations 
for  an  "ideal"  fluid: 
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In  Lagrangian  coordinates  the  latter  equations  become: 
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The  equations  of  two-dimensional  motion  result  if  one 
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deletes  all  terms  containing  z,  w,  and  c . 

The  motion  is  called  irrotational  if  it  satisfies  the 
additional  equations 


j  w-  _  \r , 
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or,  in  two-dimensional  motion, 
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(2.8*) 


r>  v  >•  '-i 

In  the  case  of  irrotational  motian  there  exists  a  potential 


function  i$(x,  y,  z,  t)  such  that 

'~*H‘  (2-9) 

It  is  a  classical  theorem  of  hydrodynamics  (cf.  Lamb  [1932, 
§§17,  33])  that,  if  the  motion  of  an  inviscid  fluid  with 
f  •  |>(p)  la  irrotational  at  any  instant,  it  is  so  thereafter. 
Ip  particular  a  potion  started  from  rest  is  irrotational. 
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If  the  equation  of  state,  the  following  Integral 

of  the  equations  of  motion  exists  for  irrotatlonal  motion: 
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where  P  •  (  f  ^  j>  and  ;\(t)  is  an  arbitrary  function  of  t. 

ft 


If  llie  fluid  Is  incompressible,  the  usual  case  in  this  article, 
j>  is  independent  of  p  and  the  integral  becomes; 


Jt 


,  ££■=  Alt). 


t  (-N  ■)  +  va  ■  ■  r 

In  this  case  one  obtains  also  from  (2.2)  and  (2.9) 


(2. 10’) 
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Even  if  the  motion  is  not  irrotatlonal,  th'- *e  still  exists  an 
integral  like  (2.10)  if  tie  motion  is  steady,  the  so-called 
Bernoulli  integral: 


'  (uSv-W)  +  Jj'fr  +  ?  =  C  (2.10") 

Here  C  is  constant  along  a  single  streamline: 

,  ig-.-r  4}, 

but  may  vary  from  one  streamline  to  another. 

There  will  be  occasion  in  the  following  to  treat  problems 
In  moving  coordinate  systems.  Let  Opr  -  y  be  e  fixed  coordinu  » 
system  and  be  a  system  moving  with  respect  to  Ox^but 

without  rotation.  Let  V;  be  the  vector ■jjfSO, the  velocity  of 
a  particle  referred  to  Oryybe  y*  and  to  be  fr.  Then 

V^ys .  We  shall  generally  want  either  to  describe  the  abso¬ 
lute  motion  v  with  respect  to  the  moving  coordinate  system 
OXj  j.  or  the  relative  motion  yr  with  respect  to  this  coordinate 
system.'  In  either  case  the  continuity  equation  remains  the 
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The  dynamical  equations  for  an  ideal  fluid  for  the  absolute 
motion  described  in  the  moving  coordinate  system  are: 
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The  dynamical  equations  for  the  relative  motion  are: 
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’*  11  ^ 

Let  us  suppose  that  the  motion  is  irrotatlonal  and  let 

be  the  velocity  potential  for  the  absolute  motion  in  the  fixed 


coordinate  system.  Let 

r*  , 


C^fy  +  ju.i^y  [  +  J  l*,  t)  ' 

Then  $>  is  the  velocity  potential  for  the  absolute  motion  in 


the  moving  coordinate  system: 

?i.\*  ,  §4=^,  '^1 

ds  »•*  DI 

The  integral  (2.10)  becomes; 


r  w. 
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then  $  is  the  velocity  potential  for  the  relative  motion: 
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and  the  integral  (2.10)  may  be  written: 
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The  more  general  equations  when  the  system  d*y^  is  also  ro¬ 
tating  will  not  be  necessary  for  this  article. 


3 .  Boundary  conditions,  at  an  Interface 

Let  us  now  suppose  that  we  are  given  two  immiscible  fluids 
with  a  cannon  boundary  surface,  S(t)  .  The  one  fluid,  with 
density  f  and  viscosity  ,  will  occupy  region  TV, (t)  ;  the 
other,  with  density  fv  and  viscosity the  region  I<W. 

Let  describe  the  surface  we  assvme 

F/  Fj 1  >G  (where  F„  -c>r/ax,  etc.). 

The  first  condition  which  the  surface  S(t)  must  satisfy 
is  a  kinematic  one.  As  the  surface  moves,  the  velocity  of  a 
point  (a,  jj^)  on  the  surface  in  the  direction  of  the  normal  to 
the  surface  is  given  by  -  j\j  r  \  .  Here  one  takes 

the  normal  in  the  direction  (fA)  F, , F^j  •  A  particle  of  fluid 
at  the  same  point  of  the  surface  at  that  instant  will  have  a 
velocity  component  in  the  direction  of  the  surface  normal  given 
by  For  S/f)  to  be  a  bounding  surface  means, 

v  F,  +  F>j\  F jv 

of  course,  that  there  can  be  no  transfer  of  matter  across  the 
surface.  Consequently  the  following  equation  must  be  satis¬ 
fied: 


u  F„  *  -  Ft  .  (3.1) 

-here  we  heve  weed  the  aseumptlon  f,’  >0  In  dropping 

the  denominators.  If  one  defines  the  "material  derivative" 
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by  Che  equation 

ff  =  “F-  +vFJ  ’•u’F!r+  F<  < 

Chen  (3.1)  is  Che  same  as 

£f  -  o  (3. 

Dt 

This  condiCion  must  be  satisfied  by  any  bounding  surface, 
whether  an  interface  or  a  rigid  boundary*. 

There  are  further  dynamical  conditions  to  be  satisfied 
at  an  Interface.  Let  us  first  consider  the  general  case  of 
viscous  fluids  with  surface  tension  at  the  interface.  The 
following  assumptions  are  made: 

1)  The  effect  of  surface  tension  as  one  passes  through 
the  interface  is  to  produce  a  discontinuity  in  the 
normal  stress  proportional  to  the  mean  curvature  of 
the  boundary  surface. 

2)  For  yiscous  fluids  the  tangential  stress  must  be  con¬ 
tinuous  as  one  passes  through  the  interface. 

3)  For  viscous  fluids  the  tangential  component  of  the 
velocity  must  be  continuous  as  one  passes  through 
the  interface. 

In  order  to  formulate  these  statements  in  mathematical 

language,  we  introduce  the  following  notation.  Let  %( * ,  £,J.) 

be  some  function  defined  in  both  1?,  and  T?v  and  let  f*o 

be  a  point  of  the  Interface  S.  Assuming  that  the  following 

limit  exists,  we  shall  write  ^  \)  a8 

*For  further  discussion  of  this  condition  see  Truesdell,  Bull. 
Tech.  Univ.  Istanbul  l  (1950),  no.  1,  71-78  (1951);  Lichten¬ 
stein,  Grundlagen  der  Hydromechanik ,  Springer,  Berlin,  1929, 
pp.  159-170,  234  ff. 


Let  the  component®  of  the  strata  tensor  be  denoted  by 


I 


C» 


5  *  VT:  *  aj- 

r<»  ra 

Consider  an  element  of  area  of  the  surface  5  at  a  point  (V,  ]  ],) 
of  S.  Let  the  unit  normal  vector  to  S  at  (*•],£}  be 
Then  the  surface  element  will  have  associated  with  It  the 
stress  vector  with  components: 

i  *  5*1  w'  +  n  j  ^  k  ^  ■*  ^  n  »  ^  ^  **  *  j  . 

Let  1?,  and  7?v  be  the  principal  radii  of  curvature  of  5  at 
.  Then  statements  1}  and  2)  are  combined  in  the  one 
equation 


k,  2 

K*  k  - 


»  T  if  T<' ,  ~ J  ■-+■  ?,"}  }, 

1?,*"  )  ^ 


r-n  -♦  cra  ^  a 


(3.2) 


[<$ „  (  »  T  -*  1? v*' )  n , 

where  T  is  a  constant  of  proportionality  depending  upon  the 

two  fluids  (and  their  temperatures,  but  this  will  not  be  con¬ 
sidered  here).  T  io  called  the  coefficient  of  surface  tension*. 

The  kinematic  condition  imposed  in  (3,1)  is  clearly  e- 
qulvalent  to  continuity  of  the  normal  component  of  the  velocity 
as  one  passes  through  S*  Consequently  the  condition  3)  above 


may  be  combined  with  this  to  give 


U  v 


v.  »  v. 


Uj  s S.  UA 


(3.3) 


*For  an  air-water  interface  T*7l,8  dynes /cm.  at  20  C,  for  mercury  - 
air  interface  T»4f.r  dynes /cm.  at  20*  C,  for  a  mercury-water  interface 
T„4I2  dynes/cm.*  for  benzene-air  T  *  IS.  9  dynes/cm.  at  2(fCJfor 
liquid  helium-helium  vapor  T*ai4 dynes/cm.  at  -270°C. 
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In  the  linearized  theory  of  viscosity  the  stress  tensor  Is 
given  by 

't't'*}***) 

-  j*[v:  *u/_)  (3.. 4) 
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f  -  Mu  „  -*■  V'^  +  Uy ')  -  If 

/'a (ur_  p  - 

If  the  fluid  is  incompressible,  the  terms  containing  the  second  vis¬ 
cosity  coefficient  h  vanish. 

The  geometric  quantity  R"1  +  R^1  la  given  tv  the  formula' 

R,+i-  ^  (3.5) 

The  sign  Is  so  selected  that.  If  It  is  positive,  the  direction  of 
Increase  of  the  normal  component  of  the  stress  vector  at  the  inter¬ 
fact  Is  In  the  direction 

f  r  C 

•  x  i  >»-  r  i 


.  (3-6) 


'  /  FA  Fa\  1  /fM\  f^F1 

In  the  case  of  a  surface  given  by  J?  f equation  (^.5)  becomes 

-jr*  4  = 

I?,  f?v  (i  •*  7,'*  7J1- 

In  the  case  of  two-dimensional  motion  this  simplifies  further  t  »  the 
well-known  formula 


i  ** 


(l (3.5") 
If  one  now  substitutes  (3.4),  (3.5),  and  (3.6)  in  (3.2),  one 

obtains  the  general  boundary  condition  at  the  interface. 

^"See ,  e  .g . ,  A  *  Duschek  and  V.  Mayer,  Lehrbuch  der  Differentialgeo- 
metrle,  Bd.  I.  Teubner,  Leipzig-Berlin,  1930,  pp.  150-152. 


The  result  Is  unwieldy  in  its  general  fora*. 

If  Che  interface  is  given  by  3.=  7^  ,  the  boundary  condition 


becomes 


ones 

tpi'f,  -iur««]7„  *-  = 

(fj  - 1  Crfv-  ♦*>)]?,  -  nr  *  l  (3  7) 

W?v  -/  ir H-  -,>]?„  *  2 1  r~',; f?;‘)  *, 

with  R,  -*RV  given  by  (3.5').  Here  fluidj  is  the  lower  and 
fluidj  the  upper  fluid.  For  two-dimensional  motion  the  equa¬ 


tions  take  the  following  fora:  ,, 

[p]  y(M)  -  ( 1  [  p  ft*)  -  [ r  (u.  - 1; ) j]  =  t  - — , 


[hi  +  j 1  NVL^]?7'»-  ^>^ij  -  r 


(3.8) 


rr:^  ■ 


(  !  -*  •J'M'" ; 

One  may  also  write  this  condition  in  terms  of  the  components 
of  the  stress  vector  normal  and  tangential  to  the  interface: 


[f]  -  ?  .  .  -r^ _ 

"  0-rvjVt) 


7  *  I1'  A ) 


2f^-‘*"^<)b/  +  [>K  f 7'v~  j) 

(,N^l|l/l  ~C' 

If  surface  tension  is  to  be  neglected,  one  obtains  the 


(3. 8’) 


resulting  boundary  condition  by  setting  T - 0  in  the  various 
equations  above.  In  this  case  equation  (3.2)  simply  states 
the  continuity  of  the  stress  vector  as  one  passes  through  the 
interface. 

If  viscosity  is  neglected,  but  not  necessarily  surface 
tension,  the  condition  on  the  stress  vector  becomes  simply 

[f]  =  t(  0.9) 

where,  of  course,  the  mean  curvature  is  still  given  by  (3.5), 


*In  tensor  notation  the  condition  is  somewhat  more  perspicuous: 

in  ft  r  <  n  tt  ,  *  >  C  r  cp?fr  A  -  I  _ 


Jfj  ,1, 


?fpl  -Du.  *.1  a*.  |  -l*lu.  tu.  .■)]'/  — — Eli -  —  t  “Fii-Rrf.ij 

r  JJ  (f.9  W  *  (F,»  QV*  ’Jitt 

a..t  (*,.A  x,j ,  J.,).  (ui1,ui,,i4)),Hi,vt.)..j  F..-S9F/3*/  .  We  have 

reirainea  from  using  tensor  notation  because  its  particular 
advantages  cannot  in  general  be  exploited  here. 


lib;  iii  Wiiillilli!  liliti  j[ ; [1  |,:i 
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The  other  boundary  condition  (3.3)  changes  more  drastically 
upon  neglecting  viscosity:  Condition  3)  stating  the  contin¬ 
uity  of  the  tangential  component  of  velocity  is  abandoned. 
The  continuity  of  the  normal  component,  i.e.  (3.1),  is  still 
retained,  of  course.  Condition  2)  concerning  the  tangential 
stress  is  satisfied  vacuously  for  an  inviscld  fluid. 

So  far  we  have  considered  the  boundary  condition  at  an 


interface  between  two  fluids.  If  the  second  fluid  Is  absent, 
the  boundary  surface  for  the  first  fluid  is  called  a  "free 
surface".  Usually  the  pressure  above  a  free  surface  Is  as¬ 
sumed  to  be  some  given  function,  say  ft  (*,),£, f)  >  position 
and  time;  in  most  cases  it  is  taken  to  be  a  constant,  either 
an  assumed  atmospheric  pressure  or  zero  The  boundary  con¬ 


ditions  concerning  the  stress  vector  at  a  free  surface  are 
slight  modifications  of  those  for  an  interface,  and  can  be 
obtained  by  setting  The  result  is  again  somewhat 

unwieldy  in  its  complete  form*  For  an  incompressible  fluid 

( p  -  p)  f„  +(v°r7  4  ^  '  [r<'  '?v  jF'  ■  (3.10) 

(P  -r>)F;*rj  Ti  f<[ 

(f  -  p)  fi  *  p  f  *  ■*})  ft  *  f  u-r- 1  !*  i  '  lr\  V  34  T f  14  ^  C  ^  • 

Here  we  have  written  n  for  ,  and  f*  W  p,  ;  p,  p,^  ate  to 

be  evaluated  at  ( A'  “]/  G 

The  case  with  which  we  shall  be  chiefly  concerned  is 
that  of  an  inviscid  fluid  without  surface  tension  and  with 


*In  tensor  notation  it  may  be  written: 


? 


i/i 


l*.-f1(tF,rj  -F.r  F,rF.t; 
-  "F  rjr  "c"  \>/v 


Here  we  have  written  P  for  JV  and 
for)i,h  .  All  variable  quantities  in  the  braces  are, 
of  course,  to  be  evaluated  at  the  free  surface  /=*0  . 
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p(4  j  p,  ,  a  constant.  In  this  case  the  boundary  condition 
reduces  to  the  single  equation 

|>  *  A-  j  J  0  *  P°  (3.11) 

on  -  O  .  If  the  motion  Is  Irrotatlonal  and  Incom¬ 

pressible,  one  may  determine  f>  explicitly  from  (2.10*)  so 
that  (3.10)  becomes 

£  ,  «  -  (3.11') 

Tt  L 

to  be  satisfied  on  F  ( *  ;,j  t)  ;0 

In  the  case  of  steady  motion  of  an  incompressible  fluid, 
the  Bernoulli  integral  (2.10”)  still  exists  even  if  the  motion 
is  rotational.  Consequently,  in  certain  two-dimensional  prob¬ 
lems  of  steady  motion  in  which  the  free  surface  is  a  streamline 
one  continues  to  have  a  boundary  condition  like  (3.10''): 

J  l-k  ♦  ”  :  •»  *  fm/f  ~  C  (3.11") 

to  be  satisfied  on  r'*  -  r  j 

■4 

4 .  Boundgr^  c  cndit  ion  son  rial  d  sur  fac  e  s 

Let  the  equation  of  the  rigid  surface  be  given  by  the 
equation  6  J,f)  -  0  .  Then  in  the  case  of  au  inviscid 
fluid  the  condition  to  be  satisfied  on  Ci  -0  is  the  same  as 
the  kinematic  condition  (3.1): 

^  A  ■+  V*  G  y  • 1  W  Gj  r'  ^  ,  (4.1) 

i.e.,  the  component  of  velocity  of  the  fluid  normal  to  the 
surface  must  equal  the  velocity  of  the  rigid  surface  in  the 
direction  of  its  normal. 

If  the  fluid  is  viscous,  it  must  stick  to  a  solid  boun¬ 
dary  and  move  with  it  without  slippage  An  equation  of  the 
form  C  (*,  j-,  ?  ,t)~  o  is  not  ruitable  for  for.  ala  ting  this 
statement  in  equations  (e  .g. ,  -A  jv*  /  -  does  not  distinguish 
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between  a  rotating  and  a  stationary  Bpher©),,  Let  eh®  surface 
be  given  in  parametric  coordinates  by:  x  r  X(  -  ^  -  Y(VaA 

j.sZ(r,j,N  ,  where  a  given  point  on  the  surface  corresponds 
to  a  given  pair  of  values  (r,\)  .  Then  the  condition  for  vis¬ 
cous  fluids  may  be  written: 

_  -7>X  .  .  aZ 

A~  Tt  '  ^  '  Tt  St  '  <4. 2) 

If  a  solid  boundary  penetrates  the  free  surface  (or  an 
Interface)  of  a  viscous  fluid,  there  will  be  some  difference 
In  treatment  of  the  boundary  condition  according  as  the  fluid 
wets  the  surface  or  not.  In  the  case  of  mercury  sloshing  In 
a  clean  glass  basin,  the  fluid  pulls  free  of  the  surface  as 
It  moves  up  and  down,  whereas  water  in  the  same  basin  will 
continue  to  adhere  to  any  part  of  the  walls  already  wetted. 
Furthermore,  if  surface  tension  is  taken  Into  account,  the 
angle  of  contact  of  the  free  surface  with  the  solid  surface 
will  enter  into  the  boundary  condition;  In  the  first  case 
mentioned  above  the  angle  may  vary  according  as  the  liquid 
Is  rising  or  falling  along  the  wall*.  Although  attempts  to 
prove  very  general  existence  theorems  for  fluid  motion  would 
presumably  take  such  complications  into  account,  they  are 
usually  neglected  in  most  solutions  of  special  problems, 
there  being  indeed  little  choice  in  the  matter, 

5.  Other  types  of  boundary  surfaces 

Geophysical  problems  sometimes  suggest  situations  in 
which  there  is  an  interface  between  a  fluid  and  an  elastic 
medium.  This  may  occur,  for  example,  in  the  study  of  the 

*See,  e.g„,  R.  S.  Burdon,  Surface  tension  and  the  spreading 
of  liquids,  Cambridge,  1949,  pp.  76-82, 
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effect  of  ocean  wave#  on  the  ocean  floor,  as  In  Lonf  tggini" 
11950]  theory  of  alcroaelama .  Other  possibilities  «  suggested 
by  wav#  motion  o«  a  body  of  water  covered  'with  ar  Ice  sheet  or 
at  an  Interface  between  two  fluids  separated  by  an  elastic 
membrane  or  plate.  In  one  series  of  Investigations  the  Ice 
sheet  has  been  assumed  broken  Into  pieces  small  with  respect 
to  the  prevalent  wave  lengths.  In  this  case  only  the  density 
of  the  Ice  layer  enters  into  the  modified  boundary  condition 

«i 

(see  Peters  [1950],  Keller  and  Goldstein  [1953],  Keller  and 
Welt*  [1953],  Shapiro  and  Simpson  [1953]).  Waves  In  a  thin 
plate  over  an  Infinitely  deep  fluid  have  been  considered 
briefly  by  Landau  and  Lifshits  [1953,  pp.  762-3],  but  with 
neglect  of  gravity  Greenhill  [1387,  p.68;  1916]  included  gravity. 

The  kinematic  boundary  condition  (3.1)  must  always  hold. 

The  dynamical  conditions  will  depend  upon  the  nature  of  the 
assumptions.  The  matter  will  not  be  further  considered  here. 


C.  Preliminary  Remarks  and  Developments 
6.  Classification  of  problems 

Moat  of  the  theory  of  water  waves  Is  concerned  either 
with  elucidating  son*  general  aspects  of  wave  motion  or 
with  predicting  the  behavior  of  waves  in  the  presence  of 
some  special  configuration  of  Interest  to  oceanographers, 
hydraulic  engineers,  or  ship  signers  Unfortunately, 
even  some  of  the  apparently  simplest  problems  have  proved 
too  difficult  to  solve  in  their  most  complete  formulation. 
Approximations  have  been  necessary,  and  in  many  cases  the 
problems  which  have  been  solved  are  those  which  could  be 
solved  by  the  approximate  methods  in  use.  An  examination 
of  the  theory  also  shows  that  many  of  the  concepts  and 
definitions  are  almost  Inextricably  bound  up  with  these 
methods  of  approximation,  following  rather  than  preced¬ 
ing  the  making  of  th?  approximation. 

The  nature  of  the  approximations  used  in  treating  a 
particular  problem  provides  a  natural  way  of  classifying 
It.  First  there  are  the  assumptions  concerning  the  prop¬ 
erties  of  the  fluid:  viscous  or  invlscid,  compressible 
or  incompressible,  surface  tension  or  not.  Althovgh  as¬ 
suming  the  fluid  to  be  invlscid,  Incompressible,  and  with¬ 
out  surface  tension  simplifies  the  equations,  they  are 
still  not  easily  manageable,  even  for  the  simplest  kinds 
of  problems.  Other  approximations  of  a  different  nature 
are  required.  These  are  in  a  sense  mathematical  approx¬ 
imations.  Their  physical  significance  is  not  in  restrict¬ 
ing  the  nature  of  the  fluid  but  in  restricting  the  charac- 
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tar  of  the  waves  and  the  boundaty  configuration.  The  kind 
of  mathematical  approximation  used  provides  another  means 
of  classifying  problems,  and  Is  the  principal  one  which 
will  be  used  in  this  article.  There  are  two  principal 
methods  of  approximation,  explained  below  in  section  10, 
the  infinitesimal-wave  approximation  and  the  shallow-water 
approximation.  Thus,  the  development  of  these  two  approx¬ 
imate  theories  and  of  the  exact  theory  constitutes  the 
bulk  of  this  article. 

7.  Progressive  waves  and  wave  velocity- -Standing  waves 

It  will  be  convenient  to  call  any  motion  of  a  fluid 
in  a  gravitational  field  with  a  free  surface  or  an  inter¬ 
face  a  wave  motion . 

If  the  velocity  components,  pressure,  and  free  sur¬ 
face  or  interface  may  be  expressed  in  the  form 

respectively,  then  the  wave  motion  will  be  said  to  be  a 
progressive  wave  traveling  in  the  direction  C*.  In  this 
esse  a  change  to  a  moving  coordinate  system  with  t 

1>z )  i  r*duces  the  motion  to  steady  notion  with  respect 

to  the  moving  coordinate  system  With  rear:  .ct  to  the  fixed 
coordinate  3ystem  the  profile  of  the  free  surface  or  inter¬ 
face  is  being  transported  vlthout  change  of  form  in  the 
direction  J *•  with  velocity  c.  It  might  seem  reasonable 
therefore  to  call  C  the  velocity  of  propagation  of  the 
progressive  wave. 

However,  Stokes  [1849;  or  1880,  pp.  202  f f . ]  has  pointed 
out  that  the  velocity  of  propagation  of  the  profile  of  the 
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free  surface  does  not  by  Itself  give  a  useful  definition 
of  wave  velocity.  Let  the  fluid  be  invlscid,  either  in¬ 
finitely  deep  or  with  a  horizontal  bottom,  and  unlimited 
otherwise.  Now  let  the  whole  fluid  in  the  progressive 
wave  described  above  be  transported  with  velocity  C  (pos¬ 
itive  or  negative)  in  the  direction  Ox.  Then  the  motion 
will  still  be  consistent  with  the  laws  of  fluid  mechanics, 
the  various  parts  of  the  fluid  will  move  the  same  rela¬ 
tively  to  each  other,  but  the  velocity  of  propagation  of 
the  profile  will  be  arbitrary,  depending  upon  the  choice 
of  C .  What  is  required  for  a  useful  definition  of  wave 
velocity  is  the  velocity  of  propagation  of  the  profile 
with  respect  to  a  coordinate  system  fixed  in  some  sense 
in  the  fluid 

In  the  case  of  an  infinitely  deep  fluid,  if  the  axes 
may  t  chosen  so  that  as  y -*--*•  the  velocity  relative  to 
these  axes  vanishes,  then  one  may  reasonably  measure  the 
profile  velocity  with  reapect  to  these  If  the  motion  far 
ahead  or  far  behind  the  disturbance  approaches  a  uniform 
velocity  (possibly  zero),  then  axes  moving  with  the  fluid 
with  this  velocity  may  be  used  When  the  disturbance  does 
not  behave  thus  (as  in  the  case  of  periodic  waves)  and  when 
the  depth  Is  finite,  there  is  no  longer  an  o’»vious  way  to 
select  a  set  of  reference  axes. 

In  order  to  put  the  problem  somewhat  differently,  let 
us  assume  that  the  wave  motion  is  given  as  a  steady  motion 
with  velocity  field  W*,  j)  and  free  surface  -}  . 

We  wish  to  find  a  moving  coordinate  system  *'=■■* y’-'frj 
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•o  that  In  some  sense  the  relative  motion  vanishes  on  the 
average.  We  now  have  the  free  surface  given  by  y '  -  t) 

and  the  relative  velocity  by  vr y'j  *  -  *.  i  . 

7  /  t 

How  Is  to  be  chosen?  Stokes  made  two  suggestions.  One 

Is  to  define  It  by  the  equation 

/,  r- 

j  r  .  .  i 


J:z.u,  rri  i>'} 

a  .  ^ 


f 


(  -  rx 


r  ! 

)  6*  J  [  -  w,J  ly  -  £  / 


(7.1) 


where  }-s~k  is  the  equation  for  the  bottom.  In  case  the 
motion  is  periodic,  with  period  >i  ,  the  defining  equa¬ 
tion  may  be  written 


>  iU) 

\  <Lx  J  /><</}) -w-le 


4 


(7.2) 


o  -  t 

If  one  notes  that  the  mean  depth  is  given  by 


t. 


> 


f  'll*)  +UJ  lx  flr  d  y 

then  one  sees  that,  with  h'  as  mean  depth, 

mc  h  -  a 


(7.3) 


where  Q  is  the  average  discharge  rate  per  unit  width.  ua  is 
thus  defined  so  that  the  average  discharge  rate  with  respect 
to  the  coordinate  system  is  zero.  is  usually  de¬ 

noted  by  C'  . 

Stokes'  other  suggestion  was  to  define  u,  by 


JU,  /)  Ji  =0 

T-*  oo  T  J  c 


(7.4) 
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If  u  Is  pet  iodic  in  x  with  period  )t  ,  one  may  write 


w « 


,  4'*> 

~~  \  i* 

>  ) j1 


(7.5) 


In  either  case,  for  the  definition  to  be  useful  ue  must  be 
independent  of  /  and  y  .  If  u  Is  bounded,  it  follows  eas¬ 
ily  that  6Vi/Im'  =  C  for  both  cases.  If  the  motion 
is  irrotational,  u;  -  ^  and  it  follows  again  that 
'bUt/dfy  *  0  if  v  is  bounded.  Wave  velocity  defined  in 
this  manner  is  usually  denoted  by  C,  For  the  two  special 
cases  considered  earlier,  the  two  definitions  coincide. 

The  definition  of  wave  velocity  In  cases  where  the 
motion  cannot  be  reduced  to  a  steady  motion  is  no  longer 
s tralght forward .  In  many  cas^.s  of  interest,  the  asymp¬ 
totic  behavior  of  the  motion  for  large  positive  or  nega¬ 
tive  /  allows  one  to  define  a  wave  velocity  in  a  manner 
similar  to  that  above.  In  more  complicated  wave  motions 
one  may  simply  follow  the  motion  of  some  special  phase 
of  the  profile,  say  a  crest.  This  provides,  for  example, 
a  definition  of  phase  velocity  for  a  cylindrical  wave. 

A  general  definition  of  standing  wave  is  somewhat 
more  awkward  to  formulate  than  that  for  a  progressive 
wave.  For  the  case  of  a  plane  wave,  the  free  surface 
y  -  K  M  must  be  periodic  in  each  of  Jr  <i«d  t  ,  with 
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wave  length  X  end  period  c  ,  say.  In  addition,  the  curves 
in  the  (*,  i )  -  plane  represented  by  ,  '.'lien- 

the  undisturbed  surface,  must,  consist  of  two  sets  of 
curves  oscillating  about  the  lines  *  •  ~  ^  and 

t  ■=  -  *  r  y  .  For  progressive  waves  the  curves 

'I  La, 1 1  -  -  consist  of  a  single  set  of  straight  lines,  all 
parallel  to  x-  ei  ~  C  .  The  prototype  for  the  standing 
wave  Is  the  surface  defined  by,  say,  y-  i  •  wuit-rf/v. 

However,  as  shown  by  both  Penney  and  Price  [1932b]  and 

by  Sekerzh-Zenkovich  [19471  ,  neither  set  of  curves 

consists  of 

straight  lines,  or  even  fixed  curves,  for  standing 
waves  of  finite  amplitude. 

There  remains  the  problem  of  establishing  that  pro¬ 
gressive  and  standing  waves  exist  under  suitable  boundary 
conditions.  For  the  exact  boundary  conditions  for  a  per¬ 
fect  fluid,  the  existence  of  progressive  waves  was  first 
established  by  Levi-Civita  [1925]  and  Nekrasov  [1921,1922]. 
The  existence  of  standing  waves  satisfying  the  exact  boun¬ 
dary  conditions  is  apparently  an  open  question. 

8 .  Energy 

Let  T ft)  be  a  region  occupied  by  a  perfect  fluid  with 
a  boundary  bit'  represented  by 

F  ].  0  ' 

the  representation  being  chosen  so  that  is 

in  the  direction  of  the  exterior  normal.  The  surface  l(t) 
moves  independently  of  the  motion  of  the  fluid.  It  is 
assumed  that  contains  no  singularities  of  V”  and 
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that  surface  tension  doe#  not  act  upon  the  surface  SfH 
at  any  tine.  The  energy  of  the  fluid  contained  In  T ft) 
Is  given  by 
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E  -  \\  \  l  * 


(8.1) 


For  irrotational  motion  of  an  lnviscld  incompreoaible 
fluid,  one  may  use  (2.10')  and  express  e  by 

c  --  SlU-f  -flfj  i- 

7  b) 

(Here  <f  has  been  redefined  so  that  MO  may  be  set  equal  to 
zero.)  One  may  now  compute  aC/n  t  by  using  the  general 


formula: 


j/j  V^j.yi  jt  ;;  !■* 

T  ft 


-  r  f  lr 

r  'it.  ^  r  v 


/  F*  ^  r  „  ■*  Cf 

j 


One  finds  [cf.  F.  John,  1949,  p.  19  ff.,  which  we  follow 
closely  here] : 

if  =  W  J f"" hini *• ir  +  ]  ir 

SI  i)  w» 

by  Green's  Theorem  and  the  equation  of  continuity.  Finally, 


c/tl 


/F,v,  trk  j  (8.2) 

;  * 


We  recall  that  -  ft  // T^*  F,%  is  the  velocity  of  $(t) 

a 

in  the  direction  of  the  exterior  normal.  Two  cases  are 
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of  special  Interest.  If  bW  Is  s  "physical"  boundary, 

*1, 

!•*•»  on*  moving  with  the  fluid,  then  the  first  summand 
vanishes  and  one  finds 


[ 


Aj 

if 


-  -  Si  p 

ZH> 


cf.  Lamb,  Hydrodynamics,  p.9, 
geometrical  '  boundary,  then 


lit  A<r 

7*  il 

eq.(5)).  If  S(h  Is  a  fixed 
r  o  an(j  one  gets 


(8.3) 


UK 

;ih 


d  T 


(8.4) 


If  one  considers  any  portion  of  $/f)  ,  then  the  in¬ 
tegral  of  (8.2)  taken  over  this  portion  and  with  a  minus 
sign  gives  the  rate  of  flow  of  energy  through  this  portion 
of  SJtl  .  in  case  a  par.  of  5/f)  is  a  physical  boundary 
which  is  fixed,  c'^/tn  *  0  and  the  flow  through  this  part 
is  zero.  The  same  conclusion  holds  for  any  portion  of  S*/f) 
that  is  a  free  surface,  for  then  p  =  0. 

If  one  has  a  progressive  wave  moving  to  the  right  with 
V  J,}  and  t*kes  S  as  a  region  in  the 
fixed  plane  X  =  x.  ,  then  the  rate  of  flow  of  energy 
through  S  in  the  positive  direction  is  given  by 

H  f  ‘  ^  *  o,  ( 

i.e.,  energy  always  flows  in  the  direction  of  the  wave. 

In  cases  where  one  Is  dealing  with  waves  generated  by 
moving  bodies,  it  is  frequently  possible  to  choose  the  re¬ 
gion  T  so  that  no  energy  is  lost  from  it,  .he  latter  being 
true  only  as  an  average  If  the  motion  is  periodic  in  time. 

As  an  example,  consider  a  body  moving  steadily  with  velocity 
C  In  the  X- direction  in  an  infinite  ocean  with  horizontal 


bottom.  In  addition  to  the  boundary  conditions  on  the  body, 
free  surface,  and  bottom,  we  assume  that  the  motion  vanishes 
(in  the  limit)  far  ahead  and  to  the  sides  of  the  body.  The 
surface  Sit)  may  the  i  be  chosen  as  a  plane  ft :  =0 

far  ahead,  another  plane  A/.-{*-cf)*r=-G  behind  the 
body,  planes  K  *  t  .  J.=  10.  on  either  side,  and  the 
bott9in  H  ,  the  wetted  surface  of  the  body  3  ,  and  the 
part  of  the  free  surface  F  included  between  the  body  and 
the  planes.  The  energy  within  this  region  is  clearly  con¬ 
stant,  and  one  easily  obtains,  with  cj>  ( f (*’ c  ^  ^  )♦ 

O  -  -  \\  ['  >7  7  -  tL  ¥,  ^  i  c[j  [J  f 

j  [  l  V  (  tf/  t  1  Jl  )■*  i  l  k]  *  ~  • 

N 

’Since  on 3  one  has  ^ /7>o  =-  e  ,  one 

find®  for  the  first  integral,  remembering  that  f\  points 
into  the  body, 

-  \\  f»  ~~  At  ~  ~  c  1;  p  ^  ^  r  , 

where  ^  is  the  force  on  the  body.  The  parts  of  the  second 
integral  over  M  vanish  as  u.-*  oo  and  similarly  for  the 

flr.t  summand  In  the  third  Integral.  The  terms  In  fjy 

give 


which,  as  a.  -*  so 


,  converges  to 

if 


One  obtains  finally 


2? 


•R  -  f  f  fi>  f  «£«■).«)] 

k  +  \  ii  i  iVfjWj.- 

Thl»  exact  formula  for  resistance  will  be  put  Into  a  dif¬ 
ferent  form  later  after  linearization  of  the  boundary  con¬ 
dition*.  Although  the  plane  a -<■*'*  C  may  be  taken  at 
any  distance  behind  the  body  without  destroying  the  val¬ 
idity  of  (0.6),  it  is  usually  convenient  to  take  it  so 
far  behind  that  asymptotic  expressions  for  cf  can  be  used. 

If  in  (8.1)  a  part  of  the  surface  5/ B  ,  say  S , /  +)  » 
is  an  Interface  with  another  fluid  with  surface  tension 
acting,  then  the  energy  is  given  by 

£  Z.  j]5  [i  f  tf <  j  dr  +  T  i(  ir  . 

Trf)  S|W 

Let  S  {¥)  be  bounded  by  the  curve  C(f)  given  paramet¬ 
rically  by  *}(*,♦),  J.KB  and  let 

fortnui®  analogous  to  (8.2)  is 

U  =  \\  +  +  $  P  Jr 


(8.6) 


(8.7) 


At 


Lsir  +T  1 


S,(H 


j 

On 


f.  S  fj. 

\  v.  J. 

3.  M 


u  <»•«> 


)T,3 


If  S, (f)  ia  a  free  surface,  then  the  boundary  condition 

l°°  ‘  p r 

where  b6  is  an  assumed  constant  pressure  implies  that*, 
there  is  no  flux  of  energy  through  S, . 

If  the  motion  is  two-dimensional,  with  S,  given  by 


[|  ,N||| Jrl-ll'l'  dll  III  f  ir||i!l|ll| 


4|.i|l.ll  IMIi.bw>|ilillilwMMWiiH..fl.MI.H>i|ll'l>'l'll  l"H  «  f|»  l44lP  kf  «d  1 1**)*  I 

liktl!  1 1  lillMMI  1HikniiiMriiM[i|l|:  ijif i ij 


. 


in 

-r  . 

vmt  e  .r.ilii,  r |t]  jmn 

wmHHillinHMjmil 

iiiiiiwwipw . i«i 


4,  tl  i  X  i  X%IV 


%  S 

f  '  0 


,  then  (8.7)  becomes 

fi  f  (**v*  ’/'■  +  fljj”  J«r  +  j  di*'i 
r,H  -  *'<» 


(8.9) 


.dA 


and  (8.8)  becomes 

if :  I  r4*‘'?-  +  7^f7Jj/l+ll’7F^i 

",tL 

'1  f  !>■*  T  4  )vt  j»  +  T  4=2^  •*  T/l”)^  »'wl 

If  5,  Is  a  free  surface,  the  Integral  over  $,  may  be 
dropped  by  suitably  redefining  |j  . 

9  a  |  - 

Expressions  for  rate  of  change  of  momentum  may  be  de¬ 
rived  which  are  analogous  to  those  for  rate  of  change  of 
energy.  With 


(8.10) 


ft 


\\\  f  * 


(9.i: 


and  otherwise  the  ..'me  notation  as  In  section  1 ,  one  finds 

*fi.  Kff  *  ?  + .  TTZfpf  1 J<r 


>  tri 


»  -  fr:  I Jr- 


(9.2; 


Here  the  first  line  of  (9.2)  Is  derived  by  a  direct  com¬ 
putation  of  ilft/Jt  with  V-  ~  grad  ^  ;  the  second  Is 


derived  analogously  to  (f  ’>.) ;  the  third  follows  directly 


ft1 . . . . ftt ;  ;SnhiN>>aHMHM’ 


by  use  of  (2.10').  Comparison  of  line  a  one  ar.d  three  gives 

the  known  relation  (Levi-Clvlta): 


Y\  z  ~  55  ( V'  T )  T  i,r • 

5  s 


(9.3) 


Note  tha,,  In  (9.2)  and  (9*3)  SM)  may  move  In  an  arbitrary 
manner  as  l°rg  as  the  region  TW  contains  no  singularities 
and  only  liuid.  If  the  boundary  is  physical,  the  terms  in 
square  brackets  vanish  in  (9.2);  if  the  boundary  is  fixed, 
then  Fj  • 

Let  $«  (t)  be  a  physical  boundary,  possibly  the  surface 
of  a  solid  body,  and  SW  a  closed  surface  containing  it. 
Applying  (9*2)  to  the  region  of  fluid  bounded  Jointly  by 
$0  and  S  ,  one  finds 


fo  =  55  (f  +  pj*)*  Jcr 

So 


-  -  55  pf ^5  ■+  y*j»  fHr  -+ JJ  f (iv'e-p*  y) J<r- 
s.  "  s' 


(9.4) 


Here  F0  is  the  hydrodynamic  force  on  Sa  and  does  not  in¬ 


clude  the  hydrostatic  force. 

If  singularities  are  allowed  in  the  region  occupied  by 
fluid,  they  may  be  enclosed  in  spheres  of  small  radius  and 
the  formula  (9-4)  applied  to  the  remaining  fluid,  with  S 
modified  to  include  the  spherical  surfaces.  If  the  singu¬ 
larities  are  isolated  sources  of  strengths  at  the  points 
d;  ,  then  by  shrinking  the  spheres  about  the  singularities 
in  a  customary  fashion  (cf.  Milne-Thomaon  [1956.  pp .448-450] ) , 
one  obtains  the  following  modification  of  (9.4): 


+yr<ivVr-sir)J3‘)  (9>5) 

$*>  -S 


where  nr;  is  -he  velocity  at  the  point  fit*  when  the  source 
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«t  that  point  la  removed.  Other  modi  fleet!  onus  may  be 
derived  for  other  type®  off  singularities . 

If  the  velocity  field  Is  such  that  t**V  ~»o  as 
r  -+ «e  for  some  t»0  ,  then  the  last 

Integral  in  (9.4)  or  (9.5)  will  vanish  as  S  is  expanded  to 
infinity,  provided  the  latter  can  be  done  without  destroying 
the  validity  of  the  formula.  In  the  case  of  a  body  moving 
In  a  fluid  with  a  free  surface,  one  cannot  expand  in  all 
directions  and  must  Include  the  contribution  of  the  last 
integral  over  the  free  surface.  However,  the  formulas  are 
still  useful  In  computing  the  force  on  an  obstacle  result¬ 
ing  from  waves. 

10.  Expansion  of  solutions  in  powers  of  a  parameter 

In  their  exact  form  even  the  simplest  problems  with 
surface  waves  are  difficult  to  solve.  If  one  neglects  vis¬ 
cosity  and  assumes  Irrotational  motion,  the  problem  Is  re¬ 
duced  to  flndlrg  solutions  of  Laplace's  equation,  which  is 
at  least  linear  In  the  unknown.  However,  the  problem  Is 
still  difficult  because  of  the  nonlinear  boundary  condition 
at  the  free  surface  or  interface.  This  lack  of  linearity 
deprives  one,  for  example,  of  the  mathematical  tool  of 
superposition  of  solutions  expansion  in  eigenfunctions 
or  use  of  Green's  functions  is  not  possible. 

In  order  to  be  able  to  treat  special  problems,  the 
equations  are  approximated  by  ones  which  are  more  tractable. 
The  two  principal  methods  of  approximation  may  each  be  treated 
as  a  perturbation  procedure.  As  was  mentioned  In  section  7, 
this  procedure  Is  not  concerned  with  the  assumptions  about 
the  nature  of  the  fluid,  for  example,  whether  or  not  vis- 


coaity  is  neglected,  but  rather  with  the  nature  of  the 
motion  and  its  generation.  An  adv/Mtage  In  using  the 
perturbation  procedure  Is  that  the  assumptions  about  the 
motion  ar*  displayed  in  such  a  way  that  it  is  clear  how 
to  obtain  approximations  of  higher  order.  The  me  t  h  o  d  ha s 
been  applied  to  water-wave  problems  by  Sekerzh-Zenkovlch 
[1947,  1951,  1952],  K.  Friedrichs  [1948] ,  Keller  [1948], 

F.  John  [1949],  Longue t -Higgins  [1953b],  Peters  and  Stoker 
[1957],  and  others.  As  used  here  the  method  is  purely 
formal,  the  nature  of  the  convergence  of  the  perturbation 
series,  whether  it  be  uniform,  pointwiae,  asymptotic  or 
what  not,  being  left  open.  However,  for  each  method  of 
approximation  it  is  possible  to  point  to  at  least  one  case 
in  which  convergence  has  been  proved:  for  the  infinitesi¬ 
mal-wave  approximation,  Levi-Civita's  [1925]  proof  of  the 
existence  of  a  periodic  wave  of  permanent  type;  and  for  the 
shallow-water  approximation,  Friedrichs  and  Hyers1  [1954] 
proof  of  the  existence  of  a  solitary  wave. 

To  a  certain  extent  the  two  methods  of  approximation 
have  different  aims.  The  infinitesi  al-wave  approximation 
fits  ini  a  *  aral  scheme  for  approximating  nonlinear  equa¬ 
tions  ana  nvndary  conditions  by  linear  ones  (see  Souriau 
[1952]  for  a  liscussion).  To  apply  it,  one  must  know  a 
particular  exact  solution  to  start  with.  In  addition,  one 
must  be  able  to  select  a  dimensionless  parameter  (or  param¬ 
eters)  ,  say  £  ,  which  helps  to  determine  the  exact  physical 
problem  and  la  such  that  the  solutions  to  the  exact  problems 
associated  with  each  value  of  £  approach  (in  some  sense) 
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the  known  exact  solution  when  t  0>  It  is  then  assumed 
that  the  various  functions  entering  into  the  problem  may 
be  expanded  into  power  series  in  &  .  The  series  are  sub¬ 
stituted  into  the  equations  and  boundary  conditions  and 
grouped  according  to  powers  of  t  •  The  coefficients  of 
each  power  then  yield  a  sequence  of  equations  and  bound¬ 
ary  conditions,  the  coefficients  of  i  giving  th«i  first- 
order  theory,  those  of  g1  the  second-order  theory,  etc. 

As  an  exact  initial  solution  it  is  usually  most  convenient 
to  take  either  a  state  of  rest  or  of  uniform  motion.  Var¬ 
ious  choices  of  £  will  be  made  in  the  applications  later. 

The  shallow-water  approximation  differs  in  that  a 
change  of  variable  involving  the  expansion  parameter  is 
made  initially.  This  introduces  t  into  the  exact  equa¬ 
tions.  When  the  power  series  expansions  are  introduced 
into  the  equations,  the  resulting  equations  of  the  sequence 
are  linear  in  quantities  of  the  same  order,  but  the  equa¬ 
tions  are  too  degenerate  to  determine  all  these  quantities 
without  recourse  to  the  equations  of  next  Higher  order. 

This  leads  to  nonlinear  equations  for  the  desired  functions, 
but  ones  of  a  type  which  have  been  intensively  investigated. 
In  this  case  the  procedure  is  perhaps  artificial  in  that  the 
perturbation  scheme  is  devised  to  lead  to  a  special  set  of 
equations  for  a  first-order  theory,  derived  originally  by 
quite  different  reasoning.  However,  in  doing  this  it  makes 
clear  the  nature  of  the  approximation  and  gives  a  systematic 
procedure  for  finding  higher-order  approximations.  It  is 
instructive,  in  this  connection,  to  read  the  usual  deriva¬ 
tion  as  given,  for  example,  in  Lamb  [1932,  pp. 254-256]  or 


ill  i  III  rd^liliiilLH.IMlillllll1  lie'll 'III 


i  run. 


33 

Stoker  [1957,  pp.  22-25]  (who  also  gives  the  one  given  here). 
10*.  The  Infinitesimal -wave  approximation 

We  i» hall  derive  the  equations  of  motion  and  the  free- 
surface  or  interface  boundary  conditions  for  this  linearized 
theory  without  identifying  explicitly  the  parameter  £  used 
in  the  expansions.  Later  on,  when  specific  choices  are  made, 
the  boundary  conditions  on  certain  geometric  boundaries  as¬ 
sociated  with  the  choice  of  t  will  be  modified  to  conform 
with  the  linearization. 

Consider  two  inccxnpressible  viscous  fluids  in  contact 
along  an  interface  represented  by  >  -  "1  (  *  Quantities 

referring  to  the  upper  fluid  have  subscript  2,  those  to 
the  low^r  fluid  subscript  1;  the  coefficient  of  surface 
tension  is  7  .  Assume  the  following  expansions  in  the  par¬ 
ameter  L  : 


t  V- 
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K  ( )<  ~  J  jJ/'  4  f:  * 

'/  *  7  r'J  +  cv7  1  ♦  ••  '  • 


(10.1) 


Substitute  these  expansions  in  equations  (2.2),  (2.5),  (3.1), 
(3.3),  and  (3.7),  rememberirg  in  addition  that  formal  ex¬ 
pansions  of  the  following  sort,  for  example,  hold: 
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Collecting  first  che  terms  independent  of 
from  (2.5)  and  (3.3) 
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Collecting  the  coefficient*  of  the  flret  power  of  t 
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If  the  upper  fluid  is  replaced  by  a  given  atmospheric 

pressure  distribution  ii  f  *  ,  ■»  t  s  )  then  the  equations  for 

1  a  f  t W 

the  lower  fluid  become  (propping  the  subscripts) 
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For  convenience  we  have  assumed  above  that  the  expansion 


for  c.  starts  with 


I 
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If  we  had  assumed  instead 
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,  we  would  have  found  from  (3.1)  and  (3.7) 


the  equations 
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(<s 


(0,1 
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and:  hence,  *},&'  *  const.  The  zero  values  of  /  in  the 


boundary  conditions  would  then  be  replaced  by  this  constant 


Taking  the  constant  equal  to  zero  means  that  we  have  taken 


the  undisturbed  interface  as  fx,  Vj-  plane. 


The  equations  above  give  the  linearized  equations  of 


motion  and  boundary  conditions  at  the  interface  or  free 


surface.  If  one  now  proceeds,  as  we  shall  not  do  for  this 


case,  to  collect  coefficients  of  £u  ,  one  may  obtain  the 


differential  equations  and  boundary  conditions  for  the 


second-order  corrections  to  be  added  tp  the  solutions  of 


the  linearized  equations,  and  so  forth  for  higher-order 


corrections.  In  general  the  resulting  equations  are  too 


unwieldy  to  be  useful 


A  special  case  of  the  linearized  equations  which  is 


of  particular  interest  is  irrotaticnal  flow  of  a  perfect 


fluid.  There  is  then  a  velocity  potential  which  we  as¬ 


sume  has  the  followirg  expansion 
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(10.5 


Condition  (2.11)  becomes 


&  0<lJ  “  O  /  **1,1, 


(10.6 


Let  there  be  two  superposed  fluids  with  velocity  potentials 
£$>,  «nd  describing  the  motion  in  each;  otherwise  the 
same  notation  as  above.  Then  condition  (3.1)  at  the  inter¬ 


face  gives  the  linearized  condition 
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and  condition  (3.9),  together  with  (2.10*) ,  gives 
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The  further  special  case  when  both  the  upper  fluid 
and  surface  tension  are  missing  will  be  dealt  with  so 
often  later  on  that  we  repeat  the  boundary  conditions  for 
it.  We  allow,  however,  a  pressure  distribution  on  the  free 
surface,  j7  f)  ss  t  p'"  k  •  ■  •  .  The  first- 

order  boundary  conditions  are 
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Eliminating  between  the  last  two  eojaations,  one  gets 
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The  boundary  conditions  for  the  second-order  corrections  are 
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not  too  long  to  write  down: 
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Eliminating  *j  and  '*  from  (10.11),  one  finds  a  counter¬ 
part  to  (10.10): 

i  w  v ■'  r/'V:  htt )  (10.12: 

Under  certain  circumstances  the  next-to-last  term  will  vanish. 

The  boundary  conditions  for  higher-order  corrections  will 
not  be  worked  out  in  detail.  However,  they  are  o£  the  form 


f'f','  =  A;\$  ,■■■.$ 


i  j  '  *  *1 


j, /i ^t"V, ♦(‘V" 


(10.13; 


where  A;  a  ,d  J34  are  functionals  of  the  functions  in 
brackets,  in  this  case  complicated  polynomials  of  the  func¬ 
tions  and  their  derivatives  evaluated  at  ^,  —  0. 

It  is  useful  to  have  the  form  of  the  linearized  boundary 
conditions  when  certain  additional  assumptions  are  made. 

First,  let  us  suppose  that  the  (x,  5, ^-coordinate  system 
is  moving  with  velocity  cl*)  in  the  A- direction  with  respect 


to  the  fixed  (/  ,  ■  ) -coordinate  system.  Then,  from  the 
equation  following  (2.15)  with  j  :  <$• 


A 

+  - 
M 


and  the  boundary  conditions  become 


—  m 


C  (>'.». "  c'i-  -  cg: 

•  (it  -I  -  (0  _  _ 

-  4 15  pf  -Cf  P  .  -o  . 


(10.14) 


If  t  is  constant  and  the  motion  is  steady  in  the  moving  co¬ 
ordinate  system, 

<f  ‘  (j?.  I  p 


and  the  linearized  boundary  conditions  are 

>•  -  «*'  =  *■»  r\  —  o 


}  +  1,  0  P  *t  r  f*  !*.  J)s0- 


•I  -  0 


(io.i5: 


If  the  motion  is  steady  with  respect,  :.o  a  mov:.ng  co¬ 
ordinate  system,  one  may  impose  a  uniform  flow  in  the  op¬ 
posite  direction  and  then  treat  the  problem  as  a  steady  one 
in  an  absolute  coordinate  system,  but  carrying  out  the  per¬ 
turbation  about  the  uniform  flow.  We  illustrate  this  for 
die  case  of  two-dimensional  irrotational  flow.  Let  cp(*,  ^j) 
and  be  the  velocity  potential  and  stream  function, 


■ 


«*«  III  a.«^44«4l 


respectively,  end  assume  expansion#  of  the  form 

i  *  -  c  y  -  2  -f  m j)  *  ‘ 1  cf 1  * 

't ( *>  ^  :  * L  'j  *  k  k 1  ( <  -j)  ■<  «'  •/ * 
*  j  ( *  >  —  1 1  ,  *■  i  'j1  ' 


(10.16) 


The  differential  equations  -  c 
become 


<  <  M'»  ■  rf>  cf 


ILf 


.  i 


4  f 


<  n 


(10.17) 


The  kinematic  condition  (3.1)  Is  replaced  by 

'df  (  ■<  ^  r  ^ 

Substituting  the  expansions  (10.16)  In  this  equation  and  In 
(3.111),  one  finds  from  the  coefficients  of  : 


-  j  ""  (*r)  ^  C  -  i  ;  -> 

Gj  j  '>  ;*\  -c  0  +  f  '  p  "(*)  =  o. 


(10.18) 


Eliminating  ■V'  and  using  the  third  of  aquations  (10,17), 
one  gets 

^  4  c  f  f 


Collecting  the  coefficients  of  <*v,  one  obtains  after 
some  manipulation 

v  j1  -gO  *+n> 

(10.20) 


C  >J  *>  -  o)  +  ^ 

here  we  have  assumed  for  simplicity  that  jo  ~  Q  , 
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10  3  .  The  shaliow-war.er  approximation 

Tills  approximation  has  been  widely  used  by  hydraulic 
engineers  In  the  study  of  open-channel  flow  and,  in  a  further 
simplification,  is  used  for  the  theory  of  tides.  In  deriving 
the  equations  from  the  exact  ones  we  shall  follow  the  method 
of  Friedrichs  [1948]  and  Keller  [1948],  However,  a  somewhat 
different  approach  to  this  approximation  due  to  Ursell  [1953] 
is  also  instructive.  Although  it  is  possible  to  carry  through 
the  derivation  while  taking  account  of  surface  tension,  this 
will  not  be  done  here.  It  will  be  assumed  to  start  with  that 
there  are  two  perfect,  incompressible  fluids  with  an  interface 
f)  ;  the  bottom  fluid  is  bounded  below  by  a  rigid 
surface  ^  =  $(*  j).  Variables  pertaining  to  the  lower  fluid 
have  subscript  1,  those  pertaining  to  the  upper  fluid  sub¬ 
script  2.  The  motion  will  be  assumed  irrotational . 

Before  making  an  expansion  in  powers  of  a  parameter,  it 
is  essential  to  make  a  change  of  variable  in  which  vertical 
and  horizontal  distances  are  stretched  by  different  amounts. 
Let  wi  be  a  scale  for  horizontal  measuremen.  nd  n  one  for 
vertical  measurement.  Define  t  ~v*/v»\x  •  Introduce  new 
variables,  i  ,  by  the  equations 

***&,  pp.  (io.21) 

Equations  (2.2),  (2.6),  (2.8),  (3.1),  (3.9),  and  (4.1) 
(with  T  =-0)  become: 


where  w,C-,w,  j", ,  (  possess  suppressed  subscripts  1  and  2  for 
the  lower  and  upper  fluids  respectively,  except  in  the  last 
equation . 

Now  assume  expansions  of  the  form 


r,  = 


X  1 


+  I  ¥. 


t 


(10.23) 


f  '  /  *  i 

substitute  in  the  equations  of  (10.22)  and  collect:  according 
to  powers  of  t .  (We  shall  henceforth  suppress  the  bars  on 
y,  5  f  f  .K  )  The  terms  independent  of  £  give  the  equations 

J  j  *  '  ' 


ir^Va"-  o  t 

1  ^  U 

fet  (•) 

tc' 

z 

II 

(10.24) 
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The  first  and  fourth  equations  give 


V  i  J  £  c 


(10.25) 


The  third  then  states  that 


I  ,  >,  /£») 

,J  ’  iW  ro  tr  u 

i  a 


‘'-Vs  v -  rj 


(10.26) 


The  terms  which  are  coefficients  of  *  give,  after 
making  use  of  (10.25)  and  (10.26), 


'«'!  (j  ,'r  ('  >  — 

u,  -t  UV  i  *  y  ^  U 


j 


/Uj  0  [0  "Ol  ro  (.?> 


t  +  W  u  ,  t UJ  •  J 


:vr 


o 


L- 


r  o , 


vK'.'r  =- . 

ff“  +  o  uv/  +•  ur'C1  u/,  -  n  .  /  f 

io  y*  ^ ^  -•/"  —  o  •f»r  > r  7 

+  l  -V  -j  t- 


(10.27) 


a  ' 


(The  equation;’'  deriving  from  lrrotational lty  and  the  contin¬ 
uity  of  pressure  will  be  brought  .in  later.)  The  first  and 
last  two  equations  of  (10.27)  together  with  (10.26)  give 


'■M 


r«‘»  ~  .|N,,  *  ‘-'.i  )  *  (U 


i  * 


-i  -6  J  i 


(10.28) 


r"  -  J *K  •/  7  /•  ■  •  /  i-  r  7  Jj 


The  third  equation  of  (10.27)  gives 

f'"=  -  fp  1  ’  (,'a  r!- 


■+  (*,*VJ*  *  /u-'  ^ 
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In  order  to  evaluate  4  ,  further  Information  Is  necessary. 

Here  are  two  cases  of  interest.  1)  If  the  upper  fluM  Is 
absent,  the  condition  ^  /,  ,  r  / )  ~0  gives 

T  =  '!!}>  4  *'&t]  (10.29) 


2)  If  the  upper  fluid  is  unbounded  above,  then,  up  to  an 
additive  constant, 


p.w=  -f.u  ♦  ‘(•Myr  +  i, 


(10.30) 


If  the  upper  fluid  is  bounded  above  by  a  free  surface 


±  x,  j.f)  =  i  a  -  •• 

the  boundary  condition  ^  /"<  ^  r 

f,rV  Vi  V*  M  ,  /  ^  with 


,  then  one  may  satisfy 

.  ;  ; .  r)  ~ 


!<•! 


tr  =  -f. 


-p  ^Jrr) 


i1 


I  J. 


(10.31) 


It  Is  clear  from  the  form  of  pr"'  why  the  shallow-water  ap¬ 
proximation  is  sometimes  called  the  hydrostatic  approxima¬ 
tion. 


The  usual  equations  for  the  first  approximation  to 
the  shallow-water  theory  are  those  In  which  only  the  lower 
fluid  is  present.  They  may  now  be  obtained  by  substituting 
(10.29)  in  the  second  and  fourth  equations  in  (10.27)  and 
(10.28)  in  the  fifth  equation.  They  are  (10.25),  (10.29), 
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,  M)  i«  _  , 


4,  ^  ^  ^  +  1  i* 


^  +--M1  +  3  7v  s  0 

?r*  l  f7,‘,-^j,  •  r  ^  rv,7*  -*jj.  -~°' 


(10.32) 


If  one  now  collects  the  coefficients  of  and  the 


remaining  coefficients  of  C  ,  one  finds  after  some  reduction 


in  ,  f ')  .  i  ^ 

4  ,  -  -*■  ,  -  ^ 


('I  ft)  I 


M,  t  "*  4  4  j  -f  u,  i  4  i/" 


*  *  f>-  /f 


hi  H  H;  jli  f*l  J'l  .  Ju  /,  — 

vr  4  u  »  -tu.  ,  v  v-  v*  P*s  /f  “  0j 


|ii  fii  f .  \  foi  '•)  (■ 1  fw  n,  c.  hi  ID  I  _  rt 


(10.33) 


lit  .  Jh 


“  J  '  WI  '  *“*  /  r»  '  ? 


rf;  t-  rf***1 


f.  o 


Some  relations  can  be  derived  immediately  from  these 


equations.  For  the  sake  of  brevity  we  introduce  the  follow¬ 


ing  functions: 

(rt  ft)  ,  (  ,(•>  Ml 

-Mw.  ■;  ), +  {'■'{  i  J, 


,  B,=  C»  -c,  +  Si. 


Equations  (10.28)  may  then  be  written 


ir/"  =  -  A  ;  +  3;  j  <  *  i- 1 


(10.28') 
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Than  the  fifth,  first,  and  third  equations  of  (10.33)  give 

r  x,  £yr). 


m  _  i  _ 


i  i  A  ■  +  ,f), 


u,  =  u 


(10.34) 


f/f  =  t  -'X  -  -  -t;  ♦  /[/ib  4^  -u,-'V 


where  we  have  suppressed  the  subscripts  indicating  the  fluid, 
The  rest  of  the  equations  and  the  boundary  conditions  are 


still  available  to  determine  the  unknown  functions,  'ife 


carry  this  out  only  for  the  case  the  upper  fluid  is  missing. 
Then  the  last  condition  in  (10.32)  becomes  \J"f  V  •' 


T  ‘  ich  allows  one  to  determine 


after 


is  found. 


The  next -to- the -last  condition  in  (10.32)  determines  .<  1  *  ), 


The  equations  for  t ,  ,  and  ^  are 


-f  UT  ' ’ * t-l ,  *lr  [  -  -*■<  'K  =  "*  f  M  "  -3 3 ,  ^ 

i  '  * 

-v  L  +  If.  *  wx  v.  :  -  t  .  -  3  i  j  > 

J  *  }  *  f  °  (10.35) 

,  d'/ 

k,  =  , 


where 
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The  solutions  to  these  equations  give  the  second  *orde\.  cor¬ 
rections  to  the  first-order  shallow-water  thec-y. 

The  equations  resulting  from  the  coefficients  of  l 
have  been  given  by  Keller  [1948]  for  two  dimensions,  but  will 
not  be  reproduced  here. 

The  equations  (10.32)  for  the  first-order  theory  are 
nonlinear.  In  the  theory  of  tides  and  seiches  it  is  custom¬ 
ary  to  simplify  further  by  linearizing  them  in  a  manner  sim¬ 
ilar  to  that  used  in  deriving  the  equations  for  the  infin¬ 
itesimal  -wave  theory.  Let  ;  =  1,  be  the  surface  of  the 
undisturbed  water  and  assume  that  one  may  make  further  ex¬ 
pansions  In  a  small  parameter  *  1  ' '  +•  u;  ‘'k  * u/* 

1  z.  v  ■<•  '  .After  some  easy  manipulations 

I 

one  finds  for  the  linearized  approximation  to  (10.32)  the 


equations 


IttO 


•;  i* 


r  \ 

O  -4 

J  t 


(10.36) 


If  the  bottom  is  flat,  the  equation  for  ¥J  becomes  the 
simple  wave  equation. 


"NOTICE;  When  Government  or  other  drawings,  specifications  o 
other  data  are  used  for  any  purpose  other  than  in  connection 
a  definitely  related  Government  procurement  operation,  the  U.S 
Government  thereby  incurs  no  responsibility,  nor  any  obligatioi 
whatsoever;  and  the  fact  that  the  Government  may  have  formul 
furnished,  or  in  any  way  supplied  the  said  drawings,  specificatio 
or  other  data  is  not  to  be  regarded  by  implication  or  otherwise 
in  any  manner  licensing  the  holder  or  any  other  person  or  corj 
tion,  or  conveying  any  rights  or  permission  to  manufacture,  use 
sell  any  patented  invention  that  may  in  any  way  be  related  the 


# 
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D.  Theory  of  Infinitesimal  Waves 


This  chapter  will  deal  with  special  solutions  of  the 
linearized  equations  derived  in  section  1,0  x.  This  approximate 
theory  has  been  very  fruitful  in  its  application  to  problems 
with  various  boundary  configurations  the  linear  character  of 
both  the  equations  and  boundary  conditions  allows  one  to  use 
easily  found  simple  solutions  to  construct  othur  solutions 
satisfying  special  boundary  conditions.  The  derivation  of  the 
equations  in  section  'Dot  suggests  the  Limitations  of  their  use 
in  physical  p:  o’u  u  m  < :  If  L  and  V  are  a  typical  length  and  ve¬ 
locity  associated  with  the  physical  problem,  then,  when  the 
perturbation  parameter  t  is  small,,  the  surface  elevation  and 
velocities  (or  their  deviation  from  a  uniform  flow)  should  be 
small  with  respect  to  L  and  V  respectively.  The  smallness  may 
riot  be  uniform,  but  the  quantities  in  question  should  approach 
zero  point-wise  with  &  except  at  singular  points. 

11.  The  fundamental  equations. 

With  few  exceptions,  this  chapter  will  be  concerned  with 
the  Solution  oi  a  problem  in  potential  Lhcoi'y.  Let  the  ( x ,  z)- 
plan©  be  at  the  undisturbed  free  surface  We  shall  be  seeking 
a  function  ,  the  velocity  potential  of  the  motion, 

satisfying  the  conditions  (cf.  10  10) 


1 1  e  ' 

on  solid  boundaries, 


(11. 


where 


4 . ? 


is  to  be  satisfied  at  all  nons ingular  points  or* 
the  fluid  in  the  region  >.  -j  and  is  the  normal  velocity  of 
thtf  solid  boundary  at  .1  given  point  p  ,,  1  ♦  is  a  given  pres¬ 

sure  distribution  on  the  Inn  sin  Inc,  in  many  problems  it  will 
be  0.  The  form  of  t lie  free  .suriate  is  given  by: 


/ 


(11.2) 


Two  special  cases  occur  frequently.  If  the  motion  is  steady 
in  a  coordinate  system  moving  with  constant  velocity  e  in  the 
x-di reel  ton,  then  with  x  y ,  z  as  moving  coordinates,  the  free- 
surface  boundary  condition  and  equation  of  the  surface  are  given 
by  (cf.  10.  15) 


■  J  1  cT J 


t, 

>  > 


1".,'  -  JrU'  '  r - T 


If  (p  and  jr  are  harmonic  functions  of 


'he  time,  i . e. 


(11.3) 
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where 


“J »  >) 


and  similarly  for  j*i  ,  then  the  free -surface  condition  and  ©qua 


tion  of  the  surface  become 


f?ty  1  -  TT  -  ' 
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(11.4) 


In  the  few  cases  where  we  consider  superposed  fluids, 
viscous  fluids  or  surface  tension,  we  shall  refer  back  to  sec- 
t,  i  o  n  1 0  f  o  r  t  he  e  q  u  a  t  i  o  n  s 

Uce  of  complex  variables.  For  two-dimensional  i rrotat iona l 
motion,  it  Is  frequently  advantageous  to  use  complex  variables. 
Le  t 

%  -  1  -  ~  i  '  '  r  '  l' 

where  and  /'  are  velocity  potential  and  stream  function, 
respectively.  (It  should  be  cleat  from  context  whether  z  ir 
being  used  for  >  *  <  >-j,  or  one  of  the  horizontal  coordinates  ) 
Since  the  equations  relating  and  Y  > 

J  -'V  i.  > 

'  'I 

are  just  the  Cauchy-Rtemann  equations,  the  function  *  t  is 
an  analytic  function  of  z  for  all  points  z  for  which  >,  and  i 
exist.  T  ,,  til  will  be  called  the  'complex  potential".  The  "com¬ 
plex  velocity1'  is  given  by 

“  ;t(bVt  =  ^  V' 

The  boundary  condition  at  the  free  surface  in  (11.1)  can 
be  expressed  in  the  following  equation  in  -f 
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A  a. 
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o  r 


(11. 5) 


The  first  equation  of  (11.3)  becomes 

=  j  r, = °-  (u'6) 

However,  equation  (10.19)  shows  that  this  may  also  be  taken  in 
the  form 


?S  j  -+  cv  f(\)  j  =  f  1»W  hr  3=0. 


.Ill 


If  one  may  express  -  ,  ^  ='  f c»s  trf +<fl(^)w  then  the  first  of 

equations  (II  4)  be cones 
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a  t 


f  0 


ft  -  ■■  i) 


ft  (*) 


(ll. 


We  note  that  in  order  to  express  in  a  manner  analo¬ 

gous  to  that  used  for  0  immediately  preceding  (11.4)  one  must 

Introduce  a  second  complex  unit  j  which  does  not  "interact” 

,  ,  ,  ,,  ,  m  IT  t 

with  i.  Thus  let  **'>  •  *,  ,  ■*)>.()).  Then  h  Pt  =  Te(  if},  2  ' 

If  t  is  an  analytic  function  satisfying  any  one  of 

f  (*)t  - 

the  cone  it  ions  (11.5)— (11.7),  then  r  >  will  also  satisfy  it. 


12.  Other  boundary  conditions 

The  boundary  conditions  given  in  section  11  will  not  ordi¬ 
narily  be  sufficient  to  ensure  a  unique  solution  to  the  pro¬ 
blems  in  which  the  fluid  occupies  an  unbounded  region.  An  addi¬ 
tional  condition  at  infinity  must  be  imposed  upon  the  potential 
function  In  certain  cases  the  proper  additional  condition  is 
fairly  clear  from  the  physical  problem.  For  example,  for  a 
body  moving  steadily  In  an  infinite  ocean  undisturbed  except 
for  the  body,  It  seems  reasonable  to  Impose  the  condition  that 
the  fluid  motion  vanish  far  ahead  of  and  far  below  the  body. 

For  the  fluid  motion  produced  by  a  stationary  but  steadily  os¬ 
cillating  body,  it  seems  reasonable  to  impose  vanishing  of  the 
mot  i or.  far  below  the  body,  but  outgoing  waves  at  infinity  on 
all  sides,  if  the  body  does  not  extend  to  infinity  in  some 
horizontal  direction,  the  so-called  "radiation  condition" 

If  the  body  is  not  bounded  in  a  horizontal  direction,  one 
may  easily  see  that  the  radiation  condition  stated  above  cannot 


be  expected  to  be  satisfied  For  example,  suppose'  t  hat  waves 
are  being  generated  by  some  type  of  oscillation  of  a  vertical 
half-plane,  say  }~c  ,  >>o  ,  in  whicn  the  oscillator/ 

notion  of  the  half-^-ane  is  independent  of  *  Then  one  will 
expect  the  generated  waves  to  behave  like  outgoing  plane  waves 
from  the  two  sides  of  the  plane  as  )(-»'-<■  ;  these  will,  not 
satisfy  the  radiation  condition  in  the  direction  0/  On  the 
other  hand,  one  might  expect  that  the  influence  of  the  edge  at 
X  r  r  would  show  up  as  waves  satisfying  the  radiation  con¬ 

dition.  The  formulation  of  proper  boundary  conditions  in  situa¬ 
tions  of  this  sort  has  been  discussed  by  Peters  and  Stoker 
[1954];  see  also  Stoker,  [  1956;  1957,  p.  109 f f.] 

In  diffraction  problems  one  customarily  prescribes  the 
form  of  an  incoming  wave  and  then  seeks  *"he  scattered,  wave. 

The  preceding  remarks  concerning  the  boundary  conditions  for 
waves  generated  by  an  oscillating  body  apply  also  to  the 
scattered  wave. 

In  more  complicated  physical  situations  it  is  not  always 
clear  what  boundary  conditions  should  be  imposed  at  infinity, 
and  errors  have  teen  made.  For  example,  for  a  body  which  is 
both  oscillating  with  a  fixed  frequency  cr  and  moving  with 
a  steady  average  velocity  c  ,  one  might  reasonably  expect  no 
motion  far  ahead  if  C  is  large,  but  a  radiation  condition  if 
C  is  small.  However,  the  formulation  of  the  boundary  condi¬ 
tion  cannot  be  completed  until  the  problem  is  partly  solved. 

The  proper  formulation  of  the  boundary  conditions  at  in¬ 
finity  can  frequently  be  obtained  by  a  method  recommended  by 


Havelock  [1917,  1949a]  and  used  also  by  Brard  [ 1948a,  b], 

Stoker  [1953,  1954],  Stoker  and  Peters  [1957],  D*  Prima  and  Wu 
[1957],  Wu  [1957]  and  others  It  consists  in  formulating  an 
initial-value  problem  for  which  the  desired  steady-state  prob¬ 
lem  is  the  1 i nil  as  t~  >  For  the  initial-value  problem 

the  boundary  condition  at  infinity  is  that  the  fluid  motion 
vanishes  everywhere. 

13.  Some  mathematical  solutions. 

Some  of  the  mathematical  solutions  to  be  derived  in  this 
section  will  provide  solutions,  without  further  modification, 
to  certain  physical  problems,  others,  although  apparently  not 
acceptable  physically,  will  provide  fundamental  solutions  which 
can  be  used  in  constructing  solutions  to  other  more  complicated 
physical  problems.  In  all  cases  the  fluid  is  assumed  unbounded 
in  a  horizontal  direction  and  either  infinitely  deep  or  with  a 
horizontal  bottom  ^  [  :  the  pressure  on  the  free  sur¬ 

face  is  taken  to  be  zero  everywhere.  The  solutions  without 
singularities  are  obtained  by  the  method  of  separation  of  vari¬ 
ables,  and  are  all  harmonic  in  f  .  It  will  not  be  necessary 
to  carry  along  the  subscripts  of  (11.4). 

L/  Separation  of  the  y-variable 

Assume  that  one  may  express  Y  by 
,J[  {  *>  i)  v  “  ^  ^  j)- 

Then  & ,  cP-  4,,  \  becomes,  after  separation, 

1  '  ’  1  Ms 

^  i  Cf  t  A  cf  o  ,  Y  '  -  A  Y  =  . 

The  two  cases  A  -  >  0  and  A  -  -  tv  ‘  <  o  lead  to  dif¬ 

ferent  solutions. 

"k 

;  •  .j|M  ,j»!.  ill)1  •'  i  ;  •!  • 


A  >  -  In  .this  case  -k  ' ■*  7  satisfies  the  wave 

equation 

Z  'S  t  <*  w  —  C 

and  Y  is  given  by 

y  =  a  e"  ■*  t  4  e.  '  '• 


If  the  fluid  is  infinitely  deep  and  -f  y  1  , 
remain  bounded  as  >-»-./■  ,  one  must  have 
Then  condition  til. 4)  requires 


is  to 


/ 


and 


'  ? 

( x  ^  is  of  the  form 


cf  !  *.  - 


( 13.  1) 


If  the  fluid  is  of  finite  depth  h  ,  the  boundary  condition 
■  =  =  requires  Y  to  take  the  form 

Y  -  a  c  c  ^ 1  - +  f ) 

and  condition  (11.4)  becomes 

yv\  tu«Y  <'  h  -  '  y 

an  equation  with  two  real  solutions,  say  t.  -vy,  In  this 

case,  one  has 


'j  ~  i'  A  *  v  )  t? 


a 


(13. 2) 


We  note  that,  if  in ,  <  ty  then  ■' /j  <  ivy4,  i  "Y  .  Also 

*  •"»  -  and  -*  »s  ii-*oo  . 

A  <  0  In  this  case  satisfies 

tf  ~  cf  =■  O 

,  „  .  I 

. . ...L.  ,i  ,  i"  i  i  i  \  ,li„i  ,  ' . . .  : !.:  1  ■  !  ■:  III  lY  i 1  I  in  r  iifol  1 

. . Ill*  Itawii.  . MMIIMM . .  “I  . . . -i;H«*  . . ,|J<.l.|.,|T  ■  '  '  "*1 . -  *i|imfW%4l! 


and  is  given  by 


ft 

Y  _  ta  tv'i  '•<  j  4  t-  i"' «"  ^  • 

Condition  (11  4)  restricts  v  further  to 

Y  ~  _  ,v  t  j  i . .  ,  -»  V  L '  '  '  <*  y 

i 

If  the  fluid  is  infinitely  deep,  requiring  — -  to  remain 

bounded  imposes  no  further  restriction.  If  the  fluid  is  of 
depth  li  ,  then  ,  1  •  t  -C  requires  m  to  satisfy  the 
equation 

l  V  i  /  e)  >  ■  »*  !  —  —  j  J Y 


an  equation  with  an  infinite  number  of  real  solutions,  -  . 

r  ,  In  this  Utter  case  one  may  conveniently  take 

f  in  the  form 


The  roots  ro  j  satisfy  r  II  f?- ijir/i,  r*  j.  <  -fan  v  ,  For  fixed 

U}  rv^  ^  ->  it  t  as  h  -»  *,  for  fixed  ^  /{  ~ 

as  k  — *  0  ,  *  r  d  ‘v. .  h  — »  ^ '  i  t,- 1)  n  as  r  ~a>  °G 

For  these  two  cases  one  finds  then  for  <-(-{*  'i  i)  the 
forms : 

infinite  depth : 


(13.3) 


finite  depth: 

cf  (x.'j,  J.)  =  L  c  «s  'ufH+li, 

U/)  Further  separation  of  variables. 

We  now  assume  and  substitute  in  each 

of  the  two  equations  for  Lf  given  above. 


A  >  o 


In  this  case  subst  itution  in  -t  •* ."-f  -  J 


gives 


X 


i  frX-ftM/  -  o  Z  +  t‘2  = 


(The  equations  obtained  by  replacing  'X  by  I  V  will  gi  •» 
the  solution  obtained  below  for  '  <  J  ,  with  *  and  inter¬ 
changed.)  The  solution  for  Z  is 


Z  —  t  tus  ■+  X  —  J3  ■  0  • 

The  solution  for  X  depends  upon  the  sign  of  r/*  - 

fal<  wX:  X  =  C  ■  /•  ‘  h’'>  , 


n  V  In-  t 

ft  -  v<r\ 

\Z  >  w 


X  '  ^  -y  ^  J  - 

X  /  ^  t-mv  -  X  ^  k,-t>'V 

A  =  «-  e  -  a  £ 


A  <.  £  Substitution  in  cf  “  >'L  ^  -O  gives 


,  Z  "  4  ^  Z  -  C  , 

which  gives  Z  as  above  and 


X 


v  /  f,v< 


(Substituting  -  L  for  would  give  the  solutions  corres¬ 
ponding  to  A  X'  with  >  and  ^  interchanged.)  We  may  accu¬ 
mulate  the  preceding  results  to  obtain  the  following  fundamen¬ 
tal  solutions : 
for  infinite  depth: 

Z  J  «  co%  +  £  >^v -4‘  jcos  H'Xosfcrt-*  f), 


g.v2  (  a  x  +  £)  cos  f'Hj  -mO  t  csto’tf-  r)  ,  ^VrvV) 


(13 
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mylMWiNMiMHfi 


e  e 


f  x;  /<\ 


(m  COS  **  ->.  +  V  S'r  rr  ^  C  .*  / 

where  -v  ■=  <rv/  ^  ; 

for  finite  depth: 

, _ ,  / — 

tail,  s  1  t  '  a  C,C>S  *  s'  '.X 


-  J>  ,/  {}' 


cos(k’+(  tosfit*  r) 


t ■  j f *■  t  f  / i-o 5  £r>  ■» j;- 


<:  j/rN*  , 


4  / f}X~-  fv\  c 


yj *<  -**#■ 


C  o  S  k  l’>« 


'  ,  //. ,  ,  j>  •  /**■ 51,1  ( 13  • 6) 

zS  'I?  •+■  jc  v  ' 

/  J 


||’  «>  v" 


where  fis.ta**  wt.ii  -  tV 


and  ** 


ffl  n  w  ^  4  5"  /5  -  ^ 


The  corresponding  solutions  for  two  dimensions  may  be 


obtained  toy  setting  f- 


and  deleting  the  second  and  third 


equations  in  each  group. 


For  either  set  of  solution.;  only  the  first  in  each  is 


bounded  for  all  values  of  the  variables  for  which  y  4  0  or 


~  [  <  -i  <  O  For  two-dimensional  motion  it  has  been  shown 

4.9^ 

by  Weinstein  [  1927|  t'  t  the  only  function  harmonic  in 


and  satisfying  (11  4)  and  =  C  for  which  both  Lf  and 

CP  are  bounded  in  -  l,  £  '/  <  O'  is  <P  -  f\  con  f- 1 u)Si *(•  ' *  *•*)  • 

Tr 

Keldysh  [ 1935 ]  and  Stoker  [1947,  pp.  7-9]  have  proved  a  similar 


•*  "i*  2m 

theorem  for  the  lower  half-plane:  If  Cf  and  %  *  '■f  ?  are 


bounded  for  y  <  as  aVj'*  ->  ^  ,  the  only  c p  satisfying 

(11.4)  and  harmonic  everywhere  in  the  half-plane  y  40  is 
f\  &  i  *  (  -*0  ■  Weinstein's  theorem  has  been  generalize 

by  John  [1950,  p  59]  to  three  dimensions:  If 

5 

satisfies  (11  4),  <-f  (  A ,  h  j)  , 


r  7 


tos  «x,  \  R  -  1  r  c 

and  is  haimonic  everywhere  in  -  r  *  5  C  ,  then  -*■ '  K,  ^ 

is  of  the  form  (L3.2)  with  <  an  everywhere  regular 

solution  of 

^  ;.  Y  t 

The  condition  at  infinity  is  necessary,  as  the  solution  derived 
below  in  1*3.8),  ^  .< ,  -  shows.  The  corres- 

ponding  theorem  for  infinite  depth  was  proved  by  Kochin  [1940], 
The  equations  for  cf  {  y, .  ^  ^  may  also  be  separated  in 
polar  coordinates  (R,  b)  f  /-  P  i%  «  r  *  K  S»  •*  at  .  We  give  only  the 
so lu t ions ; 
infinite  depth: 


Y.  '  "  K  ,!  •  f  '  *  £  -  - 


(13.7) 


where  v  -  /  and  o  is  an  Integer, 

finite  depth. 


tisi  *•;>.?)• 


'J 


-  ‘T  /'•<«►/) ‘••skt-*’*),  lii,  (13.8) 


where  r*-,  .-77,-0,  to*  ^  4~/s  =5  and  is  an  inte¬ 

ger.  Here  J  J  Y,  ,  J  „  ;  ]<Y  are  Bessel  functions  (we  use 
Watson’s  notation).  Yv  and  rY,  are  boch  singular  at  I?  r  c 


but  approach  zero  as  — 1  *  ;  u- .  and  i  „  are  both  finite  at 

R  ~  ^  ;  approaches  zero  as  <  ->-c  ,  increases  ex¬ 

ponential  ly . 

Singular  solutions 

In  this  section  we  shall  find  solutions  of  the  problems 
set  in  section  11  which  have  singularities  of  simple  type  at 
a  single  point.  We  shall  indicate  proofs  only  for  the  case 
of  simple  sources,  i  e.  singularities  of  the  form 


.  u 

i. 


We  shall  consider  first 


'<■  _  or  log  „  *'->  - 
the  case  of  a  stationary  source  of  pulsating  strength,  then 
the  case  of  a  moving  source.  Three-dimensional  problems  are 
trea  ted  first. 

Source  of  pulsating  strength  m  three  dimensions. 

Let  (a,  b,  c)  be  in  the  lower  half-space.  We  wisi:  to 
find  a  function 

n  -  'y,  ■  ‘I.  ("  -  '  r~ 


f' 


defined  for  / 

i)  l  or, 


except  at  (a,  b,  c >  and  satisfying 
except  at  (a,  b,  c),  ‘  ~  ‘  L  , 


2)  '-f;  ,  I't)  \Y  V 


/ 

3)  i(' 

whei'e  ]  is  harmonic  in  the  whole  region  /  c  O  , 


(13.9) 


4)  0rhi  cP  ■  -  C  ,  <  '  ’/  - 

/  o 


5) 


r. 


R- 


P$(  +»£f,  -o,  /?  I  -V  <¥•*',  =o. 

s  r<  S'  i 


'Till . r . 


He  re  r 


■(/.  +  v-t  ^-e'w  and  5-  X-a,  -v  1,-c) 


Condition  5) (  usually  called  the  "radiation  condition",  re¬ 


quires  the  waves  at  infinity  to  be  progressing  outwards  and 


iuposes  a  uniqueness  which  would  not  otherwise  be  present. 


However,  other  such  conditions  could  be  imposed. 


We  assume  that  a  solution  can  be  found  in  the  form 


&(*r  ■  ,1  i:  _  -  -f.(v  ...  :-<-rT  *  ^ 


v  -  t  r , 


(13. 10) 


will  be  determined  at  the  end  so  as  to  satisfy  3) 


note  the  double  Fourier  transform  in  >  and  ^  of  by 


Afi 


f  Jt  ,  J  .r  -  ■  ;  <V  v» 


/  ;■  j. 


Then  condition  l)  app,.?^  to  cfc  becomes  after  transforming 


■Tu  fJ 


u  •- 
ft  h£. 


S.  z.  A  Ah 


(13. 11) 


where  we  have  used  4)  to  discard  the  other  solution, 


From  the  known  integral 


oe  r'  _^|-|  /  b  (Xc*o  S|,A  *) 

'  *  __  r  \  £  e  /-JJ  i  ^ 

*  —  j  ~7i  j  * 


i  -  i^A 


(13. 12) 


'SiiJ.il  .'Hi  1 1 


■ . 

ijiiiiiiiiuiifiMjinil&iiii 


( 0 


one  may  compute 


-  I V  |  .  ,  z.  .?•>>£  -f I  s*  *vi  J 

e  '  c. 


(  13  13) 


Substituting  _  *  -  in  the  transform  of  condition  2) 

gives 

A  .  ■,  \  .  i  h  [ ^  £  C  i  f  4  c  &  /•  ‘  , 

-{  v  -  V 

(13.14) 


We  now  have,  formally, 


;V  S  - 


iV: 


it 


Sir’ce  the  .ntagrand  has  a  singularity  at  r  V 
it  Jt  meaningful  without  further  definition. 


pret  the  integral  as  a  Cauchy  principal  value. 


% 


f 


i  i'i  I 

i  ■'<  4  .-'-r 


_  X  4  -X  +  —  r—  ?  '  t 

T  r,  *  4  «*r  *‘v 


,  the  integral 
We  shall  i nter - 
Then 

i  J  +1 

C  v*  f  * 

c-Ji  (13  • 


15) 


where  r,v  ;  f  *  V-*  -  ‘jV  i  ' ;  *e)  .  The  second  equa- 

t ion  may  be  derived  easily  from  the  first  one  by  use  of  (13  12) 
suitably  modified.  c^:  satisfies  1),  2)  and  4);  tfa  is  harmonic 


In  the  wholes  region 


In  order  to  satisfy  5)  we  shall  first  find  the  asymptot  Lc 


form  of  <i  for  large  .  With  polar  coordinate 


<  f 


one  may  write  (13.15)  as 


>,j): 


_  *  -  +  ~  ■ 
t  r  * 


i  *'  ,  *  ""  /  f'  5  t  1  f'J  ■  :* 


b  J-rr  * 


!«'  Ji  e  l<i 


w  i  + 1?  ;v 


j'  c  b-y 


CtS^fc'k  C6S  ti)  =s>  ck 


f  j  s  V  =  <-< 


&  - -  fc  V;  ,  a 

„  p-y 


In  the  next  to  the  last  equation  we  shall  change  the  order  of 


integration,  write 


c  t  S  1?  A  k  -  f  5  i  V  n  v  «CS*<>  f<->;  -  ^  -•  h  A  1/ 


and  use  the  following  theorem  from  the  theory  of  Fourier  into 


grals  [see,  e.g.,  Bochner,  Vor lesungen  uber  Fourlersche  Inte- 
grale,  Leipzig,  1932,  ch.  I  and  j  8],  If  »  w)  is  a  differ¬ 
entiable  function  in  (:>,-■-]  ,  if  4"(  Jf3)  ,  ,  exists 

and  if  4r  *  <  /  *  and  4  (*'*  !'<  are  both  absolutely  m  tegrable 


in  [a  ,  3  t  then,  as  1? 


(  -Pf*»  — — ^  * '  ~  <L<  =  -  -f  ili- 

i  x  -  x  *  y  J. 


CyS  ff'i  *-*'■)  , 


X  -  X  a 


(tx  -*  5  /o). 


(13.16) 


113’  i  i i ,  Tr  „  u  ::n wr  'tt1* — t"- 


Homo  in  be  r  i  n  g  t  h  a  t  b  o  t  b 
fin  Is 


1  .  ,  - 1  ,  ' 

and  are  -  ,  one 


—  -  L. 


I  f  v  <  /  '  '  - 


c  ) 


The  asymptotic  exp. ins  i  <>,„  ot  this  j  n  logr.il  is  well  known  [set? 
e  g.  Erdelyi,  Asymptot it  expansions,  Dover,  New  York,  195.1, 
p  48]  and  we  may  write 


-  '4*<  /  l 


Jf-  +  o'-' 


R  / 


If  we  can  find  a  harmonic  function  satisfying  I),  2)  and 

4)  and  having  the  asyn.pt  otic  behavior 


>'y 


-  i  n  v 


-  Vv 


t  he  n 


f  r 


-•  -fi-4  ft  j 

IVvt,  * 


'  1  ”  ^ 


-  «fj 


will  be  a  solution, 
qu  i rements : 


\ 


The  following  function  fulfils  the  re- 


-j  T'  ( 1  a  1  *  \ 

i  Jo  ,  *) 


[see  the  first  equation  (13,7)  and  Watson,  Bessel  functions 
Cambridge,  1944,  p  199].  We  note  in  passing  that  J  has 
the  same  asymptotic  behavior  as 


■In  v  £ 


Vft|.+  (■) 


Il'lMI-Mi  l|l<l|l  #'H  1 1  ■'  |-|i||i||i  IIMI  lllilj  . . . 


* . . . . r . . . . . .  r? . r"w  . . . . . . . . 
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The  final  result  is 


-.1 


■i  * 


— -  +  t-v 

l  ‘C. 


7T-~ 


£f 


i  '3**  i  f  ,  -c-1 

-  v£  \jr  '  l-.SirtTl  ,  ^  “  C 


(  13.  17) 


1  ‘ 


Haskind  [1.954],  using  a  derivation  having  some  similarity  to 
that  used  below  for  the  two-dimensional  case,  has  found  the 

following  form  for  j 

r .  ,  V/"£  ,  *  V1'1 

Cpfx  -  i'T  “t  ~  iv  &  — ej,  -  2.FVC 


-  2  n  v  fe 


T 


J  (v? 

'n->  ' 


( 13.  17*  J 


It  is  sometimes  convenient  to  use  the  complex  form  for  the 
potential,  ,  with 


'-f(  V 


i  -  1 

■1  >  r 


9/ 


ro 


f<  -V 


>0 


■J>P) 


( 13 


17m) 


for  then  |\'£  gives  (13.17)  and  / tf  € 

the  source  potential  for  an  outgoing  wave  with  singularity  of 
the  form  r'S'n  'f  .  Equation  (13.17’)  may  be  written 
analogously.  By  deforming  the  path  of  integration  in  a  famili¬ 
ar  way  one  may  also  express  <J^  (  K>  in  tie  form  (cf. 

Havelock  [1942,  1955]): 


. . * . 


G 


(( ‘  {  "  r  “  /  jb  •  V  C  • :  1 ' :  * ~  J~h  '  :*  J 


,  K 
VS 


•  ’  e  |. 


,-1-*  (  »  , 


,  >  *'  r  ’  -  r 

f,  r  *  -  isve  "  v  vk . 


<13.  17*) 


In  the  analogous  problem  for  finite  depth  it  one  re¬ 


places  4)  by  4*  )  -i  \  -t  I  i  *2  and  proceeds  somewhat  s L ir. i • 

*1  * 


larly.  However,  in  order  to  satisfy  4* )  it  is  convenient  to 
look  for  a  solution  in  the  form 


4 


*  \  ~  • 


■i  . - > 


4  cfc(*'^)JcoS  at  f  J»l 


!  x  ^3  5'  i"i  tr , 


where  v  '  ~  f  .<  -  1  k  4-  '■ 

then  becomes 


J  1  1  k  +  Ir  j  I  I  ^  -  C  j  *” 


Equation  < 13 .11) 
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term  of  (13. 18)  may  be  altered  by  using  the  identities 
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John  [1950,  p.  95]  has  derived  the  fol lowing  series  for 
the  analogue  of  (13.17  ) 
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where  M ^  (  ft  ■>  are  the  positive  real  roots  of 

*"  r^l\  *v  -0.  Either  expression  may  also  be  given  in  complex 
form  as  in  ( 13 . 17M) . 

Potential  functions  satisfying  the  condition  (13.  9),  but 
with  r*  t  rt  in  3}  replaced  by  a  higher-order  singularity 
have  been  giv^n  by  Thorne  [1953]  and  Havelock  [1955],  In 
fact,  Thorne  gives  a  rather  complete  census  of  the  possible 
singular  solu-  ions  for  both  two  and  three  dimension"  and  for 


finite  and  infinite  depth.  Included  are  series  expansions  as 
well  as  integrals.  For  infinite  depth  the  general  expression 


which  includes  (13.17)  is 
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the  associated  Legendre  polynomials  defined  by 
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The  asymptotic  behavior  of  (L3.20)  is  given  by 
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It  has  been  pointed  out  by  both  Havelock  [1955]  and  MacCamy 
[1954]  that  solutions  can  be  constructed  which  vanish  much 
faster  than  this  at  infinity.  Let  the  function  of  (13.20)  be 

t  i  -  ■  : 

denoted  by  j  .  .  Then  -  u  .  ^  is  the 

following  function- 
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where  CoS  0,  -  f  i  >  (?- '/ (  (  \  '  <a-&\ ' f  -  ^  p  +  ,  r)1  .  For 

and  large  fV  these  solutions  are  t-> ”  j  if  -  and  v\  are 

both  odd,  {J  j\  if  one  is  even  and  one  odd,  and  u-(n  J 

if  both  are  even.  Although  they  have  the  form  of  standing 
waves,  they  satisfy  the  radiation  condition  because  they  de¬ 
crease  so  rapidly  with  large  K  . 

In  addition  to  the  papers  cited  above,  one  can  find  treat¬ 
ments  of  the  submerged  source  of  pulsating  strength  in  Koch^n 
[1940],  Havelock  [1942],  John  [1950,  p.  92  f f .  ] ,  where  a 


* 


detailed  discussion  Is  g  i  vc.i  lor  the  case  of  (inite  depth, 

Hasklnd  [ 1944],  and  Li.  [1952]  The  definition  of  the  improper 
integral  In  (13  15)  and  following  is  always  the  same  in 
these  different  treatments  In  some  cases  the  variable  ii  is 
treated  as  complex  and  the  path  of  integration  deflected  arouiu 
the  singularity  /;  -  v  b,  following  a  small  semi-circle  in 

the  lower  half  of  the  >;  -plane  The  radiation  condition  is 
then  automatically  satisfied  if  one  writes  4>  in  the  complex 

-  /  fj  t 

form  ife  J  (Jj>  s  f  y  a  Other  treatments  achieve  the  same 

end  by  introducing  a  "fictitious  viscosity"  <  /-  which  has  the 
effect  of  replacing  the  singularity  at  h  - v  by  one  at 
J  :  *  and  thus  placing  the  path  of  integrat  ion  below  the 

singularity.  In  the  end  one  must  find  the  limit  of  the  solu¬ 
tion  as  ^  j  The  fictitious  viscosity  has  no  relation  to 
real  viscosity  and  may  be  considered  a  mathematical  device  to 
enable  one  to  interpret  an  improper  integral  in  a  suitable  way 
(for  the  purpose  it  seems  to  be  infallible 

Source  and  vortex _of  pulsating  strength  in  two  dimensions . 

The  two-dimensional  problem  can  be  formulated  analogously 
to  (13  9) ,  and  solutions  found  in  a  similar  manner.  The  funda¬ 
mental  singularities  will  now  be  of  the  form  log  r  j  y  „  \ 

Y  r  J  ^  ^  and  r  ii. '  -  1,  i,’  *  •  .  The  results 

are  given  in  the  paper  of  Thorne  [1953]  cited  earlier.  We 
shall  follow  a  different  method  here  in  order  to  illustrate  the 
use  of  complex  variables  to  solve  sue  problems. 

We  siall  consider  simultaneously  a  source  of  strength 
and  a  vortex  of  intensity  I"  at  the  point  C  st  ft  4  /  4.  $<.Q  . 

* 
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In  the  notation  used  at  the  end  of  section  11,  we  shall  be 
coking  for  a  function  l  f ,  ||  analytic  in  ±  and  of  the  J'ovrm 
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(13.22) 


where  {  and  _  have  no  singularities  in  the 
In  addition,  and  must  each  satisfy  the 
condition  (11.7)  which  we  write 

I  | *  -  -  o)  -v  < v  x - ; o){>  -  o ,  v-  c’Vj. 


lower  half-plane, 
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Condition  4)  of  (1.3,9)  will  be  taken  in  the  somewhat  stronger 
lr  orm  t 
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where  m  and  fA  ar.  given  constants, 
The  radiation  condition  becomes: 
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Following  a  method  apparently  originally  due  to  Levi-Civita 
[see  Tonolo,  1913),  but  used  frequently  by  Keldysh  [1935], 

Kochin  [e.g,,  1939],  Stoker  [1947],  Lewy  [1946]  and  others,  we 
Introduce  the  functions 
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Then  ( 13  23 )  becomes 
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and  (13  22)  becomes:  the  two  functions 
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and  are  both  regular  everywhere  in  the  lower  ha  1  f  - 

plane.  A  function  with  -  /  M  “  -  ■'  J1  may  be 

continued  into  the  upper  half-plane  by  defir.  ig  A  f  j  = 

‘i  ’  -  <  ,  ,  )  ^  ,  the  bar  indicating  complex  conjugate.  Since 

is  regular  in  the  lower  half-plane,  the  extended  function  will 
be  regular  in  the  whole  plane.  In  addition,  one  may  derive 
easily  from  ( 13  24)  that  l  A  t  J I  <  £  i  I  "*■  D  for  suffi¬ 
ciently  large  <  *  :  then,  from  the  regularity  of  ft i  , 

such  an  inequality  holds  in  the  whole  lower  half-plane  and 
hence  in  the  whole  plane  after  reflection.  It  then  rollows 
from  a  known  generalization  of  Liouville's  Theorem  [see 
Caratheodory,  Theory  of  functions  of  a  complex  variable,  vol. 

I,  <>  168,  Chelsea,  New  York,  1954]  that  Ai(%)  “  ^  j  > 

where  6.  and  [-  are  constants.  It  follows  from  (13.27)  that 
Ca,  and  /-  are  real.  The  differential  equation 
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One  may  set 


without  loss  of  generality 
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ally  this  provides  a  proof  of  the  theorem  of  Stoker  and 
Ke ldysh  me n t i one d  earl i e r  ( s h o r 1 1 y  a f t e r  13. 6 ) ) 
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Since 
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satisfies  (13.27),  and  the  four  singular  terms  taken 


together  also  have  vanishing  imaginary  part  for 
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same  must  hold  for 
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Hence  A 
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must  have  the  same 


form  as 
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Substituting  the  resulting  expression  for 
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in  (13.26),  one  has  a  differential  equation  for 


The  solution  is 
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lower  half-plane  As  in  the  case  *  ,  we  may  set  -  - 

*  and  -  ,  must  now  be  chosen  to  satisfy  (13.25)  Making 
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One  may  easily  verify  that  this  choice  of  C 
produce  outgoing  waves. 

If  one  makes  the  change  of  variable  v/---  '  •  ‘  r~-c  '  in 

the  int tgral  term  In  f,  and  deforms  the  resulcing  path  to 
v*  ,  one  finds 
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Substituting  this  in  the  expression  for  f,  ,  one  finally 
obtains 
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Singularities  of  higher  order  may  be  found  by  different! 


ating  (13,28)  with  respect  to  y 


The  expression  foi 


]’  , 


«*  « 


may  be  put  into  a  somewhat  different  form  by  using 


r 


■“  *  «* 

. 


f 


-*i 

4  »U‘ 

i-tliw 


03 


"i  >  'i 


C 


o 


4^ 


\ 


O 


\ 


Then 


IV 

t  ' 


4  h 


»h  wi  wtnwwtiaww»t 


mH. 


.** 

< 


dF 


i  -  *  v 

. .  ' 


X  | 


P  / 


Vf1  4  V 


-jhr 


mrnmmmmmmm w—*1 

V* 


/  » 


i*HS, 


aw 


A  V  t  , 

•  '  *  rc  ""  "  /,  r  ci 


""1  V  \  l* 

w 


(13.29) 


/  /... 


One  may  continue  differentiating,  using  either  form  for  t  > , T •'  . 
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symmetrical  from  unsymme trica 1  solutions)  and  taking  the  appro 


priate  real  or  imaginary  part,  on©  finds  the  following  formulas 
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In  the  formula  for  the  logarithmic  singularity  may  be 

eliminated  and  the  coefficient  of  coi  crt  written  as 

(v.  John,  1950,  p.  100 j: 
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For  wat€»r  of  finite  depth  the  method  used  above  does  not 
work  as  conveniently  because  of  the  difficulty  of  formulating 
the  boundary  condition  on  uie  bottom,  Iv  v  It)  =  0  ,  in 

terms  of  the  function  A  {'£) .  However,  it  can  be  done,  yielding 
a  differential-difference  equation  for  {{£)  which  can  Us  solved 
by  use  of  Laplace  transforms  [cf.  Bochner,  Vorlesungen  tlber 
Fouriersche  Integrale,  Leipzig,  1932,  pp,  167-168).  The  method 
used  for  the  three-dimensional  problem  can  also  be  carried 
through  (see  Haskind  [1942b],  John  [1950],  and  Thorne  [1953]). 
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It  is  convenient  to  separate  the  vortex  from  the  source 
The  result i ng  func  ti on  s  a  re  a s  t o 1 1 ows . 
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The  remark  following  (13.18) 


concerning  the  form  of  the  last  term  of  that  formula  applies 


also  here.  The  rea]  part  of  either  of  these  gives  the  corres 


ponding  potential  fuact* 

For  the  source,  the  Integral  representation  and  the  series 


representation  analogous  to  (13,19)  are 


■i' 

'■< 


Thorne  [1953]  gives  the  potential  functions  for  the  higher-order 
singularities  and  the  function  for  the  logarithmic  singularity 
in  a  form  involving  and  '  and  hence  more  analogous  to  the 
one  i’  (13.31).  Voitsenya  [1958]  has  derived  the  complex  |*)ten- 
tial  or  a  source-vortex  situated  in  an  infinitely  deep  fluid 
of  density  £  lying  beneach  another  of  density  >  ,  , 

and  of  thickness  a  . 

Source  of  constant  strength  tn  uniform  motion:  three 
dimensions . 

IVe  shall  assume  the  source  moving  in  the  direction  0* 
with  constant  velocity  u0  .  Let  <  ■  %  \  be  coordinates 

in  a  system  moving  with  velocity  ^  „  in  direction  -.;x  and 
let  the  source  be  at  «.  i- ,  c 1  ,  <■  L  .  Then,  from  section 

11,  we  wish  to  find  a  function  Cp  ( A ,  i^,  satisfying 
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is  harmonic  in  the  region 
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For  fluid  of  finite  depth  I 


1"  1 


4)  is  rep1 aced  by  4  ) 


V|iJ.  Without  condition  5) ,  demanding  vanishing  of  the 
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motion  far  ahead  of  the  source,  the  solution  A?ould  not  be  unique 
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The  profile  of  the  free  surface  is  obtained  from  ntx,l) 
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Strictly  speaking,  the  solution  of  (13.30  will 


represent  a  sink,  i 


»■  * 


a  source  of  strength  -1.  However,  we 


shall  continue  to  call  such  solutions  sources 


A  solution  to  this  problem  may  be  obtained  by  methods  verj 


similar  to  those  used  for  the  source  of  pulsating  strength 


The  details  will  not  be  repeated,  but  can  be  found  in  Havelock 


[1932],  Sretenskii  [1937],  Kochin  [1937],  Lunde  [1951],  Peters 


¥ 


md  Stoker  [1957],  T tinman  and  Vossers  [  1955]  and  elsewhere. 


The  result  is 


-5 


* 

/> 


it*  -r 


«*  V 


A 

i-* 


if  ( jj>  ^ 

°  ■'+*'** 

_  -t-  ■  !  »  =  ^ 


; 


’  c 


-f-  ro5l"V"“V 


,1  v(}t4jsgc.\? 

i/\  o  5,V[v(x-ajs«c 


,  A  ^  ■  i 

( 13 .36) 


iiiIiIm  iji  -  i  i  r  i  *iiii  • 


ll  "iHtiH 


. I  lllhlirf  m 


IUVMI  l.iUlll 


I 


whts  re 


r’i  «-o  *• 


v 


'  *  -  a. 


■ff 


u 


1l» 

o  • 


The  potential  functions  for  higher-order  singularities  are 
unwieldy  and  will  not  be  given.  The  one  corresponding  to 
r  >  can  be  easily  obtained  by  ^-fold  differ¬ 
entiation  with  respect  to  ,  i f  one  remembers  that 


The  dipole  with  axis  in  the  direction  /  is  obtained  by 

differentiating  ( 13.30.  with  respect  to  X  and  will  be  used 
1  a  te  r . 

The  velocity  potential  for  a  source  moving  in  fluid  of 
finite  depth  has  been  calculated  by  Sretenskli  [1937]  and  by 
Haskind  [1945b]  The  form  given  bo  low  is  essentially  that 
given  by  Lunue  [1951] 
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where  K  „  -  r  •> .  is  the  real  positive  root  of 
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q:  j-i~  and  \{J  -*  J  l  s>  go 

In  the  double  integral  the  principal  value  is  necessary  only 
for  ,) j  «  J)  <  r~  ,  for  the  singularity  does  not  occur 

in  the  denominator  for  ;  i  -/  <  v  .  The  part  of  the 

double  integral  with  -3.  i  •/<*->  „  approaches  zero  as 

>.  — >  ±  co  ,  so  that  no  correction  is  necessary  in  order 

to  satisfy  condition  5).  This  is  the  explanation  of  the  lower 
limit  ’-'o  in  the  second  integral.  In  this  integral  the  deno¬ 
minator  vanishes  onLy  at  ,J  z  .  One  may  verify  that 

the  integral  is  convergent  by  noting  that 

.  ,  A.**'  2'> 

--  — - T— 

tv ;  ^  -  v  ^  5  e  c  *  o  i 

and  rewriting  it  as  an  integral  with  respect  to  ^ „  .  When 
“*  «o  ,  (13.37)  reduces  to  a  form  of  (13.36)  in  which  t , 
is  absorbed  into  the  double  integral. 

For  the  stationary  pulsating  source  the  asymptotic  form 
of  the  velocity  potential  for  large  fX  was  found  in  the  course 
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of  deriving  the  potential  function.  For  the  moving  source  of 
constant  strength  the  asymptotic  form  is  more  difficult,  to 
compute.  Since  the  form  of  the  free  surface,  r,  -;jU4  ^(rf,  3-4 
is  of  principal  physical  interest,  we  shall  discuss  the  asymp¬ 
totic  form  of  c# ,  Instead  of  -v  . 

Introduce  cylindrical  coordinates  >  j  1 
into  the  *  derivative  of  (13.36): 
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For  large  K  the  first  two  terms  taken  together  are 
0'^  )  •  Apply  the  theorem  (13,16)  to  the  integral  with 

respect  to  A  .  This  gives,  after  combining  with  the  second 
i n  rul ( 
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Since  is  symmetric  in  a  ,  we  consider  only 
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Consider  the  two  Jitegrals  containing  c^j  arid  let 

^  $  ec.1-  />  2  a  S  f  *  =*  \ 

Then  for  J  #  >  T  the  integral  becomes 

-  4  v%  . _i_i - -  .  j  ;  v  PA  A  /  ■ 

(3  ^  f  i  ^  +  o, )  -  ,  2  ^  c 

If  ,  the  lower  limit  is  cm  , 

In  either  case  one  may  show  that  the  coefficient  of  £os  vf? 
is  single-valued,  continuous,  absolutely  integrable  and  mono- 
tonically  decreasing  as  a  function  of  A  .  b.’  integration 


by  parts  one  may  then  establish  the  following  estimates  as 
[cf.  S.  Bochner,  Vorlesungen  uber  Fouriersche  Integrate, 
Leipzig,  1032,  i>  3}  : 
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If  o<  s  T7  ,  the  two  integrals  in  (13.40)  combine  to 

give 
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Consider  now  the  remaining  integral  it  (13.40),  and  let 
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The  integral  taites  the  form 
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an  equation  which  has  real  roots  when  <  |  /  P  t  i.e. 

when 

|  y  /.  <  -<  <  '  f'S 

When  -JT  <;  at  <  *r  -  /  -esi*  -  r* ,  ,  the  Fourier-integral  esti¬ 

mate  used  for  the  other  integral  may  be  applied  to  give 

f  i  *  I' V  ^  £  5  I  O'  •  . 

- - - -  ^  (  <"  • 

<  *.✓  <> 

When  <  »  <  t  ,  ^  is  a  two-branched 

function  of  m,  and  the  resulting  two  integrals  each  have 

singularities  at  one  of  the  limits.  Thus  tne  elementary  method 

of  am  lysis  used  above  can  no  longer  be  applied.  However,  a 

modification  of  the  method  above  can  be  carried  through  [see, 

0 

e.g..  Erdelyi,  Asymptotic  expansions,  Dover,  New  York,  1956, 
pp.  46-56];  the  classical  treatment  is  by  the  method  of  sta¬ 
tionary  phase  which  is  well  discussed  in  Stoker  [1957,  ch.  8]. 

The  estimates  already  derived  are  of  the  same  order  as 
the  remainder  term  in  (13.39).  Analysis  of  this  term  produces 
terms  Mch  cancel  the  terms  in  !/Rsi*»w  already  derived 
thus  oving  an  apparent  singular  behavior  near  the  X-  axis. 

The  asymptotic  form  for  the  surface  yj  ( T?  ,  oi) 
a  source  of  strength  W*  (i.e.  -**  /  r  )  is  as  follows; 
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In  order  co  have  some  idea  of  the  font  of  the  free  surface  far 
behind  the  source,  one  may  graph  the  curves 

V  k  /M,  (  ^  )  ~  ^  a  -  -  «-  >  n  ,  ?*  t\  Mi  *)  -t  -  77  r  -L  c  7r  ^  >  C  , 

showing  the  traces  of  the  wave  crests  in  the  region 

TT  -  h  *  (  i  s  —  c  <*  .  T  +*  *  •  '  ■**  - 

This  gives  the  well  known  pattern  she  n  in  Figure  la.  The 
first  equation  gives  the  transverse  waves,  the  second  one  the 
diverging  waves.  The  wave  length  along  is  i~  fv  and 

along  the  boundary  lines  s  .  The  expansion  is  not 

suitable  in  the  region  near  the  boundary  lines  y  :  *rt  . 

f  -s  ,  */n 

As  '!  -ii  ,  f-T  ,  the  term  1 1  -  1  ^  j  — '  O  and  tht 

amplitudes  become  infinite.  A  special  investigation  of  the 

-  f/ 1 

region  near  is  necessary  and  shows  /vf<  as 

leading  term;  ^  may  be  expressed  in  terms  of  Airy  functions 
(cf.  Ursell  [1959]). 

Eisentially  the  same  pattern  is  produced  by  a  moving  con¬ 
centrated  pressure  on  the  free  surface;  it  was  first  analysed 
by  Kelvin  [1906  -  Papers,  vol.  IV,  pp.  407-413].  The  asymptotic 
behavior  for  moving  pressure  distributions  has  been  extensively 
studied  (e.g.,  Hogner  [1923];  Teturo  Inui  [1936];  Peters  [1949]; 
Bartels  and  Downing  [1955]).  L  amb  [1926]  has  given  the 
asymptotic  form  of  the  surface  over  a  moving  submerged  dipole. 
The  form  of  the  surface  near  the  moving  dipole  has  been 
investigated  by  Havelock  [1928],  who  gives  traces  of  the  profile 
on  planes  x  -  const,  for  several  values  of  <x  between  j  7r 
and  ft  (the  radial  lines  of  Figures  lb  and  lc)  and  for 
v |  =  i/i  and  4.  Havelock’s  computations  were  later  used 
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by  Wig ley  [1930]  to  produce  the  contour  curves  shown  in 
Figures  lb  and  1c , 

A  similar  analysis  may  be  nr  da  tor  (13.37),  a  source 
moving  in  fluid  of  finite  depth.  For  a  moving  pressure 
distribution  the  problem  has  been  treated  by  Havelock  [ 1908] 


and  Teturo  Inui  [1936].  The  pattern  is  modified  as  follows 


If 


IS  r 


>  (  ,  the  pattern  is  qualitatively  like  that  for 


*  i  i* 


However,  the  wedge  within  which  the  disturbance  is  chiefly 


contained  has  a  wider  aperture  and  as 


the  aperture 


approaches 


ft 


radians  on  each  side  of  the  line  of  motion 


In  addition,  the  wave  length  of  the  transverse  wave  system 


Increases  and  approaches  infinity  as 


i- 


'-/hen 
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the  transverse  wave  systeu  is  missing  completely,  but  diverging 


waves  still  occur  In  a  wedge  of  aperture  varying  from 


t  o  ... 


as 


h  -  ->  < 


f 


(See  also  Ekman  [1906],  who  has  considered  the 


free  surface  over  a  dipole  on  a  flat  bottom.) 

Figure  2  from  Havelock  [1908]  shows  the  half-angle  of 
the  aperture  . 

Koch  in  [  1938c]  has  gone  furt*°r  in  this  type  of  problem 


He  has  derived  the  potential  for  s~ 


e  situated  in  a  fluid 


of  density 


and  depth 


> 


bounded  below  by  a  plane,  over 


which  is  lying  another  fluid  of  density 


.r  L»' 
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extending 


infinitely  far  upwards.  The  lower  fluid  moves  with  velocity 


/■•I 
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the  upper  with  velocity 


i 


in  the  same  direction.  He 


also  finds  the  asymptotic  behavior  of  the  solution 
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Singularities  of  constant  strength  In  uniform  motion: 
two  dimensions. 

For  submerged  sources  and  vortices  in  two-dimensional 
motion  the  complex-variables  method  used  for  the  pulsating 
source  may  again  be  applied.  For  infinitely  deep  fluid,  the 
computation  has  been  carried  out  in  this  way  by  Keldysh  and 
Lavrent'ev  [1937]  and  Kochin  [1937];  a  detailed  exposition  is 
given  in  the  textbook  of  Kochin,  Kibr.l’  and  Roze  [1948, 
ch.  VIII,  19]  .  Havelock  [  1927]  and  Sretenskii  [1938]  have 
treated  the  problem  by  different  methods.  The  complex  velocity 
pote  itial  for  a  combined  source  of  strength  S  and  vortex  of 
intensity  f  at  c  --  ci-m  b-  ,  4-  <  O  ,  is  given  by 
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The  real  velocity  potential  for,  sav,  a  submerged  source  can  be 
obtained  from  any  ol  these  equations.  The  last  one  gives  a 
form  analogous  to  (13.3G) 
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Higher-order  singularities  can  be  obtained  by  differentiating 
(13.43).  The  complex  velocity  potential  for  a  dipole  of  moment 


t*  and  axis  in  the  direction  e 
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is  given  by 
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If  in  the  last  term  of  the  first  equation  of  either  (13.411) 


<$  S 


or  (13.45),  one  makes  use  of  the  identity 
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it  ia  not  difficult  to  see  that  this  last  term  is  equivalent 
to  a  distribution  of  dipoles  on  the  ray  from  c  to  -  ©o 
parallel  to  the  x  -axis.  The  moment  density  and  axis  can  be 
determined  for  the  three  cases,  source,  vortex  and  dipole,  by 
comparison  of  the  integrand  with  the  first  term  of  (13.45). 

For  the  case  of  finite  depth  the  complex  velocity  potential 
has  been  calculated  by  Tikhonov  [1940]  and  is  also  given  by 
Haskind  [1945a]  for  both  source  and  vortex.  We  give  separately 
source,  vortex  and  dipole: 
source : 
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vortex : 
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(13.48) 


Here  C ^  =  3.  ~  '  &  4  2  '  ^  and  the  inst  summand  in  each  of 

it 

(13.46)  to  (13.48)  Is  to  be  deleted  if  W,s^/cv  i  1  ; 
is  the  positive  real  root  of  "V  -$1)  -  UUo  ,  which 
exists  only  if  v  n  >  1 

Asymptotic  form  of  these  functions  as  >  -*  -ao  is  easily 
seen  to  be  given  by  double  the  last  term  in  each  expression. 
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When  “uV  ^  j  ,  the  disturbance  .  s  only  local,  a  fact  which  cor¬ 
responds  to  the  absence  of  transverse  waves  behind  the  three- 
dimensional  source  for  >  n  <  1  . 

Kochin  ( 1937a, b]  has  derived  the  complex  velocity  potential 
when  fluid  of  density  overlies  the  fluid  of  density  f, 

containing  the  singularity.  The  lower  fluid  may  be  of  infinite 
or  finite  depth;  the  upper  one  is  taken  infinitely  deep.  Their 
velocities  may  bi  different. 

Source  of  variable  strength,  starting  from  rest  and 
following  an  arbitrary  path. 

Consider  now  a  source  whose  position  and  strength  at  time 
f  >o  are  given  by  (■'■>  ft  ^  'ft  c't))  and  ;t)  t  where  . 

Let  vv  ti  ; ■  *  - -  t  ■  -  .  The  conditions  to  be  satisfied  by  the 

velocity  potential  $  1  t)  are 

J  c ' 
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If  one  assumes  a  solution  in  the  form 
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where  i", '  .  •  *  (i*  .*  '*  ,  then  j.  must  be  h 

in  -  v  3  and  satisfy  4>,  5),  6)  and 


a  nu on  it 
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It  follows  from  the  conditions  that,  for  •” ^ /  «f.  -  constant, 


which  we  may  take  as  zero.  Let  be  the  Laplace  tr 

form  of 
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Then  is  a  harmonic  function  «  n  y<3  satisfying  4 

and  5)  for  each  /£*  and  also,  after  T.akine  use  cf  S)  the 
condition 
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Since 
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is  a  harmonic  function  in  <  o  vanishing  on  s  O  and  at 
infinity,  it  is  identically  zero.  Making  use  of  (13.12;  differ¬ 
entiated  with  respect  to  y  and  changing  the  order  of  integration 


one  obtains 
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Making  use  of  the  convolution  theorem  and  the  fact  that  (V*  a 
Is  the  transform  of  ^  f  ^ *  r  1  t  ,  one  may  find  the 
original  function  'P ,  (  <  ^  t)  : 
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-  2  f  (J*)*  [  it 


J6  /  *  Wt.x 


For  fixed  t  one  may  easily  verify,  using  known  properties  of 
the  Fourier-Besso 1  transform  [cf.  Wfctson,  Bessel  functions,  f 
14.41],  that  $,  Is  5  '/ *i,  and  hence  that  5)  is  satisfied. 


One  has  then  the  result 
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(13.49) 
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By  a  more  refined  analysis  of  the  behavior  for  large  \  [cf, 
Stoker,  1957,  pp  190  1]  one  may  establish  tint  is 
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and 


are 


i  )  as 


For  some  time  t  >  t  . 


,  one  may  write  [  in  the  form 
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Here  the  first  summand  ,  represents  the  effect  at  time  '  " 

of  the  action  of  the  source  from  A  -  to  .  The  remaining 

terms,  .[  ,  are  the  same  as  (13.49)  with  T  measured  from 

l \  ■  :-tt  i.y;  l,:  )  ,  and  show  the  effect  at  time  r  of 

the  action  of  the  source  from  t  "■  f  j  t <■  /  T  •  (This  is,  of 

course,  what  one  would  expect  from  the  linearity  of  the  problem 
and  the  fact  that  the  choice  of  /  -  l-~>  is  arbitrary.)  When 
t  =  t «  ,  [■  i  reduces  to 

*v  ( t o)  _  rA  1  ft)..  . 
hits)  r<  fto) 

Thus 

W)  =o. 

Thin  fact  provides  a  basis  for  Havelock's  procedure  in  similar 
problems,  a  procedure  originating  with  Kelvin  in  the  treatment 
of  moving  pressure  distributions.  The  idea  is  roughly  as 

,i  i  i 

*■  * 
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follows.  Divide  the  path  of  the  source  into  sr a  1 1  segments  of 
time  span  At  .  If  .If  is  small  enough,  the  effect  of  gra¬ 
vity  upon  the  fluid  motion  produced  by  the  source  during  this 
time  interval  will  be  negligible,  and  one  may  take  the  boundary 
condition  at  the  free  surface  to  be  ,  -  ~  The  distortion 
of  the  surface  resulting  from  the  action  of  the  source  during 
this  short  interval  is  found  and  the  future  behavior  of  the 
distortion  computed  while  taking  account  or  gravity.  Summing 
over  all  At  and  taking  the  limit  leads  to  the  potential 
function . 


The  expression  (13.49)  has  been  essentially  given  by 
Haskind  [ 1946b]  and  Brard  [  1948a  ]  .  Special  choices  of  jrj 
and  of  the  motion  of  the  source  lead  to  cases  similar  to  those 
treated  earlier.  Thus,  if  t )  -  .  ,  r  '  and  (a,  4,  i)  is 

fixed,  one  has  the  potential  function  for  a  stationary  source 
of  oscillating  strength,  starting  to  oscillate  at  J  . 
Carrying  out  the  t  integration  and  taking  a  limit  by  using, 
say,  the  Fourier  Integral  Theorem  (13.1  allow  one  to  derive 
(13.17).  The  radiation  condition  is  automatically  satisfied. 
The  velocity  potential  for  finite  values  of  t  may  be  written 
in  the  form 
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The  leading  terr.  n  the  asymptotic  expansion  of  the  last  summand 
gives  the  last  summand  of  (13.17), 


If  one  takes  >vf)  =  •  <*v  -  a  constant,  i  -Ji  Jjlltjf  , 

-  ts  iit  -  f,  ,  one  obtains  the  velocity  potential  for  a 
source  suddenly  brought  into  existence  at  t  ~  Z  and  moving  with 
constant  velocity  in  the  direction  C/*  [cf.  Lunde,  1951,  p. 

18].  A  limit  as  t  -i »  "o  will  give  (13.36),  the  proper  boundary 
conditions  at  infinite  being  again  automatically  satisfied. 

For  finite  {  the  velocity  potential  in  a  coordinate  system 
moving  with  velocity  in  direction  D>  (  k  ~  >-  -  soUat 

y  is  «iven  by 
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(13.51) 
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The  two  cases  just  discussed  may  be  combined  by  choosing 
1V\  |  (  ]  =  <v  t<jr.  rt  and  Ct(f)  z  4.  d*t  ,  rti2  i  'ti  -  C0  .  The 

modification  of  (13.51)  is  simple:  a  factor  c  js  <rt  must  be 
put  with  the  first  two  terms  and  a  factor  <-qs  put 

at  the  end  of  the  integral.  The  asymptotic  form  as  T ->>  03 
can  again  be  found  by  uso  of  the  Fourier  Integral  Theorems 
(13.16)  or  simple  modifications.  However,  if  the  resulting 
formula  is  written  out  as  principal-value  integrals  plus  other 
terms,  the  expression  is  very  unwieldy;  it  may  be  found  in 
Havelock  [1958].  Use  of  complex  integrals  allows  one  to  com¬ 
press  the  formula.  Let 
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Motion  of  a  source  on  a  circular  path  of  radius  D  may 
be  treated  by  taking  -  —  t  s  -  ^  t  £n  (13.49). 

For  constant  this  problem  has  been  considered  by  Sretenskil 
[ 1946a, b;  1  937],  Havelock  [1959],  and  Stoker  (1957). 

One  may  derive  a  formula  tnalogous  to  (13.49)  when  the 
source  moves  in  the  presence  of  both  a  horizontal  bottom  at 
j  -  i  and  a  free  surface.  The  derivation  may  be  carried 
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cf.  Lunde,  3951,  p. 
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where 

-  (  *  '  !  t  )J  -4-  1  -j  V  }  ^  i- :  ft  j  -h  "  C  1  M) 

Two-dimensional  formulas  corresponding  to  (13.49)  and 
(13.53)  may  also  be  derived.  They  are  as  follows,  with  the 
source  and  v  rtex  separated  for  finite  depth: 

infinite  depth:  (13.54) 
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depth  h,  soui<e:  (13.55) 
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depth  h,  vortex: 


(13.56) 
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Higher-order  singularities  may  be  generated  by  taking  derivatives 
with  respect  to  ..  .  One  may  transfer  to  moving  coordinates, 
etc.,  just  as  in  the  three-dimensional  case  (see  Havelock  [1949] 
for  (1?.54)  in  moving  coordinates).  The  velocity  potential  for 
a  steadily  moving  source  of  pulsing  strength  in  two  dimensions 
has  been  given  by  Haskind  [1954,  p.  23  f£.),  who  also  gives  the 
asymptotic  expressions  for  large  values  of  J  3  .  When 
iv  <  '/if  ,  there  exist  one  wave  far  ahead  of  the  moving  source 
propagating  in  the  same  direction  and  three  far  behind,  one 
propagating  in  the  same  direction  and  two  in  the  opposite  direc¬ 
tion;  when  t  >  ,/4  ,  there  exist  two  waves  far  behind  propagating 
in  the  opposite  direction.  The  analysis  for  finite  depth  has 
been  given  by  Becker  [1956]. 


3.,;  !?■,  .  . ijj'H.  '  ;ifM niriirinf iniiMWffjflwiniyuiiiwii 


1* *  flo-M  simple  Physical  .solutions.. 

In  this  section  we  consider  periodic  waves  In  an  ocean  of 
Infinite  horizontal  extent,  either  Infinitely  deep  or  with  a  hori¬ 
zontal  bottom.  In  canals,  and  at  an  Interface.  The  linearizing 
parameter  £  of  section  10 a-  may  be  taken  to  be  the  ratio  of 

amplitude  to  wave  length, 

14  of  .  Standing  waves  In  an  Infinite  ocean. 

It  Is  appropriate  to  the  physical  problem  to  require  that  the 
motion  remain  bounded  everywhere , 

Consider  first  two-dimensional  motion.  Then,  from 
13  n  ,  the  only  solutions  of  the  form  =  <p  cos(  cr  fc  +  t)  are  given 
by 

cod  (r+  (14,i) 

for  Infinite  depth,  and 

■-  wc(^k)cos(w,X4o.)cas(<rNt:),  hA  "J,  -v  =  0,  (14.2) 
for  finite  depth. 

The  correspond— ig  forms  of  the  free  surface  are  given  by 

M  (*,f)  =  A  c<rx(v*4*) 

and 

rj  |Xlt)  ^  4  cc5  fwi.X-nW"  (crt  +  b); 

respectively.  These  represent  standing  waves  according  to  our  defi¬ 
nition  in  section  7  We  recall  that  r**  >v 

It  la  of  Interest  to  examine  the  streamlines  and  the  paths  of 
the  Individual  fluid  particles.  The  streamlines  of  the  motion  can 
be  'laslly  found  from 
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The  streamlines  are  then 

v  ?  ‘  4  >f  t  , 

a  |  *  ■  <i  f  v  /  t  <*t  i  r  : 


(14.3) 


Sir»U  |sin  (w0«-»fc)J  =  smfc  w#(y^+J,)y  0 

for  infinite  and  finite  depth  respectively;  here  is  the  low¬ 
est  point  of  the  streamline.  If  the  fluid  Is  infinitely  deep, 
the  streamlines  are  all  congruent.  Figure  4  shows  three  of  them 
for  a  quarte  r  wave  length  and  0  -Ot  v  =  4  .  The  vertical  line 
is  also  a  streamline.  If  the  fluid  is  of  finite  depth,  the 
streamlines  vary  with  depth.  Figure  5  shows  streamlines  cor¬ 
responding  to  -  0J  -  0Sy  -  0. 9  for  &-0  ,  U  -  i  ,  ~  i- 
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The  horizontal  line  y  *  -1  and  the  vertical  line  x  «  0  are  also 
streamlines. 


Since  the  streamlines  are  time -Independent,  they  also  contain 
the  curves  for  the  trajectories  of  Individual  particles.  However, 
the  trajectory  of  an  Individual  particle  will  be  an  oscillating 
motion  of  small  amplitude  along  a  segment  of  the  streamline  passing 
through  the  point.  Thus,  In  Figure  5  the  particles  on  the  bottom 
simply  oscillate  back  and  forth  about  an  equilibrium  position,  those 
directly  beneath  a  crest,  l.e,  at  x  •  0,  oscillate  vertically,  etc. 
In  view  of  the  infinitesimal  wave  approximation  used  In  this  chapter 
the  streamlines  have  physical  significance  for  only  a  small  distance 
above  the  equilibrium  free  surface,  y  •  0. 

In  order  to  Investigate,  at  least  approximately,  the  behavior 
of  the  trajectories  more  fully,  we  may  replace  the  actual  trajectory 
by  it-  tangent  at  an  average  position,  say  (xQ,y0),  an  approximation 
consistent  with  the  assumptions  made  In  linearizing.  Then  the  equa¬ 
tions  describing  the  trajectory  become  patting  «.  ■  x  -  0) 
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for  infinite  depth,  and 
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for  finite  depth.  The  approximate  trajectories  are  then 
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for  finite  depth.  For  Infinite  depth,  the  amplitude  of  oscillation 
drops  off  very  rapidly  as  depth  or  the  fcoulllbrium  position  Increases, 
the  ratio  of  the  ampllt  ide  at  depth  yQ  to  the  amplitude  at  the  free 

VJ, 

surface  being  <2  .  The  same  ratio  for  the  case  of  finite  depth 

Is 

Smh  fV*. «  ( S  i  * 

Siwl/lw.h  -4  S X% 

Thus,  on  the  bottom,  when  yQ  -  -h,  the  amplitude  is  zero  under  the 
crests  and  maximum  under  the  nodes.  As  is  evident  from  the  equations 
of  the  trajectories,  the  path  lines  of  particles  on  the  free  surface 
are  approximately  as  In  Figure  6. 


Figure  6 

In  order  to  explain  the  apparently  Inconsistent  behavior  at  the  nodes 
one  must  go  to  a  higher  approximation  than  the  linearized  theory  used 
in  this  chapter. 

Let  us  now  consider  throe-dimensional  solutions.  The  standing- 
wave  solutions  are  of  the  form 


or 


for  finite  depth. 


for  flnlte  iepth, 

where  ^(x,z)  Is  a  solution  of 

^  x  +v> =o  °r  a2.x+kx-o> 

respectively,  regular  everywhere  In  y  <  0. 
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Two  particular  cases  are  of  especial  Interest.  The  first 
corresponds  to  separation  of  variables  in  rectangular  coordinates 
[see  (13.5)  and  (13.6)].  Ttie  solutions  are 


I, x+^.s/^+r)  (14.6) 

4,V+ V 

for  infinite  depth,  and 

=  *  c#^  c«s  ft***) 

Vt  l?vV-  1^*V  J  **o  -  v  S.QJ  (14.7) 


for  finite  depth.  The  other  solutions  result  from  separating  variables 
in  polar  coordinates  £  see  (13.7)  and  (13.8)1.  They  are 

for  Infinite  depth,  and 

*  cc^  M,^  +  ^Owf^0!?)tos(hK4ricoi(»‘t+^  arO,  1,  ,  (14.9) 


for  finite  depth.  The  form  of  the  free  surface  may  be  found  Imme¬ 
diately  from  *j  (x , z , t )  =  -$t(x,o,z,t)/g.  These  are  all  standing 
waves . 

The  streamlines  and  path  lines  may  be  found  for  these  two 
cases  with  no  special  difficulty.  For  the  first  case  for  finite 
depth  the  streamlines  ar'-  the  Intersections  of  the  surfaces 

=  C, 

|  S  lw  ^  ^  |  M0  —  C-i  . 


i 


||j||||||jy||||,|| 


(14.10) 


The  vertical  lines,  x  *■  p -rr  /  ^  ,  z  ■  JEW  w  ,  passing 

through  the  maxima  and  minima  are  streamlines.  The  points  on 
the  vertical  lines  x  -  ,  z  -  +  passing 

through  the  sattlepolncs  are  all  stagnation  points.  The  projection 
on  y  ■  0  of  the  field  of  streamlines  is  Indicated  qualitatively  by 
Figure  7. 


The  behavior  in  a  projection  on  a  vertical  plane  Is  similar  to  that 
for  two-dimensional  motion. 

In  the  second  case  above  one  may  easily  visualize  the  stream¬ 
lines  for  the  case  of  pure  ring  waves,  n  »  0.  For  finite  depth 
they  are  given  in  a  plane  a  »  const,  by 

J,  •*,  (pi *)  =•  c,  (i4.ii) 

together  with  the  vertical  lines  at  the  zeros  of  J  •  The 

behavior  of  the  curves  is  qualitatively  similar  to  that  of  the  two- 
dimensional  case. 

In  both  cases  approximations  to  the  path  lines  can  be  found 
as  In  the  two-dimensional  case. 

14  .  Progressive  waves  in  an  infinite  ocean. 

By  taking  the  proper  linear  combinations  of  the  standing-wave 
solutions  one  may  obtain  progressive  waves.  Thus,  adding 

$  “  d  ^V^"co5  V*  C  CS  Tt  •  A  £,  ,  w  v*  S'*'  f 
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one  obtains 


4 


'  i  c  3  t  * 


■t) 


(14.12) 


which  represents  a  progressive  wave  moving  to  the  right  wltn  velocity 


(14.13) 


where  ),  -  2 tt /v  is  the  wave  length.  Subtracting  yields  a  pro¬ 
gressive  wave  moving  to  the  left.  If  one  takes  the  coefficient  in 
as  Aj  and  In  ^  as  0^ ,  the  sum  may  be  written 


ey*  [fa**  i)  «cs/W-rt)  -t  (o.  -ax)co: 


1/X  +  T 


t) 


This  Is  a  superposition  of  two  progressive  waves  of  different 
amplitudes,  one  moving  to  the  left  and  one  to  the  right.  If  a,  ■  a„, 
a  pure  progressive  wave  is  obtained,  if  a2  -  O^one  obtains  again  a 
standing  wave,  as  a  superposition  of  two  progressive  waves  moving 
in  opposite  directions. 

For  water  of  finite  depth  h  the  corresponding  expressions  for 
4*  may  be  obtained  by  replacing  by  cosh  m  (y  +  h)  and  v 

by  mQ.  The  phase  velocity  is  given  by 


(14.14) 


As  h  the  velocity  approaches  that  obtained  ebove  for  deep  water. 

In  fact,  if  h/X  >  0.2,  the  velocity  is  already  within  0.1  of  the 
value  for  deep  water.  C  is  an  increasing  function  cf  )  ,  but 

cannot  increase  indefinitely  as  in  the  case  of  infinitely  deep  water, 
for  (14.14)  implies 

■  .  '  ,  I 

w  Mil  'lit  '  WW"  . . . , 
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(14.15) 


-Tha  streamlines  for  the  progressive  wave  moving  to  the  right  are 
given  by 


$;*k  **«'* 


_  r  t  ;l  -  U 


(14.16) 


for  infinite  and  finite  depth,  respectively.  At  a  given  instant 
t  these  hive  the  same  shape  relative  to  a  crest  as  the  streamlines 
for  a  standing  wave.  However,  since  they  are  time -dependent,  the 
path  lines  for  particles  do  not  lie  on  the  streamlines.  Hie  path 
lines  may  be  found  approximately  for  a  particle  with  equilibrium 
position  (xd-vo)  from  the  equations 


-  4 u  '/*. 


Hilo  approximation  Is  consistent  with  the  assumptions  made  In 
linearizing  the  boundary  condition,  as  can  be  seen  h-  assuming  a 
solution  in  the  form 

xft)--  Xs4fX,(t)+-  ,  -  (0  =  'Jo*  1  ‘J-  "j 

where  £=  ac*v/2.T|.  for  Infinite  depth  and  £z  aer/At/lufy 
for  finite  depth,  substituting  in  the  exact  path  equations,  and 
retaining  only  first-order  terms. 

For  infinite  depth  the  particle  trajectories  are  given  by 


-t  v$< 

xa  -  *vc-  e 


cosfl'Xi 


V*(vx0 


(14.17) 


.{  v*0 

The  particles  follow  circular  orbits  of  radius  flvcr  e  u  about 

v 

the  equilibrium  position  (x0,y0);  at  the  top  of  the  orbit  they  are 


I C  7 


moving  in  the  same  direction  as  the  wave.  The  orbital  velocity  is 

V  'i  m 

0.1/ &  1  ,  30  that  the  motion  dies  out  quickly  as  ly  increases 

for  example,  at  a  depth  of  one  wave  length  the  velocity  and  orbit 
radius  are  only  1/535  the  value  at  the  free  surface.  Although 
the  particles  at  the  crest  of  a  wave  are  moving  in  the  same  direction 
as  the  wave,  their  velocity  la  not  necessarily  *he  same  and  is,  in 
fact,  much  smaller  in  view  of  the  assumed  smallness  of  £  a  (a~y/c)(V/iir). 

For  finite  depth  the  orbits  are  elliptical  with  the  major  axis 
horizontal ; 


x  =  x0  “  G. 

ft  j  "  t  k  )  e  (*'  1 

J 

X  -  fl-t ), 

V-  > 

*w.  C*"1 

s '  K 

„  (  H .  »  -+  *  )  J  1  K  '  /> 

,  x  - 

The  particles  again  trace  the  orbit  in  a  clockwise  direction  except 
that  on  the  bottom  they  simply  oscillate  along  a  horizontal  segment. 
Figure  8  from  Ruellan  and  Wallet  [l950]  shows  the  path  lines  for 
a  variety  of  cases  of  superpo~ed  waves.  The  topmost  picture  shows 
the  orbits  for  a  pure  progressive  wave  moving  to  the  right.  The 
bottom  picture  is  a  superposition  of  progressive  waves  of  equal 
amplitudes  moving  in  opposite  directions,  i.e.  a  pure  standing  wave. 
The  intermediate  cases  show  superpositions  with  varying  ratios  of 
the  amplitudes.  The  Intermediate  cases  are  instructive  in  that  not 
only  path  lines,  but  also  streamlines  are  visible. 

Since  the  progressive-wave  solutions  are  steady  with  respect 
to  a  coordinate  system  moving  with  the  wave,  it  clear  that  we 
could  have  obtained  a  steady-state  solution  as  ^mall  motion 
superposed  upon  a  uniform  flow.  If  we  take  a  complex  velocity  poten¬ 
tial  in  the  form 


I 


fo? 


F(^  --  -  *  f  '  0 


(14.19) 


Then  (see  equation  11.6)  f  must  satisfy 


~?<2  *  %  *•  *  C‘-f  —C'  for  y  -  0 

*nd  el  ..her  1  f  *♦  0  as  y  — oo  or  1  ^  f  1  «  o  for  y  ■  —h .  THe 
solution  for  the  first  case,  infinite  depth,  is  given  by 


-,v> 


f  =  ae  -  At'**  i‘  ci  ^  U-  v:  5  /cv 

(14.20) 

The  solution  for  the  finite-depth  case  is  given  by 

f  -  £•  c  os  rvc  (t.  -t  >  k') 

-  (X  t  £ii  (0S  Y,  Y<  e(^  +  K)  $  r-  iv(>  5iV,ln  >V0(--n')3;  (14.21) 


where  mQ  must  satisfy 

C-"1-  Krt  o  -  ^  Fd  *  k  r/  B  l\  =  O 

The  same  relation  is  found  in  (14  14).  A  real  solution  does  not 
exist  if  c  /gh  >  1  and  in  this  case  there  is  no  wave-like  motion 
consistent  with  the  linearized  theory.  The  streamlines,  identical 
here  with  the  path  lines,  are  obtained  from 

-  c.  d  -t  -  O 

*  « 

One  may  replace  this  equation,  consistently  with  the  linearization 
assumptions  [  cf.  (lO.l^V  *  „ 

-  c  I-  •  V”  ^0)  =  ^ 

where  yQ  is  the  mean  height  of  the  streamline.  Thus,  for  finite 
depth,  they  are  given  by 


Ic 


^  —  ~  c  'J  A  ^  ^  c  }c  +  k  1  1  '*  trio  X  j  ( 14  , 2.i. 

an  easily  constructed  family  of  curves  .  In  the  foregoing  •  -e  have 
tacitly  taken  X  to  be  real.  However,  It  may  be  complex  and  thus 
Include  waves  of  different  phase 

We  note  finally  that  (14.8)  or  (14.9)  allow  one  to  construct 
waves  progressing  like  the  apexes  of  a  wheel.  However,  outwardly 
progressing  waves  can  be  constructed  only  when  the  solution  invol¬ 
ving  Yn  Is  used,  and  this  has  a  singularity  at  the  origin. 

14  ^  .  Periodic  waves  in  rectangular  canals . 

Let  us  suppose  that  the  fluid  is  contained  between  tn°  planes 
z  -  0  and  z  -  Ji  .  Then  the  velocity  potential  must  satisfy  the 
additiunal  conditions 

r  (14.23) 

This  condition  is  automatically  satisfied  by  the  two-dimensional 
waves  discussed  in  1^  ■>  so  that  they  present  no  special  interest 
here.  However,  condition  (14.4)  does  restrict  the  th* ee-dlmensional 
solutions  (14.6)  and  (14.7),  for  k 2  must  now  satisfy  (  aklng  /  «■  0). 


n  »  1,  2,  •  •• 


2  2  2  2 

Since  +  k2  «  y  or  niQ  ,  there  can  be  no  solution  periodic  in  x 
unless 


v  d. 

11  <  — 


.  *V!  s  t*. 

or  k\  <  — 1 
'V 


(14  24) 


respectively.  Hence,  for  frequencies  below  a  certain  critical 


frequency  O'  1 ,  where 


(14.25) 


for  Infinite  or  finite  depth  respectively,  there  can  exist  no  three- 
dimensional  standing  waves  in  a  canal . 

Let  us  form  a  three-dimensional  progressive  '«ave  in  a  canal 
of  finite  depth  by  adding  standing-wave  solutions: 

$(*}  ~  a  f0S^  =n-rr/j  . 

The  velocity  of  the  progressive  wave  is  given  by 


_Z.V  -  a  L  f  <  —  \ 

4(V  ~  3  ^  1  ’  **eVJV  J  V 


o 

a 


(14.26) 


As  in  iS  case  treated  above,  there  can  exist  no  three-dimensional 
progressive  waves  unless  T  >  Oj  .  However,  if  they  exist,  their 
velocity  is  higher  than  the  velocity  of  two-dimensional  waves  of  the 
same  frequency. 

One  may  define  similarly  a  sequence  of  critical  frequencies 
i.  O' g  j  *  •  • »  where 


when*  CTk  <  0"  <  k  tyres  of  three-dimensional  wave  a  are 

possible  with  n  -  1,  2, 


k. 
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14  5 .  Wave a _ at . jn . Interface . 

Let  us  now  suppose  that  two  fluids  are  present,  one  lying  over 
the  other.  Variables  referring  to  the  upper  and  lower  fluids  have 
subscripts  1  and  1  respectively.  From  (10.7)  and  (10.8)  the 
linearized  boundary  conditions  for  a  small  disturbance  are 

,  .  _  t  (14.27) 

P.  ['P.H  =  fJ£t.  +  ?■<#»,], 

both  equations  to  be  satisfied  at  the  equilibrium  position  of  the 
interface.  We  shall  consider  several  typical  problems. 

Let  the  upper  fluid  fill  the  region  y  ;>  0,  and  the  lower  fluid 
the  region  y  <  0.  We  require  of  a  solution  that 

j  jgfrdj  4i|h>o  as  y  -  *o  and  |  $  *'<*<*  4  a.  I  — 1  0  <3  y  -*  +  *°  • 

In  looking  for  a  standing-wave  solution,  one  may,  following  14 oc  , 
takf 

-  «,£  ;  ^  cos^  +  f)  ; 


where  the  relation  between  a.  and  a„  and  m  and  cr  is  to  be  deter- 

X  c 

mined  by  (14.27),  and  ip  satisfies 

tf  t  cf  =  0. 


The  first  equation  (14.27)  gives  immediately  that 


a1  +  iig  »  0. 


The  second  one  gives  the  relation 

,  f  t  *  fl 


fi  +  ft 


M 


(14.28) 
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The  equation  of  the  Interface  may  be  obtained  from  (10.8) : 


Since  a,  ■  -a0  ,  there  la  a  discontinuity  In  i*  (and  w  if  the  motion 
la  three-dimensional)  as  one  crosses  the  interface. 

The  special  choices  of  cp(x,z)  made  in  14yj  may,  of  courss, 
also  be  made  here.  In  particular,  one  may  make  progressive  and 
standing  waves.  If  one  forms  two-dimensional  progressive  waves  at 
the  interface,  one  finds  for  the  velocity 


yL  -  _fi _ i  S 


(14.29) 


If  one  assumes  the  fluids  bounded  above  and  below  by  planes 
y  *  h*  and  y  ■  -h,,  respectively,  a  similar  calculation  shows 


<Tl  - 


ft  -  h. _ 

fX  CO it>  mVij  -V  fi  CotJi  K,, 


(14.30) 


It  is  clear  from  (14.281  and  (14.30)  that  these  solutions  exist 
only  if  f2  <  f  x.  The  case  will  be  discussed  later. 

A  more  complicated  problem  of  this  type  is  the  following 
(cf.  Lamb  [  1932,  §  23l] ,  Greenhill  £ iO07j ) -  Suppose  there  is  a 
solid  horizontal  bottom  at  y  »  -h,  an  interface  at  y  -  -d  and  a  free 
surface  at  y  *  0.  Then,  in  addition  to  (14.27)  at  y  -  -d, 
and  g  must  satisfy 

V°  r-°  •  K'-°  ai  '~L 

If  one  seeks  solutions  of  the  form 
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2  (^2  tOS»i  S  •  ^  k  \ 

4>1  -  (*-,  cask  -ft,)  tp  ( x,  ]_)  c~os(<rt-+  r)  t 

substitution  In  the  varirjs  boundary  conditions  yields  the  fol¬ 
lowing  relation  between  CT  and  nr 

£  dtk  r»t*  •>  ^  k  -  d )  +  pj_ 

-  — P^fcot’k  Md  -+  c  ofh  r/  4r  Jj  +  —0,  (14  31) 

If  fv  <  f,  ,  one  may  establish  that  there  exist  two  positive 
solutions  for  t"1'  for  a  given  m,  so  that  two  possible  fre¬ 
quencies  are  pcssib’e  for  a  given  wave  pattern  If  the  bottom 
fluid  is  taken  Infinitely  deep,  one  replaces  c  at*  *y  by  1 

in  (14  31)  and  the  two  solutions  simplify  to 


<T,  —  ^  tA 


*1 


f‘>  ~  f  l 
P,  ccfk  w  d 


(14,32) 


The  first  solution,  <Tt  ,  is  the  same  as  would  be  obtained  if 


the  two  fluids  were  identical  (and  there  is  no  discontinuity 
in  u  and  w  at  the  interface);  the  second,  cx  ,  is  interpreted 
below,,  The  inequality  o\v  <  (T  ^  holds  in  general,  and 

one  may  establish 

V  -  v 

~  ( 14 o  33) 


If  one  computes  the  ratio  of  the  amplitude  of  the  disturbance 
at  the  interface  to  that  at  the  free  surface,  one  finds,  no  matter 
whether  h  is  finite  or  not, 

4* o 5 h  ^ a -3.^  S  i  ^ k  *  (14, 34* 

An  examination  of  the  roots  of  (14.31*  shows  that  the  ratio  (14  34) 


is  negative  for  the  smaller  of  the  two  roots  and  positive  for  the 
larger.  Thus,  In  the  solution  associated  with  the  smaller  roots, 
a  maximum  of  the  disturbance  at  the  Interface  is  associated  with  a 
minimum  of  that  at  the  free  surface,  and  vice  versa.  On  the  other 
hand,  with  the  larger  root  the  maxima  and  minima  go  together.  For 
the  values  given  In  (14.32),  the  ratio  becomes 


- 1*  A 


ti 


*  f,  'ft  '  (14.35) 

respectively.  We  note  tnat,  although  the  first  ratio  is  <  1,  the 
second  is  in  absolute  value  >  1  if  J\  (/•*  e*c }  >  ft  >  ,  a  condition 

satisfied  if  ^  Is  only  slightly  greater  than  f,  .  In  fact,  the 
ratio  may  become  very  large. 

For  a  given  wave  length  and  amplitude  of  the  wave  at  the  free 
surface  one  may  also  compare  the  amplitudes  of  the  two  different 
modes  of  motion  at  the  interface  If  A*  is  the  amplitude  associated 
with  the  frequency  C; ,  then  for  the  case  h  «•  <*>  one  finds 


A, 


ft.  *■  f  t  i  —  f a ^ U  (v.J 
which  may  be  either  less  than  or  greater  than  1. 

It  is  of  some  interest  to  examine  somewhat  further  the  solution 
associated  with  the  smaller  root  cr^  of  (14.31).  Then,  since 
a^/b^  *  gm /<r  ,  the  inequality  (14.33)  implies  that  there  exists 
an  h0  with  0  <  h0  <  d  such  that 


%  ** 


tcv. k  ^  J  <  1- 


and  that 


<p,  -  yf*F-  c 


. . v  H'41*1  i»j-’i>i44 W  If  <*W^**^ 


cos l  c*s(crt  +  r)  . 


Ur 


Thus  the  part  of  the  top  fluid  between  y  -  0  and  y  -  -h0  behaves 
as  if  there  were  a  solid  boundary  at  y  »  -h0 ;  and,  of  course,  the 
fluid  between  y  «~h0  and  y  «-h  as  if  it  were  between  solid  boundaries 
If  one  has  selected  solutions  for  cf  which  can  be  combined  to 
form  a  progressive  plane  wave,  then  one  may  conclude  that  the  velo¬ 
city  C£  •  U~: j/m  associated  with  this  mode  of  motion  has  an  upper 
bound: 


In  fact,  when  h  ■  «•,  one  may  verify  immediately  from  (14.32)  that 


J- 


V  e- 


Tt  " 


<■  A  d 
/  ' 


\  "  i  * 


j  1  «-/d  «  l 

Thus  for  h  -  *  a  progressive  wave  travelling  faster  than  c2max  will 
consist  of  only  the  one  mode  of  motion,  i  e  the  one  associated  with 
(T  j •  If  c  <  c2max'  there  may  be  two  modes  of  motion  excited. 

This  fact  is  associated  with  the  phenomenon  of  "dead-water"  resis¬ 
tance  of  ships  (see  Lamb  'l91ta],  Ekman  [1904J  ,  Sretenskli  [  1934]  ) . 

For  superposed  fluids  one  i.-ay  also  find  solutions  analogous  to 
(14.20)  and  (14.21).  Let  us  suppose  that  the  first  (upper)  fluid 
flows  to  the  left  with  mean  velocity  and  the  second  with  mean 
velocity  .  We  wish  to  find  the  possible  steady  periodic  profiles 
of  the  Interface,  assuming  as  usual  that  the  disturbance  is  small. 
The  complex  velocity  potential  for  each  fluid  J s  taken  in  the  form 


hh'i--c4  +  U$'j  ,  FJj)  =  -e.p  +  U}). 


(14  Jo) 

The  conditions  to  be  satisfied  at  the  mean  common  boundary,  y  «  0, 


are : 
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lI'"  ^  -  ci'i  U  j 


c*’  fte  ^  j  -  ^  ;4?  *  *  ci  \  }  . 


(14.37) 


If  each  fluid  extends  Infinitely  far  vertically,  then 


f,  =  £ 


j  fL  -  £ 


give  a  steady-state  solution  if 


al//c.l  "  “a2^C2 


h*  rr 


(fi-Pij 


P,  “► 


~  > 


(14.38) 


where  a2  Is  the  complex  conjugate  of  a2  •  If  the  upper  fluid  Is 
bounded  by  y  ■  h^  and  the  lower  by  y  ■  -h^,  then  the  solution  Is 


f,  -a,  cos  *  K  -  Qi  tcs  “^0, 


where,  letting  ak  -  +  1  /3  It  -  1,  2, 

— '  Si  nil  IV'  In, .  ”  -  ” “  iv  In,  /—  loih  tr,  L  5  0—  Coil i  **  h : 

C,  ’  ;  c,  'tv 


(14.39) 


( h  -  fi) 


C,L  cotti  ■+  Pi  Cj1 


In  either  case  the  equation  of  the  Interface  is  given  by 


H  o' 


Sretenskii  [l952b]  has  considered  a  three-dimensional  analogue 


of  the  above  problem  in  which  the  direction  of  flow  of  one  of  the 


. . . . in . is . in . . . . . m . . IS 


. f! irraw — irmr—aa  * « 
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fluids  makes  an  angle  &  with  that  of  the  other.  Thus,  take 
velocity  potentials  of  the  following  form: 

c|?j_  ( * 1  ft ,  ^  -  -  C  L  (  x  t  o  i  *  J  i  ^  ^  J  +  *  ■  j  i  J )  » 

v  (14.40) 

<f>,  («,  ^  =  -  CjV  -V  <fl  f  *  t'3V 

The  following  are  the  boundary  conditions  at  the  interface  y(x,z) 


for  small  disturbance: 


(?L  (*•<>•])  +  )-  0,  <f1?  *+  ci  7*  =0 

^  (f<  ~  ^a.  ‘  *Ji  r  fl  x  fiC2.L^i*  Siw«P], 


(14.41) 


For  a  solution  in  the  form 


4\  -  (\<  e 


-  'V  Tcos(4,k 

A  * j  ^ 


the  following  relations  mu3t  hold 


Ai  _  _£l  CqS»U  -t  Jr 

/i.  ]_  ^ 


fl  +  <»sJ  +*ls„J)1'  =  S’"  (ft  -fl)- 


(14.42) 


These  reduce  to  (14.38)  for  -5  »  G,  =*  m^  as  they  should.  The 


equation  for  the  interface  is 


3  =  -*1  ~~  S"»(V  +  M 


(14.4.3) 


Srefcenskli  studies  the  properties  of  the  solution  in  more  detail. 

As  a  further  extension  of  the  preceding  cases  one  ma;y  consider 
a  time -dependent  disturbance  at  the  interface  between  two  fluids 
flowing  at  different  velocities.  This  will  be  treated  in  the 
section  on  stability  of  motion. 
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* 1  natural  general ization  of  the  two-fluid  system  is  the 
t.  -fluid  system  (see  Greenhill  (1887))  and  then  the  heterogeneous 
fluid  with  density  given  as  a  series 

f  (  *  7  },r)  -  f'o('j)  ■+  *  -y/yu  -  £vru  - 

If  one  assumes  a  similar  expansion  fo*’  p  and  expansions  for  u,  v, 
w,  and  starting  with  £  ,  one  may  derive  easily  the  linearized  equa¬ 
tions.  These,  discussion  of  same  periodic  solutions,  and  references 
to  the  literature  may  be  found  in  Lamb  [  1932,  \  235] .  Groen  [1948] 
has  shown  that  the  period  for  simple  harmonic  motion  in  the  linear¬ 
ized  problem  is  a  monotonic  increasing  function  of  the  wave  length 
starting  with  the  minimum  /- ^  ^ j  *  for  *>  40.  Groen 

[1950]  discusses  properties  of  internal  waves  in  an  expository  way 
and  gives  further  references  to  the  more  recent  literature.  For  some 
pertinent  theorems  about  waves  in  heterogeneous  fluids  see  section  32y3 

15‘  Gloup  velocity  and  the  propagation  of  disturbances  and  of 

energy. 

In  the  last  section  we  considered  periodic  waves  at  a  free 
surface  or  Interface.  In  this  section  we  wish  to  consider  waves 
of  a  given  but  fairly  general  initial  form  and  study  the  way  in 
which  aey  propagate.  Although  this  will  entail  writing  down  the 
solution  to  a  particular  initial-value  problem,  this  is  of  only 
incidental  interest,  the  chief  interest  being  in  the  history  of 
the  form  of  the  free  surface  or  interface.  Initial-value  problems 
as  such  be  treated  in  more  detail  later  on.  In  fact,  the 

remarks  below  apply  equally  well  to  other  initial-value  problems, 
for  example,  an  initial  distribution  of  velocity  on  the  surface. 

What  is  essential  is  the  resolution  of  the  subsequer  t  motion  into 
a  set  of  waves  moving  to  the  right  and  of  ones  moving  to  the  left, 
as  in  (15.2). 


. . . .  m  Mini  *w*ri*p "  »♦**""*  "f  *  ■  1  «**»  * 
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The  property  of  the  fluid  and  Its  boundaries  which  Is  moat 
important  for  this  investigation  Is  the  functional  relation  between 
the  frequency  c  and  the  wave  number  k,  The  earlie-  parts  of  this 
chapter  have  shown  that  considerable  variation  Is  passible  in  the 
form  of  T(k).  The  two-fluid  example  with  both  free  surface  and 
interface  gave  a  doubly  valued  function,  A  multiply  valued  function 
could  have  been  obtained  with  more  layers  However,  each  branch, 
or  the  branch.  Is  a  decreasing  function  of  k,  approaching  zero  as 
k  -*■  «.  when  surface  tension  is  taken  Into  account  [see  section  24]  , 
the  form  of  <r(k)  for  large  k  cnanges;  it  then  becomes  an  increasing 

x  /? 

function,  behaving  like  k  '  If  h  is  large  enough,  cr(k)  decreases 

initially,  i.e.  for  k  <  k  ,  reaches  a  minimum  at  k  and  then  in- 

m  m 

creases;  if  h  is  small  enough,  c T  k)  Is  everywhere  increasing.  It 
will  be  convenient  to  extend  the  definition  of  cr(k)  to  negative 
k  by  setting  cr(-k)  -  Q"(k). 

15a  .  The  propagation  of  an  initial  elevation. 

Let  us  suppose  that  at  time  t  -  0  the  free  surface  is  given 
by  y  a  'l(x,0)  and  that,  the  fluid  is  at  rest.  How  does  the  free 
surface  behave  subsequently?  One  may  conveniently  think  of  this 
as  an  initial  humping  up  of  the  fluid  near  one  point,  but  this  is 
not  essential.  We  shall  also  suppose  that  ^(x,0)  Is  sufficiently 
restricted  to  allow  a  Fourier-integral  representation.  In  part  of 
what  follows  we  shall  also  assume  it  to  be  square  integrable,  i.e. 
the  total  available  energy  is  finite^and  on  occasion  that  x*j  is 
square  Integrable.  Let 
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where 


-  ■=  J 


„  i-  , ,  l.  .  C  ;  m  \  -  -L  •}  I  m  2  )  t  K  * 

l"’*  1 1  il  *  *  A*  -I*!)  —  rr  ,  < 


F(h)  >  =  j  *£f*'*i‘ Sf*)J. 

^  r  J  *  L 


We  3hall  call  E(k)  the  spectrum  of  7  (  X,  0  ) .  Note  that  E(-k)  - 
the  complex  conjugate  of  E(k)  [we  change  notation  temporarily  in 
order  to  avoid  conflict  with  the  notation  for  averages  Introduced 
below]  . 

A  formal  solution  for  ^  and  may  be  written  down 

immediately: 


^  f».7.rt  =  -  i.  tr"'  "*11 


—  -  S'*. .  -**}  s-  -f  << 


~'Y/n.C^  &  V*  rt 

•  •  >  1  *  I  - 

=  i;(  -*•  Ji<! 

•7  ( t J  ■=  J  [cawejfe*  -<  SM*  *  ‘of  i  * 


(15.2) 
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=  J  e  '^X£^)cos  C“t  M  r  j|  Eft)[e 
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f  ' 
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y - <7^1  -  i|^f*  »  rt) 
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Here  Y(y)  »  cosh  k(y+h)/slnh  kh  for  a  single  fluid  of  depth  h, 

Y(y)  «  e|k,ysgn  k  for  infinite  depth  s the  peculiar  modification  of 
Y  for  h  ■*  uj  la  necessary  for  k  <  0) »  However,  more  general  situa¬ 
tions  are  allowable  In  which,  for  example,  ^(x,  H  describes  an 
interface.  The  choice  of  an  expression  for  <p  has  been  based  upon 
the  kinematic  boundary  condition  $  (x,o in  order  not  to 
exclude  the  possibility  of  surface  tension.  For  simplicity  we  alio 
restrict  ourselves  to  single-valued  qt's.  For  more  complicated 
problems,  such  as  the  two-fluid  problem  with  both  free  surface  and 

1 
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interface  discussed  in  §>145,  the  freedom  to  fix  both  ^(x,  0)  and 
*|2(x  ,  0  )  independently  requires  the  determination  of  two  spectra 
for  each  surface  with  relations  between  them  set  by  (14.34).  The 
remarks  below  will  still  apply  to  motion  resulting  from  each  spec¬ 
trum  separately.  Finally,  we  note  that  a  statement  concerning  spe¬ 
cific  conditions  to  be  satisfied  by  fj(\,C)  for  the  case  of  a 
single  free  surface  may  be  found  in  a  paper  by  Kampe  de  F&rlet  and 
Kotik  [1953]  . 


It  is  c)rar  from  (15.2)  that  one  may  express  v^(x,t)  as  a  sum 
of  two  functions,  one,  say  *^)R(x,t),  representing  a  superposition 
of  waves  moving  to  the  right,  the  other,  riL  ,  waves  moving  to 
the  left.  We  consider  only  R  since  similar  remarks  apply  to  ^ 

•  r  /i,  t 

with  x  replaced  by  -x.  The  spectrum  of  R  Is  given  by  jcln)e 
so  that  clearly  T(k)  plays  an  Important  role  in  the  change  of 
shape  of  17  R.  Since  each  harmonic  component  in  ^R  is  moving  to 
the  right  with  velocity  T(k)/k,  and  since  this  is  not  a  constant 
in  the  cases  we  have  been  considering,  the  different  components  will 
move  with  dlff*  rent  velocities  and  we  shall  expect  ^  R  to  change  its 
shape  with  time^ even  though  moving  as  a  whole  to  the  right. 

In  order  to  get  some  idea  of  the  overall  motion  it  is  reasonable 
to  try  to  compute  an  average  position  of  *jR(x,t)  and  find  how 
this  moves.  One  must  first  decide  how  to  define  the  average  position. 
One  possibility,  which,  as  we  shall  see  presently,  is  unsatisfactory 
is  to  use  itself  as  the  weighting  function,  i.e.  to  define 

J.W  =  /  i 


when  this  exists.  An  easy  computation  shows  that 


i.e.  the  average  motion  Is,  on  this  definition,  independent  of  the 
form  of  ~('<)  except  near  k  ■  0.  For  deep-water  gravity  waves 
7'fo)  m  m;  f  >r  depth  h,  1/  gh  ,  the  maximum  velocity  (see 

equation  (14.15,.).  In  conformity  with  the  above  one  may  define  the 
"spread"  of  the  hump  to  te 

/  40 

j  [*  ~  A,  ft)]'  %{**>*« 
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A  computation  shows  that  this  remains  constant  In  tine,  when  It 
exists.  This  definition  of  average  is  unsatisfactory,  as  could  have 
been  expected  Inasmuch  as  the  weighting  function  can  become  negative. 
We  note  In  passing  that 

_  ** 

an  expression  of  conservation  of  mass. 

Another  possible  weighting  function  without  this  shortcoming, 
but  still  allowing  ease  of  computation,  is  *J £( X  ,  t  ) .  We  note 
first  that 

S  6  V  =  J  7r  Mci*  -  i,7rj^  e  fc)  l  h 

•  9C  **■ 

Let  us  define  two  averages,  one  for  functions  of  x: 
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and  one  for  functions  of  k: 


Then,  assuming  that  the  various  quantities  in  question  exist,  cne 
finds,  using  well  known  theorems  on  Fourier  transform?  fsee,  e.g., 

r  1 

Bochner  and  Chandrasekharan,  1949,  ch.  IV,  S*  2  , 


Xe  ft)  -  *R 


C'j  +  ^  t 


(15.3) 


Thus,  on  this  definition  the  average  position  of  *1 R  moves  to 

the  right  with  constant  velocity  and  the  hump  spreads  according 

2 

to  a  quadratic  law.  We  note  that  the  coefficient  of  t  la  positive 
except  if  t'  is  a  constant,  when  it  vanishes.  It  may  become  in¬ 
finite,  and,  in  fact,  dees  so  for  infinitely  deep  water  if  the 
gravest  modes  are  present,  i.e..  If  dx  /  0.  The  coefficient 

of  t  vanishes  if  t'  Is  constant  or  if  ufc^tVC]  Is  constant; 

the  latter  will  occur  if  X,  C  )  is  either  symmetric  or  antisym¬ 
metric  about  some  print  xo,  but  this  does  not  exhaust  all  possibili¬ 
ties.  The  sign  of  this  term  does  not  seem  to  be  determined,  so  that 
the  spread  of  the  hump  may  conceivably  decrease  before  starting  to 
increase . 

Investigations  of  the  motion  of  the  average  position  of  the 
hump  and  of  its  spread  give  only  a  rather  crude  picture  of  Its 
behavior.  By  other  methods  outlined  below  one  may  obtain  further 
insight  into  tne  motion. 

We  begin  by  applying  the  analysis  of  the  average  motion  to  that 
part  of  h R  resulting  from  only  a  narrow  band  in  its  spectrum.  Let 
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We  shall  call  this  a  wave  packet.  The  average  position  satisfies 

^"9  ♦'e.SJ  “  *  T  J 

where  a*7  (ft ,  i)  la  now  the  average  of  ?'(k)  over  the  narrow  band 

[  ,  he  ^fc]  *  The  narrower  the  band,  the  closer  T'(P(.,y.j  la 

to  <T'Y^c)»  assuming  the  latter  continuous.  A3  a  limiting  case  we 

shall  say  that  the  wave  packet  resulting  from  an  Infinitesimal  band 

about  k  moves  with  velocity  •  It  13  customary  to  call  crYrt) 

o 

the  group  velocity.  This  Is  the  tame  as  the  phase  velocity  vfK/k 
only  if  rzak  .  A  wave  packet  will  spread  with  passage  of  time 
unless  the  two  velocities  are  equal,  for  (15.4)  Is  applicable  to 
the  wave  packet  with  the  restricted  definition  of  average.  As  might 
be  expected,  the  smaller  the  width  of  the  band,  the  smaller  the 
coefficient  of  t2  and  the  smaller  the  rate  of  growth.  However,  as 
we  shall  see  below,  the  initial  spread  may  be  wide  for  a  narrow  band. 

The  wave  packet  (15.5)  may  also  be  investigated  by  a  different 
method.  Let  us  expand  T(k)  in  the  first  few  terms  of  a  Taylor 


series  about  kQ  : 


(15.6) 


We  may  then  write 


iKfcfe  |  i* 

■+  L  lexf  —  j  (15.1 

,  IP*  *  «" M  (a  .  0-tfUh  A.,*) 

1  * 


,  one  finds 


Using  the  inequality  i  a  -  ■  >  -  M 

io|  <  -  t  i'  no*  i  l'k)\  '*>-*  |  r  '  *'l  (15.8) 

^  \h.ha\<l  •  k  U.Mt 

The  remainder  can  thus  be  made  smaj.1  by  taking  £  or  t  small 
enough.  However,  once  £  is  fixed,  R  will  eventually  become 
large  as  t  Increases.  Let  us  suppose,  however,  that  t  and  L 
are  small  enough  so  that  the  first  term  determines  the  main 
features  of  the  motion.  The  first  factor  represents  a  periodic 
wave  of  wave  number  k0  moving  with  its  phase  velocity  CT(ko)/k0. 

The  second  factor,  determining  the  amplitude  of  the  first,  rep¬ 
resents  a  profile  being  translated  to  the  right  with  velocity 
cr^ko).  Thus  one  may  say  that  the  gross  outline  of  the  surface 
is  moving  to  the  right  with  the  group  velocity.  One  may  see 
this  more  clearly  If  one  assumes  £  small  enough  so  that  we 
may  take  E(k)  as  constant  ovir  tie  band  width.  Then 

siv>  (x-<r'(J<0)t)z_ 

and  ^0,£)  appears  approximately  as  in  Figure  9. 

Here  the  dotted  enveloping  curves  represent  +  j  M  and  move  to  the 
right  with  velocity  T/(kQ) ,  whereas  the  inscribed  solid  curves  rep¬ 
resent  the  first  factor  and  move  to  the  right  with  phase  velocity 
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<f(ko)/k0.  The  whole  moves  as  a  fixed  pattern  only  if  the  two 
velocities  are  equal.  Otherwise,  assuming  ;y'(k  )  <  0“(k  )/k  , 

the  inscribed  curves  will  progress  through  the  wave  packet,  gradually 
disappearing  at  the  right.  For  a  very  narrow  band  the  packet  will 
spread  wide  before  Its  first  zero  on  either  side  of  the  maximum. 

A  disadvantage  of  this  last  analysis  is  that  it  becomes  les3 
and  less  accurate  as  t  becomes  ^arge.  However,  there  exists  another 
approximation  to  ^R(x,t)  for  larga  values  of  t  which  helps  to  com¬ 
plete  the  picture.  This  ultimate  behavior  of  ^R  can  to  some  extent 
be  predicted  from  the  analysis  of  the  average  motion  of  a  wave  band. 

If  we  think  of  ^  R  as  made  up  of  the  contributions  from  a  number  of 
narrow  wave  bands,  we  know  that  each  contribution  is  moving  with  the 
average  group  velocity  of  the  band.  Thus  after  some  time  we  shall 
expect  that  these  various  contributions  will  have  separated  from  one 
another,  with  the  bands  about  the  gravest  modes,  which  travel  fastest, 
having  progressed  the  furthest.  This  prediction  will  be  confirmed. 

Wuat  is  needed  for  this  final  approximation  is  an  asymptotic 
expansion  for  large  t.  It  is  convenient  to  express  ij  R  in  the 
slightly  altered  form 


In  W) 


(15.9) 


and  to  consider  it  as  depending  upon  the  two  parameters  x/t  and  t 
Then  for  each  value  of  x/t  we  3hall  give  an  expansion  for  large  values 
of  t.  For  a  derivation  of  the  expansi  .  we  refer  to  Stoker  £l9S7,$6.Q 
or  Erdelyl  [  1956,  §  2.9J  . 
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Let  the  functions  kJx/t),  r  «  1,  2,  . n,  be  defined  by 

P 


VUK)**/*.  (15.10) 

l.e.  we  allow  the  possibility  of  several  roots.  In  the  situation 
of  interest  to  us  there  will  be  either  one  or  two  roots,  or  none. 
The  asymptotic  expression  for  R  is  then  given  by 

>,  (15.11) 
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where  the  first  summation  is  over  all  values  of  r  for  which 

0-"  (k, ,)  /  0  and  the  second  over  all  k  for  which  t  “  (k  )  »  0  but 
jr*  r  i 

<T";  (k  )  ft  0;  further  terms  would  be  necessary  for  values  of  r  for 

which  both  vanish  but  this  will  not  occur  in  our  examples.  If  some 

1 

k  «  Q,  then  the  corresponding  term  must,  be  multiplied  by  W  .  For 

i’  <- 

a  value  of  x/fc  for  which  no  solution  to  (15,10/  exists,  it  is  easy 
to  show  by  a  change  of  variables  in  (15.9),  say  u  >«  kx/t  -  cr(k) , 
and  integration  by  parts  that  ^_{x,t)  **  0  (t  "* )  . 

Let  us  examine  in  some  detail  the  implications  of  one  term  of 
(15.11),  say  r  «  1,  for  the  motion  of  rj  if  several  terms  are 
present  for  a  given  value  of  x/t  one  must  superpose  the  resultant 
motions . 

If  x/t  la  held  constant  while  t  increases,  then  clearly  one 
must  set  x  -  cr,(k1)t,  l.e.  we  are  examining  v)  R  from  the  standpoint 

of  an  observer  moving  with  group  velocity  (r'(k^).  Since  the  coef¬ 
ficient  of  the  harmonic  term  is  times  a  function  of  which 
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Is  being  held  constant,  the  gross  outline  of  *j  R  will  appear  con¬ 
stant  In  form,  but  decreasing  in  amplitude  because  of  {  However, 

lust  as  in  the  analysis  of  (Li  7).  there  is  a  harmonic  ov  wave  number 
ft  moving  through  the  gross  outline  with  phase  velocity  ofk,)  / . 
The  amplitude  of  the  gross  outline  is  proportional  to  E(*(.) ,  bite 
also  depends  now  upon  o'1''  ^  j  in  contrast  to  the  situation  for 
small  t  according  to  (15  7) 

If  the  value  of  */i  is  such  that  -O  ,  then  one  must 

examine  a  term  from  the  second  summat'on  in  (15  11)  It  is  evident 
that  the  Interpretation  is  the  same  except  that  cr  "  occurs  in  place 
of  cr"  and  that  the  amplitude  decreases  more  slowly  because  of  the 
f 1/3 ,  '.Tils  situation  can  happen,  for  example,  in  the  case  of  gravity 
waves  in  wafer  of  depth  h  f  o**  x  - 1  f^T  Then  4,^  syT )  *  ^  > 

cr"(o)  -  O  ,  ar.d  This  ai-so  occurs  for  combined 

gravit  y-capi  i  lary  waves  when  the  curv*»  J  7 .  /? )  has  a  minimum 

The  approximation  ■ 15  ll)  tc  *  will  obviously  be  very  poor 

1  K 

for  a  value  x/t  such  tha'  cr'VVf  j  is  near  tc  zero  for  some  r  unless 
t  is  extremely  large.  It  is  shewn  in  Jeffreys  and  Jeffreys,  Methods 
of  mathematical  physics  (.3rd  ed  Cambridge,  1956,  §  17  09]  how  an 
Airy  function  may  be  used  to  modify  the  relevant  term  in  the  second 
summand  to  give  a  usefu.  asymptotic  expansion  for  r  near  a  zero 
of  O'"  (see  also  Chester  Friedman  and  Ursell,  Proc  Cambridge 
Philos,,  Soc,  5  3  <195 n,  599  611  ) 

TSSm 

If  x/t  is  fixed  at  a  value  for  which  (15  10)  has  no  solution, 
then  for  an  observer  moving  with  this  velocity  the  disturbance  of 
the  surface  is  very  small,  for  it  has  been  dying  out  as  t  .  The 
first  term  of  the  expansion  may,  of  course,  be  computed  as  indicated 
above.  This  situation  will  occur  for  a  disturbance  ir,  water  of 
depth  /,  if  x/t  T  «  It  will  also  occur  when  surface  tension 

is  taken  into  account  for  x/t  < 

. HUmwHM**** **»«*it 1  ll,w '**i*W' 1 1 '  ** *  w  _  ** " J ^ "lirT| ~  'i"1 ^ uI.Tr'  *"’rT " "TrfTfflTH " 11 W1 


The  asymptotic  expansion  (15.11)  may  also  be  used  In  a  different 

fashion.  Let  us  fix  our  attention  upon  one  value  of  x  and  let  t 

increase.  Then  x/t  will  decrease  and  the  value  kj^x/t)  associated 

with  the  point  x  at  a  given  moment  will  also  change;  for  pure  gravity 

waves  It  will  Increase.  The  observer  stationed  at  x  will  then  observe 

waves  of  continually  increasing  wave  number  (decreasing  wave  length) 

moving  by  with  phase  velocities  appropriate  to  their  lengths.  The 

amplitudes  at  a  given  instant  will  depend  upon  the  first  factor.  The 

gross  outline  of  the  waves  will  pass  the  observer  at  the  group  velocity 

appropriate  to  the  wave  number  present  at  the  moment,  and,  of  course, 

-1/2 

the  amplitude  is  decreasing  as  t  '  .  In  the  case  of  a  disturbance  on 
water  of  depth  h,  if  the  observer  is  initially  far  from  the  hump,  then 
even  for  large  enough  values  of  t  for  the  asymptotic  expansion  to  be 
valid  the  value  of  x/t  may  be  greater  than  /gh.  Then  the  observer 
will  see  practically  no  disturbance  until  the  gravest  modes  begin  to 
reach  him.  We  note  again  that  he  must  anticipate  the  arrival  of  a 
given  wave  number  by  its  group  velocity,  not  phase  velocity,  for  it  is 
the  former  which  controls  the  amplitude.  In  the  case  of  combined 
gravity-capillary  waves,  when  t  is  large  enough  one  will  have 
x/t  <  and  the  disturbance  will  be  negligible. 

It  is  also  possible  to  find  an  asymptotic  expansion  for 
r(x,c)  for  x/t  fixed  and  large  x.  It  turns  out  to  be  the  same  as 
(15.11)  with  0  (t“2/3)  replaced  by  ^(x"2/3).  This  expansion 

allows  one,  so  to  speak,  to  take  snapshots  of  the  right-hand  end  of 
j  at  different  instants  of  time.  If  we  fix  t  and  let  x  increase, 
x/t  increases  also  and  k^(x/t)  increases.  Thus  the  wave  length 
decreases  as  one  moves  to  the  fight;  the  observed  amplitude  will 
depend  upon  the  first  factor.  For  gravity  waves  on  water  of  depth  h, 
if  x  is  large  enough,  x/t  >  \iTgFi  and  the  disturbance  will  be  small 
of  order  x“*. 
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Finally,  we  jse  the  asymptotic  expansion  to  Investigate  the 
motion  of  a  particular  phase  of  •'J  R{x,t),  say  a  zero,  for  large  t 

3uch  a  point  will  be  determined  by 


a(x,t)  -  k1x  -  crCk^Jt  m  const.. 


where,  as  usual,  *  k^(x/t);  solving  for  /  gives  x 
may  find  x(t)  from 


x(t).  One 


t  -k',  _  £**•). 


Thus  a  particular  pnase  travels  with  the  phase  velocity  of  the 
harmonic  com;,  snent  associated  with  It  at  the  moment.  ..owever,  If 
the  group  and  phase  velocities  are  different,  it  is  then  moving  at 
a  different  velocity  from  a  point  just  keeping  pace  with  waves  of 
a  given  wave  number.  In  particular,  for  gravity  waves  it  is  moving 
faster,  hence  moves  into  region  of  lower  wave  numLer  and  higher 
velocity  and  is  accelerating.  A  computation  of  x  bears  this  out: 


i 


*  '  TJ 


,T“  1*.  J  ' 


for  this  Is  always  positive  for  gravity  waves.  The  right-hand  side 
is,  of  course,  a  function  of  x  and  t.  For  deep-water  gravity  waves 
the  function  x(t)  may  easily  be  found  from  the  earlier  equation: 


x*  T7  %  -  *<r'fw  *  T, 


,  Ug.  Xv 

*  r  7* 


or 


m-  £v. 


. . . . . 


Hence  x  -  g/2a  and  for  large  t  the  acceleration  is  constant.  If 

the  depth  is  finite,  the  computation  is  no  lo.  ger  simple,  although 

it  is  possible  to  show  that  x(t)  varies  from  x(t)  -  t  /  gh  for  a 

2 

phase  associated  v.  th  k  *  0  to  x(t)  *  At  for  a  phase  associated 
with  very  large  k. 

Figure  10  is  taken  from  a  paper  of  Kelvin's  [l9C ’J  .  and  shows 
the  computed  values  of  *|(x,t)  for  an  initial  displac  -ment  given  by 
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U 


2V‘  0 


'vj 


1/2 


\l  -  (HXl 


),/L] 


and  for  ^  ,  1 ,  |-,  4,  8  (the  units  have  been  chosen  so 

that  g  m  4).  The  description  of  tne  behavior  of  ^  R  { x ,  t )  outlined 

in  the  preceding  paragraphs  can  be  easily  verified  qualitatively 
by  inspection  of  the  successive  snapshots  of  h  R(x, t;).  Green  [1909] 
has  shown  that  if  one  estimates  the  wave  length  at  ary  maximum  as 


double  the  distance  between  the  two  including  zeros,  then  the  posi¬ 


tion  is  very  close  to  that  which  would  be  estimated  by  using  the 

group  velocity  fcf.  Havelock,  191**,  p.  37j  • 

L  -1  1;  also  1958b, 

Figure  11  from  a  report  by  J.  E.  Prins  [_195qJ  shows  measured 
time  histories  taken  at  various  distances  from  tne  center  of  an 


initial  rectangular  hump  of  length  2L  and  height  Q  in  water  of  depth 
h  for  specific  values  shown  in  the  figure.  In  general,  the  features 
of  the  motion  described  above  were  well  verified  by  this  experimental 
investigation. 

We  assemble  here  the  expressions  for  cr(k)  and  kc r'/cr  for 

a  number  of  cases  of  water  waves. 
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1)  Deep-water  gravity  wave  at 


. I 


6)  Waves  at  a  free  surface  of  a  fluid  of  depth  h  with  both  gravity 
and  surface  tension  acting: 


df  k  >  0 , 

In  cases  1)  to  4)  T  is  always  negative),  In  case  5)  it  crosses  the 
k-axls  at  k  -  f j  ^  and  becomes  positive.  In  cases  i) 

to  4)  cr'<  7/k  for  k  >  0.  In  case  5)  7  '<  Or/k  for  0<k<  / g  f  /T  ; 

then  7'  crosses  cr/k  at  the  minimum  of  the  latter  and  thereafter 
remains  larger.  (Note  that  t'  always  passes  through  a  stationary 
value  of  7 /k ,  passing  from  beneath  to  above  in  going  through  a 
minimum,  and  the  reverse  at  a  maximum.)  We  shall  not  discuss  6)  in 
detail.  For  h  >  hc  -  /3T/2  J5  g\  CT/k  has  a  minimum  for  some  kj( 

0  <  kQ  <  /pg/r  *  and  cr'  a  minimum  to  the  left  of  this.  For 

h  <  h,  ,  CT/k  is  an  increasing  function,  starting  at  /  gh  for  k  ■  0, 

and  t'  Is  also  increasing,  rj'  >  CT/k  for  k  >  0,  7 '(o')  -  /gh  . 

Figure  12  shows  graphs  of  a  ,  J / k  and  for  1),  3),  5),  and  6) 

(the  scales  were  cno3eri  for  convenience). 

One  may  also  tf.ce  )  «  2  7?  ,/k  as  the  independent  variable,  and 
then  express  the  ptase  velocity  c  and  group  velocity  U  as  functions 
of  )\  .An  easy  computation  shows  that 

=c'u' 

Tills  equation  has  a  simple  Interpretation  in  the  geometry  of  the  curve 

hr  c(Jl),  as  was  shown  by  Lamb  [  1932,  p.  38?.]  :  For  a  given 


value  of  A  ,  U  la  the  intercept  on  the  vertical  axis  of  the  tangent 
to  the  c ( A  )  curve  at  the  point  ,c(,\  )).  One  value  of  Umay 
correopond  to  more  than  one  value  of  A  ,38,  for  example.  In  the 
case  of  gravity-capillary  waves,  [see  Havelock,  1914,  $11.] 

IS,-?  .  The  propagation  of  energy. 

It  seems  lntultlv  -ly  clear  that  as  long  as  the  right-moving 
part  of  an  initial  hump  keeps  Its  Integrity  the  energy  associated 
with  the  motion  will  in  some  sense  move  with  the  hump.  We  wish  to 
consider  in  what  aenae  this  Is  true.  We  limit  ourselves  In  the 
following  discussion  to  a  single  fluid  of  depth  h,  where  h  may  be¬ 
come  Infinite.  However,  surface  tension  may  act  upon  the  free  sur¬ 
face  . 


We  first  introduce  the  notion  of  energy  density  for  a  given 
value  of  x.  It  will  be  convenient  to  separate  potential,  kinetic  and 
surface  energy.  Let 


V^-  n  fa  -  i  fi?r< . 


-  iT 

be  the  densities  of  potential,  kinetic  and  surface  energies,  respec¬ 
tively,  where  here  <p  Is  the  velocity  potential  corresponding  to  ^  R 
These  functions  may  now  treated  in  the  same  way  as  ^  ^ 
was  In  15  oL  .  We  may  ask  for  the  average  position  of  the  distribu¬ 
tions  of  the  several  densities.  They  are  defined  by 


1 


*  :-0  ~  30 

/.l 

>  rjfii  f  •*  -<0 


(15.13) 


respectively.  Since  all  three  densities  are  non-negative,  one 
avoids  the  difficulty  met  with  in  defining  the  average  position  of 
*7  R-  In  fact,  it  is  obvious  that  the  definitions  of  xR  and  xv 
coincide,  so  that  the  conclusions  com ernlng  xR  can  be  applied  immedi¬ 


ately  to  xv(t).  In  particular, 

*v(t)  *  *v(°)  4  T  ^  '  (15.14) 

Consider  now  xT(t).  First  we  note  that,  from  Green's  Theoi’em, 

[  T =  if  [  **  + 


~  «o 


V  — —  tA-*  V 


From  the  assumed  square-lntegrability  of  7  R,  the  limit  vanishes. 
Use  of  the  identity  x($>^  +  $  y)  »  (x$  )x£x  +  (x$)y$y  -  $  x 


and  Green's  Theorem  gives 

/■  °° 


•  00 


x  TU.^Jx  a  J  f  j  x  #Cx,»,t)^{x,o,tijx  -  JjrJff 

o 

+  lx*  *?)  j]i), 

x!~»  06  ”  ^ 

where  again  the  last  two  Units  vanish.  A  similar  computation  shows 
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Collecting  these  results  we  have 


rrcx.fWx  -  if  j 

-  *c 

(,“*  tmvu  = 


>•  w 

Vro+u*  - {ff  s&xaM >,o,tw*+if I ) 

J  V  ,  •  rjO  ”  Oft  ~  U 


(15.15) 


Since  from  (15.2),  .  >, 

,  r*V„\  J  i, 

7rI*^  r  iLff* 


(15.16) 


one  finds  easily 


^v(x,0it)-  f  )  <r^)EIK)i 

Y  *-06 


- ;  -  *■+) 


One  may  now  apply  again,  as  in  section  15<X  ,  theorems  on  Fourier 
transforms  to  obtain 


mill  I  Hill  l 


CO  •*  V  . 

(  ixs  -■*?  ^  f  j  co*i»  *k , 

«  0O 

'M  °**  *L 

j  X  ^/,f) jx  -  ■—TTfj"  *  ^ ''Vic 7|':-i-  cstii^i  Jft  +^iTftj  fa) cr^i ■—  <o>^#k  'U , 
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f°V  T(*,tu*  =  fV T>  •)«>  +  i*t t  j  tcW^c'f'toU* *k  J*  (15-17) 

-  oO  -  **  ** 

+  7  ^  j  t  ><■  t  *fk)  r'W;  ~  c L 

If  one  uses  the  definition  introduced  earlier  for  average  of  a 


function  of  k,  one  now  finds 

x  m  --  yTh)  -  i 

4  k 


(15.18) 


and  a  rather  unwieldy  expression  for  [x  -  x^t)]2  ,  similar  in  char¬ 
acter  to  (15.4).  We  note  that  if  we  are  dealing  with  pure  v  avity 
waves,  so  that  or2  -  gk  tanh  kh,  then  formulas  (15.17)  simp.  j 
considerably  and  become  identical  with  those  for  V.  In  this  case 
the  potential  and  kinetic  energies  are  equal  and  propagate  with  the 
same  velocities. 

We  may  now  carry  out  similar  calculations  for  S(x,t).  The 
corresponding  formulas  follow 

f  JM<!>  =  Jttt 

-  oo  -  oc 

OKJ  CO  »  M, 

^  xM*M*  =  J  i  r'lmmiMM, 

_  dCi  -o£ 


.  AJ  Oj  -K> 

J  X^ySMi*  =  *%)E(k\ 

*h  -  oo  «oe 

^  i t: t r  ^ v,-l  f we*?#) m , 


(15.19) 
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and 


=  Xslt)  +  t 


( 15 *20] 


and  again  a  formula  for  [x  -  xs(t)]  2  similar  In  character  to  (15.4 
One  should  note  that  the  total  potential,  kinetic  and  surface 
energies  associated  with  *JR(x,t)  each  remain  constant  in  time.  If 


). 


T  /  0,  then  the  mean  positions  cf  the  three  energy  densities  propa¬ 


gate  with  different  velocities,  each  velocity  being  an  average,  in 
some  sense,  of  .  If  one  considers  a  wave  packet  (15.5),  then 
as  the  width  Zg  of  the  band  of  wave  numbers  approaches  zero  the 
velocity  of  propagation  of  the  Individual  energy  densities  will 
each  approach  a,/(kQi),  the  group  velocity. 

Consider  now  the  total  energy  density. 


Making  use  of  the  form  of  <7( k), 

3-ViO  -  Ifh  +  Tk'/f  )U'\  -U, 

one  flnd3 

j  A*  "PjjS  +f  foK  ^  4 

-  am 

-  i  7T  f  f  +  T4*J  £lt)E'r* )J  *, 
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At  any  instant  t  half  of  the  total  energy  la  kinetic  energy  and  the 

other  half  la  divided  between  potential  and  surface  energy. 

There  la  another  way  of  considering  the  energy  transported  by 

surface  waves  which,  at  first  glance,  is  different  from  the  preceding 

treatment.  Consider  a  fixed  plane  x  «  constant.  Theij  from  the 

results  in  §  8  one  may  compute  the  rate  at  which  energy  is  being 

transported  through  this  plane,  the  so-called  energy  flux.  Let  us 

denote  it  by  ?(x,t).  After  appropriate  linearization,  formula  (8.10) 
gives 


t' 


^ 1 5  *  22 J 


The  expression  for  the  flux  has  an  advantage  over  the  expres¬ 
sions  for  mean  positions  considered  above  in  that  no  strong  restric¬ 
tions  upon  are  required  for  it  to  exist.  In  fact.  It  can  be 
computed  for  a  single  harmonic  wave 


1 


=.  /\  sin 


( 1 5  •  24r  ) 


With 


d>  - -a- ?  i!^l 

Si*k  ^k 


one  finds  by  a  straight  forward  calculation 


y*t)~AXT-hr  ttf’fkj'-irti+tff  ^f,+  s'" 


Averaging  over  a  wave  length  (or  over  a  period,  it  makes  no  difference 
which),  one  finds 


(15,25) 
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J  /\X  (of  +  T  Al)crYA). 


Thus  the  group  velocity  enter®  again  In  connection  with  energy 
propagation,  even  though  no  "group"  is  present.  The  energy  density1 
and  average  energy  per  wave  length  for  (15.24)  are 

£i*,t) ;  a1/  j  f? 


(15.26) 


£,.  -  i  (jf  + 


If  one  is  dealing  with  a  composite  wave,  averaging  over  a  wave 
length  Is  possible  only  if  the  resulting  wave  is  periodic.  However, 
even  without  this  restriction,  one  may  compute  both  the  average  flux 
and  average  energy  per  unit  length  from 

X  -  I'Z 

(_  (15.27) 

c  :  J.  | 
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Then  if  a  composite  wave  propaf  ting  to  the  right  is  given  by 


oO 


with 


(x,t)  =  Z-  cJe  1  ] 
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(15.28) 
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wh ere  £ j  =  f.j  *  5  if  « [  '  *  j)  4.^  ■  -  *  j  ,  X;  s  s'ftj  1  =  -  7j  ,  one  finds 
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and 
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In  order  to  obtain  these  relatively  simple  formulas  in  which  the 


contributions  from  the  Individual  harmonics  are  isolated.  It  la  essen¬ 
tial  that  the  averages  be  taken.  Otherwise,  for  £  (x,t)  or  $(x,t) 
one  obtains  a  complicated  double  summation,  and  the  role  of  the  group 
velocity  is  not  apparent. 

A  similar  analysis  may  be  carried  through  for  the  right-moving 
Initial  hump  (15.16).  However,  an  average  of  either  3»  or  ££  computed 
according  to  (15.27)  would  vanish.  Instead  we  take  the  total  flux 
and  total  energy,  respectively; 

oo 

7«.,„  =  -L  ^ (15.31 ) 

The  resulting  formulas  are  analogous  to  (15.29)  and  (15.30); 
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If  the  last  result  la  applied  to  a  narrow  wave  band,  such  as  (15.5), 
then  one  flrd#  the  limiting  relationship 

-h±l  = 

€<%° 

In  the  first  method  of  treating  the  propagation  of  energy, 
l,e,  in  terms  of  the  motion  of  the  mean  position  of  the  energy 
density,  It  was  not  surprizing  that  cr'  should  appear,  for  it  is  s 

familiar  property  of  Fourier  transforms  that  taking  the  derivative 
of  the  transform  is  associated  with  multiplying  the  function  by  the 
variable.  Thus ,  if 

$  fh)  “  ^  i (V>  a  e 

then 

C  °°  .  r  *  'hx  j 

&  -  J  *  X  e  ’ 

Jvt 

In  ti.e  cases  considered  above  the  transform  contained  S  as  a 

factor,  and  the  derivative  contained  cr/t  in  one  summand.  However, 

the  appearance  of  Q"'  In  the  formulas  'or  or  3j.otai  seems 

in  some  ways  coincidental:  One  makes  a  calculation,  and  after 

gathering  and  manipulating  terms  discovers  that  a  certain  combination 
contest!  ns 

of  them  indeed  *  v'  That  this  is  not  really  coincidence  is 
Indicated  by  the  following  theorem  for  the  case  (15.16): 

SO 

A  [  *  £Cx,i)kx  -  fftoJt’  (15.53) 
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It  may  be  proved  as  follows.  From  the  definition  of 
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Integrating  the  second  and  third,  tenrm  by  parts  and  taking  account 


of  the  assumed  behavior  of 


and 


* 


at  <f  one  finda 
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*1  V 


hU 


Since  d  + 

*  XX 


first  integral  as 


0,  one  may  express  the  third  summand  In  the 


Hence  the  first  Integral  may  be  written 

t  1  ^ 

which  vanishes,  since  the  term  in  brackets  is  Just  the  dynamical 
boundary  condition  at  the  free  surface.  The  second  integral  above 
is  Just  k7Lnt.«-i  t  so  that  (15.33)  is  proved. 

bwvaJL  1  '  r 

A  similar  line  of  reasoning  allows  one  to  establish  the  follow 
lng  relation  between  £  and 

<5 _  _  o>  1( *,  ^  (15.34) 
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